Load Frequency Control of Pumped Storage Power Station Based on LADRC

: The pumped storage power station has the characteristics of frequency-phase modulation, energy saving, and economy, and has great development prospects and application value. In order to cope with the large-scale integration and intermittency of renewable energy and improve the ability of pumped storage units to participate in power grid frequency modulation, this paper proposed a load frequency control (LFC) strategy for pumped storage units based on linear active disturbance rejection technology. Firstly, based on the operating characteristics of the pumped storage power station, the LFC model of the two-area reheat steam turbine under nonlinear conditions such as governor dead zone and generation rate constrains was established. Secondly, a second-order linear active disturbance rejection control (LADRC) was designed. The feasibility and control performance of the proposed LFC system were quantitatively analyzed through simulation. The results show that the LADRC has better control e ﬀ ect and stronger robustness than fractional-order proportion integration di ﬀ erentiation (FOPID) and traditional proportion integration di ﬀ erentiation (PID) controller. Finally, the pumped storage power station was added, and it was found that it has better correction performance under both generating and pumping operations, which greatly improved the dynamic response of secondary frequency modulation.


Introduction
With the gradual depletion of traditional primary energy sources such as coal and oil, energy issues have become increasingly prominent and have become the focus of current power systems. Green renewable energy has obvious volatility and intermittent. Due to these characteristics, when large-scale renewable energy is integrated into the power grid, it may adversely affect power quality, safety, and stable operation of the power grid [1,2]. Therefore, effective load frequency control (LFC) [3] plays an indispensable role in the safety of power grid and users, power quality assessment and power monitoring.
The majority of scholars at home and abroad have proposed different advanced control methods for LFC, such as load frequency predictive control with wind power connection, sliding mode control, robust control based on linear matrix inequality (LMI) and delay marginal estimation, fuzzy control, LFC based on differential game theory, etc. But at present, there are few researches on LFC of the unit models are studied, and a two-area reheat steam turbine LFC model considering nonlinear factors is established. A second-order LADRC is designed and compared with traditional PID and fractional-order PID (FOPID) [19] controllers. According to the characteristics of the two-stage work of pumped storage power station, the effects of generating and pumping operations on system frequency regulation are studied. Simulation results show: (1) The LADRC shows stronger control effect and better dynamic performance. (2) After joining the pumped storage power station, the frequency stability of the area has been significantly improved.
The rest of the paper is organized as follows. In Section 2, the system model is described in detail, including the unit models, the area control error models, and the nonlinear links in the unit. In Section 3, a model of pumped storage power station is built. In Section 4, a LADRC model is designed. In Section 5, the simulations are performed and the results are analyzed. Conclusions are presented in Section 6.

The Unit Model
The unit models include prime movers and governors. Prime movers are steam turbine and hydraulic turbine [11,20].
The reheat steam turbine is selected for steam turbine, and its model is given in Equation (1): The hydraulic turbine model is given in Equation (2): where K r is the ratio of the steam generated in the high-pressure cylinder segment to the total power of the turbine; T r is the time constant of reheater; T t is the time constant of steam chamber and the main intake steam volume; T w is the water starting time; s is the Laplace transform operator. The model of steam turbine governor is given in Equation (3): The governor with transient slope compensation given in Equation (4) is often used in traditional hydraulic turbines, which can effectively suppress the strong hysteresis of hydraulic turbines. The new hydraulic turbine governors in this paper adopt the Digital Electric Hydraulic Control System (DEH), which is given in Equation (5) [12]: where R is the difference adjustment coefficient of the hydraulic turbine; T R is the reset time constant; T g is the governor time constant; R T is the transient rate of decline; K p , K i , K d are the proportional, integral and differential gains of the DEH; f is the system frequency.

Area Control Error
The balance between the active power output of the generator set and the load of the grid is related to the stability of frequency. At present, AGC is widely used in secondary frequency modulation of power grid to adjust the output of generator units in the region, so as to ensure the stability of system frequency and content the demand of load change. In this paper, the strategy of calculating area control error mainly adopts the following two methods [21].

1.
Flat Frequency Control (FFC): This mode is mainly applied to the main system or independent system of the integrated system. The system frequency deviation coefficient β is defined in the Equation (7).
In order to make the system frequency error and Area Control Error (ACE) zero, system control signal takes ACE as input; R i as adjustment coefficient; and load damping coefficient D i is 8.33 × 10 −3 p.u.MW/Hz; β can be obtained from Equation (7).

2.
Tie-Line and Frequency Bias Control (TBC): TBC is the most commonly used grid control mode. In Equations (6)- (8), ACE is the area control error; f is the system frequency deviation when a disturbance occurs; and P tie is the exchange power deviation of tie-lines.

Nonlinear Links
The nonlinear links in the LFC system mainly including dead zone of governor, GRC and transmission delay. In this paper, the original governor link is replaced by the governor link with a dead zone after linearization, and a limit device is added to the output variation of reheat steam turbine to realize the speed constraint of the generator, which is more in line with the operating condition of the turbine unit in the actual power system. At the same time, the ACE is taken as the control target, and it is taken as the input of the system error controller to ensure that ACE is zero when the system is in a stable state. Thus, the delay caused by the control signal in the transmission process is considered.

Dead Zone of Governor
The governors in the power system are all equipped with a certain dead zone, whose purpose is to reduce the frequent action of the governors caused by the frequency deviation in the power grid, so as to protect the governors and extend their service life. Due to the hysteresis and nonlinearity of the governor, the dead zone is generally described by the gap characteristic link. Figure 1 shows the characteristics of the governor's dead zone, k is the slope. Assuming that the input of this nonlinear element is X and the output is Y, the relationship between the input and output is given in Equation (9) [22]. infinity, N1 and N2 tend to be constant. Substituting Equation (13) into Equation (11) and performing a Laplace change, combining with Equation (3), the transfer function of the governor with a dead zone after linearization can be obtained:

Generation Rate Constrains
GRC, also known as Unit Ramp Control (URC), refers to the unit output is limited by a certain maximum rate of change [23]. In the power grid, the given maximum change value is generally set as 0.0017 p.u.MW/s. In this paper, a limit device is added to the output variation of reheat steam turbine, and the limit device value is also 0.0017 p.u.MW/s. The reheat steam turbine model considering GRC is shown in Figure 2.  When the function Y has a sinusoidal relationship, there is:

Tie-Line Model
where A is the amplitude of the sinusoidal input signal and ω 0 is the frequency of the sinusoidal input signal. Linearize Y with the descriptive function method, then expand F X, .
X with the Fourier series, and take the first three terms: Among them, the various coefficients are: By integrating the Equation (12), it can be obtained that each coefficient satisfies the following relationship: where h is the dead zone setting value of the governor. When A = h, N 1 = N 2 = 0. When A/h increases, N 1 increases monotonically, and N 2 also increases after falling to A/h = 2. When A/h approaches infinity, N 1 and N 2 tend to be constant. Substituting Equation (13) into Equation (11) and performing a Laplace change, combining with Equation (3), the transfer function of the governor with a dead zone after linearization can be obtained:

Generation Rate Constrains
GRC, also known as Unit Ramp Control (URC), refers to the unit output is limited by a certain maximum rate of change [23]. In the power grid, the given maximum change value is generally set as 0.0017 p.u.MW/s. In this paper, a limit device is added to the output variation of reheat steam turbine, and the limit device value is also 0.0017 p.u.MW/s. The reheat steam turbine model considering GRC is shown in Figure 2. GRC, also known as Unit Ramp Control (URC), refers to the unit output is limited by a certain maximum rate of change [23]. In the power grid, the given maximum change value is generally set as 0.0017 p.u.MW/s. In this paper, a limit device is added to the output variation of reheat steam turbine, and the limit device value is also 0.0017 p.u.MW/s. The reheat steam turbine model considering GRC is shown in Figure 2.

Tie-Line Model
The tie-line is the transmission line that connects different grid areas. Ignoring the line loss, the flow power of the tie-line from area i to j can be expressed as [11]: where P tieij is the flow power of the tie-line; V i , V j are the voltage at both ends of the tie-line; θ i , θ j are the voltage angle at both ends of the tie-line; X ij is the contact line circuit reactance; P * i is the power of area i. The power micro-increment model of tie-line can be obtained from Equation (16), as given in Equation (17): where θ is the voltage angle micro-increment of the tie-line; P tieij is the power micro-increment of the tie-line; f i , f j are the frequency variation values of area i and j. The transfer function model of tie-line power deviation can be obtained by Laplacian variation of Equation (17), as given in Equation (18): T ij is defined as the contact line synchronization coefficient, and its calculation formula is given in Equation (19):

Two-Area LFC Model
The model of regional LFC system is obtained by combining the common component transfer functions described in detail in the previous section with the actual operation mode of the power grid. The LFC model with two single areas is connected with tie-lines to obtain the LFC model of the two-area reheat steam turbine with nonlinear links, as shown in Figure 3.

Two-Area LFC Model
The model of regional LFC system is obtained by combining the common component transfer functions described in detail in the previous section with the actual operation mode of the power grid. The LFC model with two single areas is connected with tie-lines to obtain the LFC model of the two-area reheat steam turbine with nonlinear links, as shown in Figure 3.   In the figure, ACE is used as the input, and the system frequency deviation is used as the output; P Di and P Dj are the load disturbance values exerted by area i and j; B is the Frequency deviation coefficient; a ij is referred to as the synchronous power coefficient, which can be obtained from Equation (20).
where P * i and P * j are the rated power of the generator set in control area i and j respectively, and the parameters of tie-line are set is given in Equation (21).
The transfer function relationship between generator and system connection in power grid is: f D , P g is the change of generator output power; P d is the change of load power; D is the load damping constant, which represents the load change caused by the change of unit frequency; H is the unit inertia constant defined by the standard value.

Model of Pumped Storage Power Station
Based on the fact that the main power plants in the power grid are all thermal power generation, the LFC model of pumped storage power station in this paper is based on the LFC model of reheat steam turbine in two-area, supplemented by the links of pumped storage power station. A two-stage pumped storage power model is established which can be divided into two basic operating conditions: generating operation and pumping operation. The generating operation is suitable for the peak load period. The hydraulic turbine operation causes the water flow in the upper reservoir to generate electricity in the lower reservoir, thus providing output for the power system. The pumping operation is suitable for the load trough period. At this time, the system relies on the surplus electricity to pump water from the lower reservoir, so as to reduce the output peak-valley difference of the power system, thus reducing the frequency of start and stop of the pumped storage units and the output variation, and maintaining the normal and stable operation of the units.
In this paper, the generating operation model of pumped storage power station is simplified as a hydraulic turbine unit, and the hydraulic turbine model simulation is used as an external power supply of the hydropower station. In pumping operation, it is equivalent to a load, considering the amount of power to support the grid when it stops pumping. Among them, the unit primary frequency regulation dead zone, ACE dead zone limit, maximum load adjustment limit and unit capacity limit are not considered. However, we consider some of the balance and constraints of the pumped storage power station: Water balance of upper reservoir: Balance of total generating water and total pumping water during the planning period: Water balance of lower reservoir: Flow limit under generating operation: Flow limit under pumping operation: Capacity limit of upper reservoir: Capacity limit of lower reservoir: where Q is the water flow and V is the storage water volume of the reservoir. When the operation conditions of the pumped storage power station meet the above constraints, the two-stage pumped storage power model is established considering the two working states of generating and pumping operation of the pumped storage power station, is given in Equation (30). Structure diagram as shown in the dashed part of Figure 3.
where P pd is the power transmitted to the power grid when the pumping operation stops; G d (s) is the transfer function model of hydraulic turbine governor under the generating operation, as given in Equation (4) or Equation (5); G t (s) is the transfer function model of the hydraulic turbine, as given in Equation (2). In the Figure 3, a = {1,2 . . . n}, n is the number of pumped storage units working under pumping operation.

Controller Design
The structure of the ADRC main components include: tracking differentiator (TD), extended state observer (ESO), and nonlinear state error feedback (NLSEF) [24]. However, the nonlinear form of ADRC needs to adjust multiple parameters, and these parameters affect each other. Based on the idea of ADRC, a linear auto disturbance rejection control (LADRC) method is proposed in [25], taking the controller bandwidth as the debugging parameter of the control performance. This method simplified the calculation, made it easy to set and maintain, and had the same control performance as traditional ADRC.
LADRC is composed of three important components: linear extended state observer (LESO), linear feedback control law, and disturbance compensation. Figure 4 shows the second-order LADRC structure adopted from [26], where r and y are the reference input and output of the system; G p (s) is the controlled object; k P and k d are the linear feedback control rate parameters; d is the external disturbance of the system; and b 0 is the estimation of the input gain b of the system in the controller. The principle and function of each structure are introduced below.
in Equation (2). In the Figure 3, a = {1,2…n}, n is the number of pumped storage units working under pumping operation.

Controller Design
The structure of the ADRC main components include: tracking differentiator (TD), extended state observer (ESO), and nonlinear state error feedback (NLSEF) [24]. However, the nonlinear form of ADRC needs to adjust multiple parameters, and these parameters affect each other. Based on the idea of ADRC, a linear auto disturbance rejection control (LADRC) method is proposed in [25], taking the controller bandwidth as the debugging parameter of the control performance. This method simplified the calculation, made it easy to set and maintain, and had the same control performance as traditional ADRC.
LADRC is composed of three important components: linear extended state observer (LESO), linear feedback control law, and disturbance compensation. Figure 4 shows the second-order LADRC structure adopted from [26], where r and y are the reference input and output of the system; Gp(s) is the controlled object; kP and kd are the linear feedback control rate parameters; d is the external disturbance of the system; and b0 is the estimation of the input gain b of the system in the controller. The principle and function of each structure are introduced below.

Linear Extended State Observer
LESO is the core part of the whole LADRC, which can observe the estimated values Z1 and Z2 of the state variable and the total disturbance Z3 of the system according to the input and output data without relying on the specific mathematical model of generating disturbance.

Linear Extended State Observer
LESO is the core part of the whole LADRC, which can observe the estimated values Z 1 and Z 2 of the state variable and the total disturbance Z 3 of the system according to the input and output data without relying on the specific mathematical model of generating disturbance.
The differential representation of the controlled second-order object G p (s) is given in Equation (31): ..
where u is the controlled object input; y is the output; . y, ..
y are the first and second order differentials of the output; w is the unmeasured disturbance of the system; b is the input gain; g is the total influence of the system's own perturbation and external disturbance, and Equation (31) can also be expressed as: ..
where l = g + (b − b 0 ) × u is the total disturbance of the system. Parameter b has a small range of change, and the approximate estimated gain value b 0 is used instead of b. At this time, the estimation error of b 0 on b is also regarded as part of the disturbance. When b 0 = b, it will not cause too much error to the observation state. Generally, take b 0 greater than b. Let y, x 3 = l, l be defined as an expansion state of the system, the total disturbance of the system is expanded to a state variable x 3 , so that the original second-order system is expanded into a new third-order linear system. Assuming that l is differentiable and .
x 3 = h is defined, Equation (32) can be expressed in the state space as follows: l. The estimated values of y, . y and l can be obtained by using the LESO: y and l, while β 1 , β 2 , β 3 are the gain parameters of the LESO to be determined to control the accuracy of the estimated values. Since the accuracy of the estimated value of the LESO is highly correlated with its bandwidth, let ω 0 be the bandwidth of the LESO, according to the work in [27], the observer gain parameters can be selected according to Equations (35) and (36): Then the characteristic equation of LESO can be configured into the form of Equation (37): The parameterized gain coefficient of LESO is expressed as: Therefore, ω 0 is the only parameter to be adjusted in the LESO. This configuration can not only ensure the stability of the system, but also give a better transition process, which makes the calculation easy. The range of values takes into account the controllability and noise resistance of the system. After ω 0 is determined, β 1 , β 2 , β 3 can be determined correspondingly according to Equation (38) [28].

Disturbance Compensation
As shown in Figure 4, the disturbance compensation link is set as follows: where u 0 is the virtual control amount. If the LESO can achieve z 3 = l, the Equation (39) is brought into the Equation (32), and .. y = u 0 is obtained, that is, the controlled object is converted into a double integrator series structure for easy control.

Linear Feedback Control
NLSEF calculates the control signal u 0 according to the system state error. The following proportional higher order differential control rate is used: Reference [29,30] analyzed the stability of LADRC in uncertain time-varying systems, and introduced that LADRC can achieve the following two goals: One is the convergence of the closed-loop system, and the other is that the dynamic performance of the closed-loop system can be optimized by adjusting the bandwidth of LESO. With closed-loop stability, the closed-loop equation of the control system can be obtained as given in Equation (41) by introducing Equation (40) into the double integrator series object Equation (39) after the disturbance compensation.
.. y + k d . y + k p y = k p r (41) By using the relationship between bandwidth and control performance, the characteristic polynomial of closed-loop Equation (41) is expressed in the form of ω c as given in Equation (42) if ω c is the controller bandwidth.
The controller parameters can be expressed as: ω c is the only parameter to be set in the controller, which simplifies the calculation of controller parameters. The value ranges between the response speed and stability of the system [27]. The values of k d and k P are determined according to the relationship between Equation (43) and ω c .

Simulation of Load Frequency Control in Single and Two-Area
In this paper, first of all, the LFC models of single-area and two-area reheat steam turbine units are established by MATLAB/Simulink (MATLAB R2018b, Math Works, Natick, MA, USA, 1984). FFC mode is adopted for single-area ACE, and TBC mode is adopted for two-area ACE. According to the model parameters of reheat steam turbine and generator power system in [15], the LFC parameters of reheat steam turbine unit in the two-area in this paper are obtained as shown in Table 1, and the control simulation module of regional reheat steam turbine unit is determined to reflect the corresponding frequency changes. The unit parameter values in the two-area are set consistently and the parameters of tie-line are given in Equation (21). A disturbance signal of 0.01 p.u.MW is added to area i at 1 s. The frequency variation chart of the single and two-area is shown in Figure 5. p.u.MW is added to area i at 1 s. The frequency variation chart of the single and two-area is shown in Figure 5.  The simulation results are shown in Figure 5. It can be seen from the simulation diagram that the system is disturbed at 1s, and the frequency decreases rapidly. The maximum frequency variation of the two-area is nearly 50% less than that of the single-area. The regulation time of the two-area is shortened by 5 to 6 s, which can restore to the stable state more quickly, but the frequency oscillation frequency increases. It can be seen that when a disturbance occurs in region i, the interconnected region j can provide power support to region i through tie-lines, helping the disturbed region to restore the frequency to a stable state more quickly. Compared with isolated network operation, the safety and stability of large power grid model are improved. The simulation results are shown in Figure 5. It can be seen from the simulation diagram that the system is disturbed at 1s, and the frequency decreases rapidly. The maximum frequency variation of the two-area is nearly 50% less than that of the single-area. The regulation time of the two-area is shortened by 5 to 6 s, which can restore to the stable state more quickly, but the frequency oscillation frequency increases. It can be seen that when a disturbance occurs in region i, the interconnected region j can provide power support to region i through tie-lines, helping the disturbed region to restore the frequency to a stable state more quickly. Compared with isolated network operation, the safety and stability of large power grid model are improved.

Simulation of Two-Area LFC Considering Nonlinearity
A two-area reheat steam turbine LFC model with nonlinear links is established as shown in Figure 3 (without dashed line part). The parameters of the unit model, tie-line and regional disturbance are set up in accordance with the previous section. The controller adopts traditional PID, FOPID, and LADRC. According to the parameter design of the controllers in [26], the cut-and-trial method is used to determine and select the LFC parameters suitable for the pumped storage power station in this paper, so as to effectively suppress the disturbance and restore the frequency stability. Its control parameters are given in Table 2, and the simulation time domain is 60 s. The frequency variation chart of the disturbance zone under different controllers is shown in Figure 6. Table 2. Basic parameters of controller.

Controller Parameters Value Parameters Value
Traditional proportion integration differentiation (PID) linear active disturbance rejection control (LADRC) The simulation results are shown in Figure 6. It can be seen from the figure that, after taking nonlinear factors into consideration, the amplitude of frequency variation of the system greatly increases to more than 3 times of the original value, and the frequency recovery speed also decreases. The reason is that the speed constraint of steam turbine in nonlinear factors restricts the rapid output of the units, resulting in the lack of power cannot be compensated in time, and the frequency recovery is slow. Compared with the three different controllers, LADRC has the best control effect, with the smallest frequency oscillation range and the shortest recovery time. The amplitude of oscillation decreases to approximately 0.016 Hz, and the recovery time is shortened nearly 20 s. The second is FOPID, and the worst is the traditional PID. The correctness and validity of the LADRC method proposed in this paper are verified.

Simulation of LFC under Random Disturbance
Considering the influence of random factors in the actual situation, disturbances are mostly unmeasured and uncertain clutter. Therefore, considering that the random disturbance time is between 10 s and 25 s, the disturbance amplitude is [-0.01 p.u.MW, 0.01 p.u.MW], and the simulation time domain is 50 s, the performance of LADRC, FOPID and traditional PID controller under uncertain disturbance is observed. The simulation results are shown in Figure 7. At this time, The simulation results are shown in Figure 6. It can be seen from the figure that, after taking nonlinear factors into consideration, the amplitude of frequency variation of the system greatly increases to more than 3 times of the original value, and the frequency recovery speed also decreases. The reason is that the speed constraint of steam turbine in nonlinear factors restricts the rapid output of the units, resulting in the lack of power cannot be compensated in time, and the frequency recovery is slow. Compared with the three different controllers, LADRC has the best control effect, with the smallest frequency oscillation range and the shortest recovery time. The amplitude of oscillation decreases to approximately 0.016 Hz, and the recovery time is shortened nearly 20 s. The second is FOPID, and the worst is the traditional PID. The correctness and validity of the LADRC method proposed in this paper are verified.

Simulation of LFC under Random Disturbance
Considering the influence of random factors in the actual situation, disturbances are mostly unmeasured and uncertain clutter. Therefore, considering that the random disturbance time is between 10 s and 25 s, the disturbance amplitude is [−0.01 p.u.MW, 0.01 p.u.MW], and the simulation time domain is 50 s, the performance of LADRC, FOPID and traditional PID controller under uncertain disturbance is observed. The simulation results are shown in Figure 7. At this time, the frequency is affected and oscillates constantly. However, LADRC has the shortest regulation time and the fastest recovery speed, and tends to be stable after 25 s. It can be seen from the oscillation amplitude that it is smaller than the other two controllers. It can be seen that LADRC proposed in this paper still has excellent control performance under unpredictable and irregular disturbances. rapid output of the units, resulting in the lack of power cannot be compensated in time, and the frequency recovery is slow. Compared with the three different controllers, LADRC has the best control effect, with the smallest frequency oscillation range and the shortest recovery time. The amplitude of oscillation decreases to approximately 0.016 Hz, and the recovery time is shortened nearly 20 s. The second is FOPID, and the worst is the traditional PID. The correctness and validity of the LADRC method proposed in this paper are verified.

Simulation of LFC under Random Disturbance
Considering the influence of random factors in the actual situation, disturbances are mostly unmeasured and uncertain clutter. Therefore, considering that the random disturbance time is between 10 s and 25 s, the disturbance amplitude is [-0.01 p.u.MW, 0.01 p.u.MW], and the simulation time domain is 50 s, the performance of LADRC, FOPID and traditional PID controller under uncertain disturbance is observed. The simulation results are shown in Figure 7. At this time, the frequency is affected and oscillates constantly. However, LADRC has the shortest regulation time and the fastest recovery speed, and tends to be stable after 25 s. It can be seen from the oscillation amplitude that it is smaller than the other two controllers. It can be seen that LADRC proposed in this paper still has excellent control performance under unpredictable and irregular disturbances.

Stability and Robustness Analysis of Different Control Methods
In this section, we evaluate various control methods in the form of tables under different error criteria. Table 3 comprehensively analyzes the simulation results of Figure 7, compares the influence of three control methods on frequency control results under load disturbance, and evaluates the stability and robustness of the three methods. As shown in Table 3, we select the four dimensions of disturbance duration, maximum amplitude, overshoot and mean square error (MSE) to evaluate the controller error. It is obvious that LADRC has the shortest disturbance time, the smallest maximum amplitude and MSE among the three. Compared with traditional PID, LADRC and FOPID have overshoot, but the overshoot is within the allowable range. Moreover, the MSE of LADRC is quite different from that of other controllers, compared with the traditional PID is 6.19, LADRC is only 3.85, almost double the relationship, which shows that LADRC has better control accuracy. Therefore, the anti-disturbance capability of LADRC is more remarkable. The disturbance time of FOPID is only 2 seconds slower than that of traditional PID, which is nearly the same, and its maximum amplitude is higher, but the MSE of observation is much smaller than the traditional PID. The comprehensive analysis shows that the FOPID data has a low degree of variation and a small error. Although the control effect is higher than the traditional PID, compared with LADRC, there is still a large disparity.

Generating Operation
Based on the LFC of reheat steam turbine, this paper builds the pumped storage energy LFC model under generating operation according to the dashed line in Figure 3. DEH is adopted as the governor of hydraulic turbine. The LADRC designed in this paper is adopted as the controller. According to the model parameters of hydraulic turbine in [15], the LFC parameters of the hydraulic turbine unit in this paper are shown in Table 4. It is determined that the frequency control simulation parameters of generating and pumping operation in this paper can reflect the corresponding frequency changes, so as to achieve the analogue simulation of frequency modulation function of pumped storage power station. The simulation time domain was set to 60 s, and a disturbance of 0.015 p.u.MW occurred in the system at 1 s. The frequency variation chart of the disturbance zone when the pumped storage power station works in the generating operation is shown in Figure 8, and the exchange power chart of the tie-line is shown in Figure 9.   Since the pumped storage power station working in the generating operation has been added to the hydraulic turbine, which is almost free from the restriction of GRC and can rapidly increase its output to make up for the power loss caused by the disturbance of the system, thus accelerating the frequency stability. Figure 8 intuitively shows that the maximum value of grid frequency decline is significantly reduced, and the recovery time is shortened by nearly 10 s. It can be seen that the frequency control of the disturbance zone after the pumped storage power station is more robust to interference, and the reliability of the LADRC method is further verified. From the change of the exchange power of the tie-line, after the pumped storage energy is added, the exchange  Since the pumped storage power station working in the generating operation has been added to the hydraulic turbine, which is almost free from the restriction of GRC and can rapidly increase its output to make up for the power loss caused by the disturbance of the system, thus accelerating the frequency stability. Figure 8 intuitively shows that the maximum value of grid frequency decline is significantly reduced, and the recovery time is shortened by nearly 10 s. It can be seen that the frequency control of the disturbance zone after the pumped storage power station is more robust to interference, and the reliability of the LADRC method is further verified. From the change of the exchange power of the tie-line, after the pumped storage energy is added, the exchange power of the tie-line increases, then decreases, and its recovery speed is obviously accelerated. Due Since the pumped storage power station working in the generating operation has been added to the hydraulic turbine, which is almost free from the restriction of GRC and can rapidly increase its output to make up for the power loss caused by the disturbance of the system, thus accelerating the frequency stability. Figure 8 intuitively shows that the maximum value of grid frequency decline is significantly reduced, and the recovery time is shortened by nearly 10 s. It can be seen that the frequency control of the disturbance zone after the pumped storage power station is more robust to interference, and the reliability of the LADRC method is further verified. From the change of the exchange power of the tie-line, after the pumped storage energy is added, the exchange power of the tie-line increases, then decreases, and its recovery speed is obviously accelerated. Due to the accelerated frequency recovery speed, the power of the tie-line inevitably produces certain overshoot, but the overshoot is within the allowable range.

Pumping Operation
The LFC model of the pumped storage power station under pumping operation is established according to the dashed line in Figure 3. The pumping power is set at 0.05, P pd = 0.05, a = 1. The parameter setting is consistent with the last section of the generating operation. The frequency variation chart of the disturbance zone is shown in Figure 10, and the exchange power chart of the tie-line is shown in Figure 11.  When load disturbance occurs in the pumping operation, the pumped storage power station can stop pumping after receiving the signal from the grid, cut off part of the load quickly, and support the power consumed by pumping to the grid to maintain the stability of the grid frequency. From the simulation results, it can be seen that the pumped storage power station working under the pumping operation obtains certain power support in a short time after the grid disturbance, and the frequency fluctuation range is shortened and restored to the stable state more quickly and smoothly. For the exchange power of the tie-line, the pumping operation reduces the value of the exchange power, significantly reduces the number of oscillations and amplitude. The recovery rate is more than twice that of without pumped storage energy, which is an effective supplement to the LFC of the power grid.

Conclusions
In this paper, an LFC method for a pumped storage unit based on LADRC is proposed. Firstly, the LFC of reheat steam turbine with nonlinear factors is modeled, and according to the working mode of the pumped storage power station, the power plant LFC model under generating and  When load disturbance occurs in the pumping operation, the pumped storage power station can stop pumping after receiving the signal from the grid, cut off part of the load quickly, and support the power consumed by pumping to the grid to maintain the stability of the grid frequency. From the simulation results, it can be seen that the pumped storage power station working under the pumping operation obtains certain power support in a short time after the grid disturbance, and the frequency fluctuation range is shortened and restored to the stable state more quickly and smoothly. For the exchange power of the tie-line, the pumping operation reduces the value of the exchange power, significantly reduces the number of oscillations and amplitude. The recovery rate is more than twice that of without pumped storage energy, which is an effective supplement to the LFC of the power grid.

Conclusions
In this paper, an LFC method for a pumped storage unit based on LADRC is proposed. Firstly, the LFC of reheat steam turbine with nonlinear factors is modeled, and according to the working mode of the pumped storage power station, the power plant LFC model under generating and pumping operations is established. The following conclusions are obtained from the simulation When load disturbance occurs in the pumping operation, the pumped storage power station can stop pumping after receiving the signal from the grid, cut off part of the load quickly, and support the power consumed by pumping to the grid to maintain the stability of the grid frequency. From the simulation results, it can be seen that the pumped storage power station working under the pumping operation obtains certain power support in a short time after the grid disturbance, and the frequency fluctuation range is shortened and restored to the stable state more quickly and smoothly. For the exchange power of the tie-line, the pumping operation reduces the value of the exchange power, significantly reduces the number of oscillations and amplitude. The recovery rate is more than twice that of without pumped storage energy, which is an effective supplement to the LFC of the power grid.

Conclusions
In this paper, an LFC method for a pumped storage unit based on LADRC is proposed. Firstly, the LFC of reheat steam turbine with nonlinear factors is modeled, and according to the working mode of the pumped storage power station, the power plant LFC model under generating and pumping operations is established. The following conclusions are obtained from the simulation results: (1) After the regional combination, the anti-interference ability of the system is strengthened, and power support can be obtained from the uninterrupted area of the interconnection, so as to better maintain the safe and stable operation of the power grid. (2) The LADRC method proposed in this paper can achieve better control effect and better robustness than FOPID and traditional PID control when controlling pumped storage power stations with nonlinear links. (3) It is effective and feasible to introduce the pumped storage power station into the LFC, which greatly improves the dynamic characteristics of the LFC of the power grid.