# Memetic Salp Swarm Algorithm-Based Frequency Regulation for Power System with Renewable Energy Integration

Authors: Fang Zeng; Hongchun Shu
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

As the penetration of renewable energy to power grid increases gradually, to ensure the safety and stable operation of power system, it is necessary for renewable energy to participate in the secondary frequency regulation of power system. Therefore, this paper proposes an optimal control model of renewable energy participating in the secondary frequency regulation to solve the dynamic power distribution problem. Besides, memetic salp swarm algorithm (MSSA) is used to solve this complex nonlinear optimization problem, so as to quickly obtain high-quality power distribution schemes under different power perturbations and maximize the dynamic response regulation performance of the entire regional power grid. Finally, based on the improved IEEE standard two-area model, the established model is verified and the performance of the applied algorithm is tested by comparing the traditional engineering allocation method and other intelligent optimization algorithms.


---

## Body

## 1. Introduction

In recent years, to cope with climate change and sustainable energy development, the penetration of renewable energy connected to power grid has increased rapidly [1–3]. Different from traditional hydrothermal power, wind power and photovoltaic (PV) power with relatively mature technologies are greatly affected by meteorological conditions, so their power fluctuations are highly random [4–6]. As the regional power grid usually cannot fully absorb renewable energy such as wind or solar energy, it is easy to disconnect the wind farm or PV station from the power grid, which also greatly increases the pressure of frequency regulation of the system [7–10]. In this situation, the traditional configuration of hydropower plant and thermal power station as the main frequency regulation resources has been difficult to meet the high quality of the system’s dynamic frequency regulation needs [11–13].Therefore, the development of new high-quality frequency regulation resources has become one of the main means to relieve the pressure of regional power grid [14–16]. Compared with traditional hydropower plant and thermal power station, wind farm and PV station have faster response speed as well as climbing speed and rapid power fluctuations can be balanced more quickly [17–19]. Therefore, for regional power grids with high renewable energy penetration, especially in low-load areas, and wind farm, PV station can be used on power point control method to control it in below the operation condition of maximum power point, with a certain reserve capacity to participate in the secondary frequency regulation [6, 20, 21], called automatic generation control (AGC) [22–24].In general, AGC mainly consists of two parts [25–27]: (1) based on the real-time acquisition of frequency and power deviation of the tie line, a centralized controller, such as proportional-integral (PI) controller, is used to acquire the approximate actual power fluctuation of the system, and then the total power instruction of the regional power network is issued; (2) according to the power allocation algorithms, the total power instruction is assigned to each unit participating in the frequency regulation [28–30]. Literature [31] proposes a sliding mode controller for multisource AGC system using teaching and learning-based optimization algorithm. Literature [32] proposes a lifelong learning-based complementary generation control of power grids with renewable energy sources. Literature [33] proposes AGC of a multiarea multisource hydrothermal power system. However, although the abovementioned literature realizes the control of multisource participating in AGC, the modeling is relatively simple, considering only the climbing response characteristics of the units and not other dynamic response characteristics of different frequency regulation resources, which will affect the overall control effect of the system and easily make the system deviate greatly from the optimal operating point. Therefore, this paper proposes a multisource optimal collaborative control method for wind and solar renewable energy based on its dynamic response characteristics to participate in secondary frequency regulation.In essence, AGC optimal control is a complex nonlinear optimization problem [34–36]. In practice, power distribution is often not optimized but arranged according to adjustable capacity ratio and climbing speed, which cannot meet the optimal control requirements of the system [37–40]. On the other hand, traditional mathematical optimization methods (such as interior point method [41]), although fast in solving problems, have poor global searching ability and are prone to fall into local optimal solutions. In comparison, genetic algorithm (GA) [42–44] and other metaheuristic algorithms [45–47] have higher application flexibility and better global search capability, but their solving speed is slow and cannot meet the needs of AGC online control for large-scale regional power grids [48–50]. Therefore, MSSA with a faster convergence speed is used to solve the problem. Compared with the original salp swarm algorithm (SSA), memetic salp swarm algorithm (MSSA) employs multiple independent slap chains to simultaneously implement the exploration and exploitation [51]. Besides, MSSA also has low dependence on the mathematical model. To verify the validity of the proposed method, this paper used the improved IEEE standard two-area model for simulation test and analysis.The remaining of this paper is organized as follows: Section2 develops the optimal control model for automatic generation control. In Section 3, MSSA is described. Comprehensive case studies are undertaken in Section 4. Section 5 summarizes the main contributions of the paper.


## 2. Optimal Control Model for Automatic Generation Control

### 2.1. Control Framework

Generally, AGC has the following three control modes: (a) flat frequency control, (b) flat tie-line control, and (c) tie-line bias control. Also, tie-line bias control was used in this paper. The two-region interconnection power network is depicted in Figure1, in which the AGC control process of each region includes two parts: controller control and optimal power distribution. A PI controller is adopted to convert the power deviation and frequency deviation of the real-time receiving tie line into regional control deviation as input and output to the real-time adjusting power of the regional power grid ΔP. Then, ΔP is assigned to each AGC unit through the power allocation algorithm. Different from the traditional power system with hydropower plant and thermal power station as the main reserve capacity, wind farm and PV station no longer need to be disconnected from the grid but participate in the power AGC regulation of the system.Figure 1
Illustration of coordinated control of multisource for AGC under two-area framework.


### 2.2. Dynamic Response Model of Units

The establishment of dynamic response model of the units is mainly in order to more accurately describe the units after the received power adjusts instruction power dynamic response process. For different types of units, the dynamic response model not only regulates the upper and lower limits of capacity, climbing rate, frequency regulation delay, and other common parts [52] but also includes transmission links with their own energy conversion characteristics. Besides, AGC often uses the frequency domain model to describe the dynamic response process of the units, as shown in Figure 2. Td represents the secondary frequency regulation delay of the units; Gs represents the power response transfer function of the units, as shown in Table 1. T1∼T9 are the known parameters of the transfer function, respectively [53]. Therefore, according to the inverse Laplace transform of the frequency domain transfer function, the real-time output of the regulated power can be calculated through the input power, as follows:(1)ΔPioutt=L−1Giss1+Tdis⋅∑k=1Ne−ΔT⋅k−1s⋅Diink,(2)Diink=ΔPiink−ΔPiink−1,(3)ΔPioutk=ΔPioutt=k⋅ΔT,where i represents the ith AGC unit; K represents the kth discrete control period; ΔPiin and ΔPiout represent the input regulation power command and real-time output of regulated power of the ith AGC unit, respectively; and ΔT represents the control period of AGC, with a value of 1 to 16 seconds.Figure 2
Dynamic response models of different types of AGC units. (a) Traditional hydropower plant and thermal power station. (b) Wind farm and PV station.
(a)(b)Table 1
Dynamic response transfer functions of different AGC units.
TypeTransfer functionNonreheat turbine1/1+T1sReheat turbine1+T2s/1+T3s1+T4s1+T5sHydropower1−T6s1+T7s/1+0.5T6s1+T8sWind and solar renewable energy1/1+T9s


### 2.3. Optimization Mathematical Model

A performance index is assigned to quantify the total power response deviation, which is defined as the sum of the absolute deviation value of the regulating power command value and the power output value of all units, as follows:(4)minf=∑j=kN∑i=1nΔPiinj−ΔPioutj,where N represents the number of control periods and n represents the number of AGC units.In addition to considering the dynamic response transfer process of the units, some constraint conditions, such as power balance constraint, units capacity constraint, and climbing constraint, should also be considered in the power distribution process, as follows:(5)ΔPink=∑i=1nΔPiink,(6)ΔPink⋅ΔPiink≥0,i=1,2,…,n,(7)ΔPimin≤ΔPiink≤ΔPimax,i=1,2,…,n,(8)ΔPiink−ΔPiink−1≤ΔPirate,i=1,2,…,n,where ΔPin represents total power regulation command and ΔPirate represents the maximum climbing speed of the ith AGC units.


## 2.1. Control Framework

Generally, AGC has the following three control modes: (a) flat frequency control, (b) flat tie-line control, and (c) tie-line bias control. Also, tie-line bias control was used in this paper. The two-region interconnection power network is depicted in Figure1, in which the AGC control process of each region includes two parts: controller control and optimal power distribution. A PI controller is adopted to convert the power deviation and frequency deviation of the real-time receiving tie line into regional control deviation as input and output to the real-time adjusting power of the regional power grid ΔP. Then, ΔP is assigned to each AGC unit through the power allocation algorithm. Different from the traditional power system with hydropower plant and thermal power station as the main reserve capacity, wind farm and PV station no longer need to be disconnected from the grid but participate in the power AGC regulation of the system.Figure 1
Illustration of coordinated control of multisource for AGC under two-area framework.


## 2.2. Dynamic Response Model of Units

The establishment of dynamic response model of the units is mainly in order to more accurately describe the units after the received power adjusts instruction power dynamic response process. For different types of units, the dynamic response model not only regulates the upper and lower limits of capacity, climbing rate, frequency regulation delay, and other common parts [52] but also includes transmission links with their own energy conversion characteristics. Besides, AGC often uses the frequency domain model to describe the dynamic response process of the units, as shown in Figure 2. Td represents the secondary frequency regulation delay of the units; Gs represents the power response transfer function of the units, as shown in Table 1. T1∼T9 are the known parameters of the transfer function, respectively [53]. Therefore, according to the inverse Laplace transform of the frequency domain transfer function, the real-time output of the regulated power can be calculated through the input power, as follows:(1)ΔPioutt=L−1Giss1+Tdis⋅∑k=1Ne−ΔT⋅k−1s⋅Diink,(2)Diink=ΔPiink−ΔPiink−1,(3)ΔPioutk=ΔPioutt=k⋅ΔT,where i represents the ith AGC unit; K represents the kth discrete control period; ΔPiin and ΔPiout represent the input regulation power command and real-time output of regulated power of the ith AGC unit, respectively; and ΔT represents the control period of AGC, with a value of 1 to 16 seconds.Figure 2
Dynamic response models of different types of AGC units. (a) Traditional hydropower plant and thermal power station. (b) Wind farm and PV station.
(a)(b)Table 1
Dynamic response transfer functions of different AGC units.
TypeTransfer functionNonreheat turbine1/1+T1sReheat turbine1+T2s/1+T3s1+T4s1+T5sHydropower1−T6s1+T7s/1+0.5T6s1+T8sWind and solar renewable energy1/1+T9s


## 2.3. Optimization Mathematical Model

A performance index is assigned to quantify the total power response deviation, which is defined as the sum of the absolute deviation value of the regulating power command value and the power output value of all units, as follows:(4)minf=∑j=kN∑i=1nΔPiinj−ΔPioutj,where N represents the number of control periods and n represents the number of AGC units.In addition to considering the dynamic response transfer process of the units, some constraint conditions, such as power balance constraint, units capacity constraint, and climbing constraint, should also be considered in the power distribution process, as follows:(5)ΔPink=∑i=1nΔPiink,(6)ΔPink⋅ΔPiink≥0,i=1,2,…,n,(7)ΔPimin≤ΔPiink≤ΔPimax,i=1,2,…,n,(8)ΔPiink−ΔPiink−1≤ΔPirate,i=1,2,…,n,where ΔPin represents total power regulation command and ΔPirate represents the maximum climbing speed of the ith AGC units.


## 3. Memetic Salp Swarm Algorithm

### 3.1. Inspiration

Salps are marine creatures that resemble jellyfish in body structure and movement. They are usually joined end to end to form a chain, also known as a salp chain. The leader is located at the first end of the chain and has the best judgment of environment and food source. The remaining salps are followers, who follow the previous one in turn. This movement mode is conducive to the rapid coordinated movement and feeding of the salps group. Literature [54] established a mathematical model of salps chain in 2017 and proposed SSA to solve a series of optimization problems.


### 3.2. Optimization Framework

This paper is combined with cultural genetic algorithm, aiming at the shortcomings of SSA algorithm to improve, and defined as MSSA. The optimization process is as follows: the culture of each salp is defined as a solution to the optimization problem. All salps in the community are divided into different populations in the unit of salp chain, and each salp chain has its own culture and independently searches for food sources. At the same time, the culture of each salp affects and is influenced by other individuals and evolves with the evolution of the population. When the population evolution reaches a certain stage, the whole community will exchange information to realize the mixed evolution of the population until the convergence condition of the optimization problem is satisfied.The optimization framework of MSSA is shown in Figure3, which mainly includes the following two operations, as follows [54]:(a)
Local search in each chain: each salp chain will implement a local search to improve the ability of global exploration and local exploitation;(b)
Global coordination in virtual population: each salp is taken as a virtual population, and the population will be regrouped into multiple new salp chains. Hence, a global coordination can be achieved.Figure 3
Optimization framework of MSSA.


### 3.3. Local Search in Each Chain

The salp chain can be divided into two roles, i.e., the leader and the follower. It is worth noting that the leader is responsible for directing the entire salp chain to the food source, following each other. For themth salp chain, the position of the leader can be updated, as follows [55]:(9)xm1j=Fmj+c1c2ubj−lbj+lbj,c3≥0,Fmj−c1c2ubj−lbj+lbj,c3<0,where the superscript j represents the jth dimension of the searching space; xmij is the position of the leader in the mth salp chain; Fmj denotes the position of the food source, i.e., the current best solution obtained by the mth salp chain; and ubj and lbj are the upper and lower bounds of the jth dimension, respectively; c1=2e−4k/kmax2, where k is the current iteration number and kmax is the maximum iteration number, respectively. Besides, c2 and c3 are the random numbers, respectively, and c2, c3∈0,1 [56].In addition, the position of the followers can be updated, as follows [55]:(10)xmij=12xmij+xm,i−1j,i=2,3,…,n,m=1,2,…,M,where xmij is the position of the ith salp in the mth salp chain; n is the population size of each salp chain; and M is the number of salp chains, respectively.


### 3.4. Global Coordination in Virtual Population

To achieve a global coordination, the virtual population will be regrouped into different salp chains according to the salps’ fitness values, as shown in Figure4. Specifically, all salps are sorted according to the order of fitness from large to small. Finally, the salps will be divided into Mth salp chains, and the distribution rules are as follows: the first salp regroups into the first chain, the Mth salps into the Mth chain, the M+1th into the first chain, and so on. The update rules of the Mth salp chain are described as follows [57]:(11)Ym=xmi,fmi|xmi=Xm+Mi−1,fmi=Fm+Mi−1,i=1,2,…,n,m=1,2,…,M,where xmij is the position vector of the ith salp in the mth chain; fmi is the fitness value of the ith salp in the mth chain; F is the fitness value set of all the salps following the descending order; and X is the corresponding position vector set of all the salps, i.e., a position matrix, as follows:(12)X=x11x12x21x22⋯x1d⋯x2d⋮⋮xn×M1xn×M2⋱⋮⋯xn×Md,where d is the number of dimensions, and each row vector represents the position vector of each salp.Figure 4
Regroup operation of virtual population.In addition, the overall flowchart of MSSA is given in Figure5.Figure 5
The general procedure of MSSA.


### 3.5. MSSA Design for AGC System

In order to ensure that the initial solution is a feasible solution, this paper forms the initial solution according to an engineering method called PROP method. In other words, the initial solution is obtained by the same adjustable capacities’ ratio.On the other hand, this paper applies the penalty function method to the fitness functionFitx to satisfy the constraint conditions (5)–(8), as follows:(13)Fitx=fx,if constraints are satisfied,fx+MZu−Zulim2,otherwise,where M is the penalty function factor, and its value is usually a relatively large positive number; Zu is the uth constraint; and Zulim is the limit of the uth constraint.


## 3.1. Inspiration

Salps are marine creatures that resemble jellyfish in body structure and movement. They are usually joined end to end to form a chain, also known as a salp chain. The leader is located at the first end of the chain and has the best judgment of environment and food source. The remaining salps are followers, who follow the previous one in turn. This movement mode is conducive to the rapid coordinated movement and feeding of the salps group. Literature [54] established a mathematical model of salps chain in 2017 and proposed SSA to solve a series of optimization problems.


## 3.2. Optimization Framework

This paper is combined with cultural genetic algorithm, aiming at the shortcomings of SSA algorithm to improve, and defined as MSSA. The optimization process is as follows: the culture of each salp is defined as a solution to the optimization problem. All salps in the community are divided into different populations in the unit of salp chain, and each salp chain has its own culture and independently searches for food sources. At the same time, the culture of each salp affects and is influenced by other individuals and evolves with the evolution of the population. When the population evolution reaches a certain stage, the whole community will exchange information to realize the mixed evolution of the population until the convergence condition of the optimization problem is satisfied.The optimization framework of MSSA is shown in Figure3, which mainly includes the following two operations, as follows [54]:(a)
Local search in each chain: each salp chain will implement a local search to improve the ability of global exploration and local exploitation;(b)
Global coordination in virtual population: each salp is taken as a virtual population, and the population will be regrouped into multiple new salp chains. Hence, a global coordination can be achieved.Figure 3
Optimization framework of MSSA.


## 3.3. Local Search in Each Chain

The salp chain can be divided into two roles, i.e., the leader and the follower. It is worth noting that the leader is responsible for directing the entire salp chain to the food source, following each other. For themth salp chain, the position of the leader can be updated, as follows [55]:(9)xm1j=Fmj+c1c2ubj−lbj+lbj,c3≥0,Fmj−c1c2ubj−lbj+lbj,c3<0,where the superscript j represents the jth dimension of the searching space; xmij is the position of the leader in the mth salp chain; Fmj denotes the position of the food source, i.e., the current best solution obtained by the mth salp chain; and ubj and lbj are the upper and lower bounds of the jth dimension, respectively; c1=2e−4k/kmax2, where k is the current iteration number and kmax is the maximum iteration number, respectively. Besides, c2 and c3 are the random numbers, respectively, and c2, c3∈0,1 [56].In addition, the position of the followers can be updated, as follows [55]:(10)xmij=12xmij+xm,i−1j,i=2,3,…,n,m=1,2,…,M,where xmij is the position of the ith salp in the mth salp chain; n is the population size of each salp chain; and M is the number of salp chains, respectively.


## 3.4. Global Coordination in Virtual Population

To achieve a global coordination, the virtual population will be regrouped into different salp chains according to the salps’ fitness values, as shown in Figure4. Specifically, all salps are sorted according to the order of fitness from large to small. Finally, the salps will be divided into Mth salp chains, and the distribution rules are as follows: the first salp regroups into the first chain, the Mth salps into the Mth chain, the M+1th into the first chain, and so on. The update rules of the Mth salp chain are described as follows [57]:(11)Ym=xmi,fmi|xmi=Xm+Mi−1,fmi=Fm+Mi−1,i=1,2,…,n,m=1,2,…,M,where xmij is the position vector of the ith salp in the mth chain; fmi is the fitness value of the ith salp in the mth chain; F is the fitness value set of all the salps following the descending order; and X is the corresponding position vector set of all the salps, i.e., a position matrix, as follows:(12)X=x11x12x21x22⋯x1d⋯x2d⋮⋮xn×M1xn×M2⋱⋮⋯xn×Md,where d is the number of dimensions, and each row vector represents the position vector of each salp.Figure 4
Regroup operation of virtual population.In addition, the overall flowchart of MSSA is given in Figure5.Figure 5
The general procedure of MSSA.


## 3.5. MSSA Design for AGC System

In order to ensure that the initial solution is a feasible solution, this paper forms the initial solution according to an engineering method called PROP method. In other words, the initial solution is obtained by the same adjustable capacities’ ratio.On the other hand, this paper applies the penalty function method to the fitness functionFitx to satisfy the constraint conditions (5)–(8), as follows:(13)Fitx=fx,if constraints are satisfied,fx+MZu−Zulim2,otherwise,where M is the penalty function factor, and its value is usually a relatively large positive number; Zu is the uth constraint; and Zulim is the limit of the uth constraint.


## 4. Case Studies

The proposed methodology is tested on the IEEE load frequency control model. It is worth noting that 1 AGC unit in region A is increased to 5 units, as shown in Figure1. Besides, Tables 2 and 3 show the main parameters of response transfer and power regulation of the units, respectively, and the control period of AGC is set to 4 seconds. Also, the response performance is compared to that of PROP, GA, and SSA. To achieve a fair comparison, the population size is set to 10 and the maximum number of iterations is set to 100. It is worth noting that if the parameters are not chosen properly, the convergence time will be too long or the local optimum will be trapped. Besides, ode23 was selected as the solver, and the sampling rate was set to 0.001 s.Table 2
Parameters of dynamic response transfer functions of AGC units.
UnitsTypeParameters of transfer functionG1Thermal power stationT2 = 5, T3 = 0.08, T4 = 10, T5 = 0.3G2LNG unitsT2 = 2, T3 = 0.05, T4 = 5, T5 = 0.2G3Hydropower plantT6 = 1, T7 = 5, T8 = 0.513G4Wind farmT1 = 0.01G5PV stationT1 = 0.01Table 3
Main parameters of power regulation of AGC units.
UnitsTd (s)ΔPrate (MW/min)ΔPmax (MW)ΔPmin (MW)G1603050−50G2201830−30G3515020−10G41—15−5G51—10−10


### 4.1. Convergence Research

Figure6 shows the convergence curves of MSSA under different power distribution instructions. It can be found that MSSA can solve the optimal solution with high quality after 30 steps of iteration, and the subsequent iteration only makes slight reduction adjustment, which also indicated the fast convergence of the algorithm. In order to better evaluate the quality of the optimal solution of different methods, Table 4 shows the comparison of the convergence results of different methods, in which each indicator unit is MW. It shows that PROP performs power allocation according to the ratio of the adjustable capacity of the unit, so the power output of thermal power units with larger regulating capacity is relatively high, and it will lead to a large power deviation. On the other hand, GA, SSA, and MSSA can significantly reduce power deviation after their respective optimization operation, and MSSA has the best performance.Figure 6
Convergence curves of MSSA. (a)ΔPin = 80 MW. (b) ΔPin = −50 MW.
(a)(b)Table 4
Comparison on convergence results obtained by different methods.
ΔPinMethodΔP1inΔP2inΔP3inΔP4inΔP5inPower deviation80PROP32.0019.2012.809.606.40676.95GA15.5223.1817.9212.710.68529.86SSA10.6125.7417.2614.6711.72471.48MSSA14.2620.8419.5215.1210.26402.451−50PROP−23.81−14.29−4.76−2.38−4.76479.83GA−4.26−23.24−10.02−3.86−8.62274.26SSA−8.65−20.62−9.85−4.29−6.59307.48MSSA−4.48−20.69−9.14−5.58−10.11234.254


### 4.2. Online Optimization

In order to evaluate the online optimization performance of MSSA, in the area of A, a step power perturbationΔPL = 80 MW has occurred. The online optimization results of MSSA and PROP are compared as shown in Figure 7. Compared with previous optimization, the power deviation obtained by MSSA is smaller, and the overshoot of the total power instruction can be avoided. In addition, wind farm and PV station can recover the power system disturbed by high power quickly in the initial stage of power disturbance because of their fast response speed.Figure 7
Result comparison of online optimization by MSSA and without optimization whenΔPL = 80 MW. (a) Total power regulation curve. (b) Curve of ACE. (c) Curve of CPS1. (d) Frequency deviation Δf. (e) Active power of tie line (MW). (f) Power output curve of the units.
(a)(b)(c)(d)(e)(f)To further test the optimization performance of different algorithms, Table5 presents the online optimization results of different algorithms (the optimal value is marked in bold), in which area control error (ACE), Δf1, and CPS1 are the average values in the simulation. Besides, the power deviation is the total deviation in the simulation time. Accuracy is used to measure the fitting degree of the actual adjustment output and the adjustment instruction curve in the simulation time. It can be found that, compared with PROP, the other three methods can significantly reduce power deviation, thus significantly improving the system’s dynamic response performance. Compared with GA and SSA, the online optimization results of MSSA are better, which is because memetic computing framework can observably improve the ability of exploration and exploitation.Table 5
Result comparison of online optimization obtained by different methods.
ΔPLAlgorithms|ACE| (MW)Δf1 (Hz)CPS1 (%)Accuracy (%)Power deviation (MW)80PROP7.76780.0356194.6995.84630.05GA7.42090.0323195.8698.25261.32SSA7.50950.0321195.8498.17267.52MSSA7.39950.0316195.5798.84254.21−50PROP5.05410.0235197.7894.99477.46GA4.60460.0203198.2598.21164.26SSA4.60440.0202198.2998.23161.48MSSA4.60400.0196198.3698.41149.96


## 4.1. Convergence Research

Figure6 shows the convergence curves of MSSA under different power distribution instructions. It can be found that MSSA can solve the optimal solution with high quality after 30 steps of iteration, and the subsequent iteration only makes slight reduction adjustment, which also indicated the fast convergence of the algorithm. In order to better evaluate the quality of the optimal solution of different methods, Table 4 shows the comparison of the convergence results of different methods, in which each indicator unit is MW. It shows that PROP performs power allocation according to the ratio of the adjustable capacity of the unit, so the power output of thermal power units with larger regulating capacity is relatively high, and it will lead to a large power deviation. On the other hand, GA, SSA, and MSSA can significantly reduce power deviation after their respective optimization operation, and MSSA has the best performance.Figure 6
Convergence curves of MSSA. (a)ΔPin = 80 MW. (b) ΔPin = −50 MW.
(a)(b)Table 4
Comparison on convergence results obtained by different methods.
ΔPinMethodΔP1inΔP2inΔP3inΔP4inΔP5inPower deviation80PROP32.0019.2012.809.606.40676.95GA15.5223.1817.9212.710.68529.86SSA10.6125.7417.2614.6711.72471.48MSSA14.2620.8419.5215.1210.26402.451−50PROP−23.81−14.29−4.76−2.38−4.76479.83GA−4.26−23.24−10.02−3.86−8.62274.26SSA−8.65−20.62−9.85−4.29−6.59307.48MSSA−4.48−20.69−9.14−5.58−10.11234.254


## 4.2. Online Optimization

In order to evaluate the online optimization performance of MSSA, in the area of A, a step power perturbationΔPL = 80 MW has occurred. The online optimization results of MSSA and PROP are compared as shown in Figure 7. Compared with previous optimization, the power deviation obtained by MSSA is smaller, and the overshoot of the total power instruction can be avoided. In addition, wind farm and PV station can recover the power system disturbed by high power quickly in the initial stage of power disturbance because of their fast response speed.Figure 7
Result comparison of online optimization by MSSA and without optimization whenΔPL = 80 MW. (a) Total power regulation curve. (b) Curve of ACE. (c) Curve of CPS1. (d) Frequency deviation Δf. (e) Active power of tie line (MW). (f) Power output curve of the units.
(a)(b)(c)(d)(e)(f)To further test the optimization performance of different algorithms, Table5 presents the online optimization results of different algorithms (the optimal value is marked in bold), in which area control error (ACE), Δf1, and CPS1 are the average values in the simulation. Besides, the power deviation is the total deviation in the simulation time. Accuracy is used to measure the fitting degree of the actual adjustment output and the adjustment instruction curve in the simulation time. It can be found that, compared with PROP, the other three methods can significantly reduce power deviation, thus significantly improving the system’s dynamic response performance. Compared with GA and SSA, the online optimization results of MSSA are better, which is because memetic computing framework can observably improve the ability of exploration and exploitation.Table 5
Result comparison of online optimization obtained by different methods.
ΔPLAlgorithms|ACE| (MW)Δf1 (Hz)CPS1 (%)Accuracy (%)Power deviation (MW)80PROP7.76780.0356194.6995.84630.05GA7.42090.0323195.8698.25261.32SSA7.50950.0321195.8498.17267.52MSSA7.39950.0316195.5798.84254.21−50PROP5.05410.0235197.7894.99477.46GA4.60460.0203198.2598.21164.26SSA4.60440.0202198.2998.23161.48MSSA4.60400.0196198.3698.41149.96


## 5. Conclusions

In this paper, a multisource optimal collaborative control method for power system with renewable energy participation in secondary frequency regulation is proposed. The main contributions can be summarized as follows:(i)
A novel AGC control model is established for the power system with high renewable energy penetration to improve the dynamic response performance of the system;(ii)
The study verified the effectiveness of MSSA for AGC. It can not only meet the online regulation requirements of AGC but also quickly obtain high-quality regulation schemes with high convergence stability, and the dynamic response performance of the entire regional power grid is significantly improved.Besides, electric vehicles will be considered to participate in AGC in future studies.



---

*Source: 6661793-2020-12-14.xml*