Control Strategy of Flywheel Energy Storage System Based on Primary Frequency Modulation of Wind Power

: As a form of energy storage with high power and efﬁciency, a ﬂywheel energy storage system performs well in the primary frequency modulation of a power grid. In this study, a three-phase permanent magnet synchronous motor was used as the drive motor of the system, and a simulation study on the control strategy of a ﬂywheel energy storage system was conducted based on the primary frequency modulation of wind power. The speed and current double closed-loop control strategy was used in the system start-up phase, and the power and current double-closed-loop control strategy were used in the power compensation phase. The model reference adaptive control was used to accurately estimate the speed and position of the rotor. The system compensates for the wind power output by using a wind turbine in real-time and conducting simulation experiments to verify the feasibility of the charge and discharge control strategy. At the same time, it can be veriﬁed that the ﬂywheel energy storage system has a beneﬁcial effect on wind power frequency modulation.


Introduction
Under General Secretary Xi Jinping's important instructions to reach peak carbon dioxide emissions by 2030 and achieve carbon neutrality before 2060, many provinces in China have successively introduced new energy configuration storage plans. New energy and energy storage projects are rapidly spreading across the country. As of October 2021, China's cumulative installed capacity of renewable energy power generation exceeded 1 billion kilowatts, accounting for 43.5% of the country's total installed power generation capacity. In the future, new energy will maintain rapid development momentum. It is estimated that the installed capacity of wind power and solar power will reach 1.2 billion kilowatts or more in 2030, which will surpass coal power and become the main portion of installed capacity [1]. Whether it is the large-scale and intensive development of centralized new energy or the nearby consumption of distributed new energy, the support of energy storage technology is an inseparable component.
The flywheel energy storage system (FESS) has a large capacity, high energy conversion rate, high instantaneous power, and high-frequency charge and discharge characteristics. It has broad application prospects in grid frequency modulation, uninterrupted power supply, and kinetic energy recovery and reuse. At the same time, the FESS belongs to mechanical energy storage, which can fit perfectly with new energy configuration energy storage. Although most domestic FESSs are only commercial demonstrations, there will be significant room for improvement in the near future.
Many adaptive control methods have been studied, improved, and applied. For example, Ref. [2] presents a characteristic model-based all-coefficient adaptive control law on an experimental test rig for high-speed FESS suspended on active magnetic bearings.
Ref. [3] presents an adaptive controller to provide failure compensation concerning a coaxial rotor helicopter to address the problem of instability due to rotor damage in a coaxial rotor helicopter. Ref. [4] presents an augmented subcutaneous model of type 1 diabetic patients by estimating the model parameters using the nonlinear least square method using the physiological data. Ref. [5] presents the coordinated control of the flywheel matrix system, which presents a neural-adaptive PID control algorithm to control each flywheel individually. Ref. [6] presents an immersion and invariance manifold adaptive controller to keep the constant DC-link voltage of FESS discharge in a wide speed range. Ref. [7] presents a controller, which is designed for the bidirectional Buck-Boost converter control in the FESS to improve the performance of FESS. Ref. [8] presents a stable, intelligent adaptive controller based on a recurrent fuzzy neural network for an induction machine-driven flywheel system. Ref. [9] presents an adaptive cut-off frequency for systems with multiple energy storage system units to realize several objectives simultaneously, including the stateof-charge limiting. Ref. [10] presents a novel adaptive control-based strategy for accurate simulation of solar panels under any given physical condition and for any general static or dynamic load. Ref. [11] presents an adaptive integral backstepping control scheme for tracking control of a 2-Degree of Freedom helicopter, which gives robustness to parameter uncertainties and unmodeled dynamics. Ref. [12] presents a model-free adaptive powerlevel control law for a nuclear superheated-steam supply system with forced primary circulation, which is free from physical and thermal-hydraulic parameters.
In this study, the FESS was used for the primary frequency modulation of the power grid. A three-phase permanent magnet synchronous motor (PMSM) was used as the drive motor, and a single FESS was used to compensate for wind power output by a wind turbine in real-time. The corresponding charging and discharging control strategy was proposed, and a simulation experiment of the FESS for wind power compensation was carried out.

The Overall Structure of the System
The FESS is usually composed of a drive motor, flywheel rotor, support bearing, and control system. As the core part of the system, the drive motor is a two-way motor: when the FESS is charging, the motor acts as an electric motor to drive the flywheel rotor to rotate; when the FESS is discharged, the motor acts as a generator for external power supply. The utilization rate of permanent magnet synchronous motors is the highest in the FESS, mainly because of its low rotor loss, high speed, and high operating efficiency. The material of the flywheel rotor has always been the focus of research, and the common ones are aluminum alloy, steel, glass fiber, etc. Currently, the support bearing part primarily uses magnetic suspension bearings in the high-speed field, such as high-temperature superconducting magnetic suspension bearings. Moreover, a good control system determines important factors, such as the rapid response and stability of the FESS [13][14][15]. The stored energy of a single FESS is directly determined by the moment of inertia and speed of the flywheel rotor, as shown in Equation (1): Usually, ω max = 2ω min , So available: where J is the moment of inertia of the flywheel rotor (in kg·m 2 ), ω is the speed of the flywheel rotor (in rad/s). Although the FESS was originally developed for kinetic energy recovery and reuse, its advantages in grid frequency modulation have gradually been developed in recent years. Figure 1 shows the structure diagram of the FESS used for the primary frequency regulation of wind power. It can be seen that the output power of the wind turbine and the FESS forms a complementary mode and finally discharges stable power to the grid. years. Figure 1 shows the structure diagram of the FESS used for the primary frequency regulation of wind power. It can be seen that the output power of the wind turbine and the FESS forms a complementary mode and finally discharges stable power to the grid.

Modeling of Pmsm
For a three-phase PMSM, the relevant variables of the motor are usually directly used to establish a mathematical model. The following assumptions were made: 1. The current of the motor is a symmetrical three-phase sine wave. 2. The magnetic saturation of the motor core was ignored.
3. The eddy current and hysteresis loss in the motor was ignored.
Then, the mathematical model is established under the three-phase static coordinate system [16,17], and the three-phase voltage is: where R is the equivalent resistance of each stator winding of PMSM (in Ω). The flux equation in the three-phase stationary coordinate system is: where L is the self-inductance between the magnetic chains (in H), M is the mutual inductance between the magnetic chains (in H), and is the spatial electrical angle of the angle between the N pole and the A-axis (in rad).

Modeling of Pmsm
For a three-phase PMSM, the relevant variables of the motor are usually directly used to establish a mathematical model. The following assumptions were made: 1.
The current of the motor is a symmetrical three-phase sine wave.

2.
The magnetic saturation of the motor core was ignored. 3.
The eddy current and hysteresis loss in the motor was ignored.
Then, the mathematical model is established under the three-phase static coordinate system [16,17], and the three-phase voltage is: where R is the equivalent resistance of each stator winding of PMSM (in Ω). The flux equation in the three-phase stationary coordinate system is: where L is the self-inductance between the magnetic chains (in H), M is the mutual inductance between the magnetic chains (in H), and θ e is the spatial electrical angle of the angle between the N pole and the A-axis (in rad).
The electromagnetic torque equation under the three-phase static coordinate system is: It can be seen from the above formula that the mathematical models under the threephase stationary coordinate system are very complicated, so they are usually transformed into mathematical models under the two-phase rotating coordinate system, as shown in the following Equation: The electromagnetic torque equation is: The Equation of motion is: where P is the differential operator, ω e is the electrical angular velocity (in rad/s), ω m is the mechanical angular velocity (in rad/s), n p is the number of pole pairs, L d and L q are the inductances (in H), ϕ f is the permanent magnet link (in Wb), T e is the electromagnetic torque (in N·m), T L is the no-load torque (in N·m), and B is the damping coefficient.

Modeling of Position Sensor-Less Control
Owing to the special design of the flywheel rotor in the FESS, it is impossible to install a position sensor in the PMSM. Therefore, position sensor-less technology was used to estimate the speed and position of the flywheel rotor. This study uses a model reference adaptive control (MRAS) to estimate the speed and position of the flywheel rotor.
The MRAS used in FESS mainly comprises a reference model, an adjustable model, and an adaptive law. The basic idea of its estimation is to use the actual motor model as a reference, input the d-axis and q-axis voltage and motor current output into the adjustable model, and feedback the speed calculated by the adaptive law to the adjustable model to reduce error. Therefore, the speed and position of the rotor can be accurately estimated [18][19][20]. A specific working principle diagram is presented in Figure 2. The MRAS used in FESS mainly comprises a reference model, an adjustable model, and an adaptive law. The basic idea of its estimation is to use the actual motor model as a reference, input the d-axis and q-axis voltage and motor current output into the adjustable model, and feedback the speed calculated by the adaptive law to the adjustable model to reduce error. Therefore, the speed and position of the rotor can be accurately estimated [18][19][20]. A specific working principle diagram is presented in Figure 2. In this diagram, means the variable's real value. Among them, for the linear forward channel, it is only necessary to satisfy the zero and pole of the transfer function in the left half plane of the S. For the nonlinear feedback channel, it is necessary to satisfy the POPOV superstability theorem. According to the working principle diagram of the MRAS, its control structure diagram can be obtained in Figure 3, and in this circumstance, the stability of the system model can be further calculated.

Adjustable model
Adaptive law  In this diagram, Among them, for the linear forward channel, it is only necessary to satisfy the zero and pole of the transfer function in the left half plane of the S. For the nonlinear feedback channel, it is necessary to satisfy the POPOV superstability theorem. According to the working principle diagram of the MRAS, its control structure diagram can be obtained in Figure 3, and in this circumstance, the stability of the system model can be further calculated.
The MRAS used in FESS mainly comprises a reference model, an adjustable model, and an adaptive law. The basic idea of its estimation is to use the actual motor model as a reference, input the d-axis and q-axis voltage and motor current output into the adjustable model, and feedback the speed calculated by the adaptive law to the adjustable model to reduce error. Therefore, the speed and position of the rotor can be accurately estimated [18][19][20]. A specific working principle diagram is presented in Figure 2. In this diagram, means the variable's real value. Among them, for the linear forward channel, it is only necessary to satisfy the zero and pole of the transfer function in the left half plane of the S. For the nonlinear feedback channel, it is necessary to satisfy the POPOV superstability theorem. According to the working principle diagram of the MRAS, its control structure diagram can be obtained in Figure 3, and in this circumstance, the stability of the system model can be further calculated.

Adjustable model
Adaptive law  where the i abc and u abc are obtained from the PMSM. Meanwhile, the current and voltage in the three-phase stationary coordinate system are converted into the current and voltage in the two-phase rotating coordinate system through the '3s/2r' module. The reference model is the model of the drive motor itself (refer to Equations (9)-(12)), and its voltage equation is changed to the following form: And assuming a new set of variables as: The equation of the adjustable model with angular frequency as the adjustable parameter is obtained as follows: The above operation is repeated to get the final adjustable model: According to the POPOV superstability theorem, to stabilize the error system, the following two conditions must be met simultaneously: 1.
The transfer function of the system, which is The nonlinear time-invariant link must satisfy the POPOV integral inequality.
where ε is a limited number of the normal.
After bringing C and D into G(s), it is easy to conclude that G(s) is a positive definite. As for the POPOV integral inequality, assuming f 1 , f 2 and h(t) as: And considering the adaptive law of the motor speed as: Substituting Equation (19) into Equation (17) yields: Split the Equation (20) , it can be deduced that To conclude, µ ≥ −ε 2 . The estimated rotor speed can be expressed as: where k p and k i are the proportional and integral coefficients, respectively. The estimated rotor position can be expressed as:

Modeling of Power Compensation Control Strategy
In this study, a hidden-pole three-phase PMSM and a vector control strategy with i d = 0 were used. With this method, only the quadrature axis current component in the stator current needs to be controlled to realize the charge and discharge control of the FESS.
In the FESS start-up phase, the control system uses the control strategy of the outer speed and inner current loops, and the rotor speed quickly increases to the charging and discharging starting speed (3000 rpm). After the system starts, the speed of the flywheel rotor remains unchanged when entering the power compensation phase.
With the FESS in the power compensation phase, the control strategy of the power outer and current inner loops was used to perform real-time peak shaving for the power curve output by a wind turbine. When wind power is too great, the FESS uses the excess power to charge up; contrarily, the FESS compensates for fluctuating and/or extremely low wind power by discharging power more smoothly to the grid [21][22][23][24][25]. The FESS flow chart for the real-time compensation of wind power output by a wind turbine is shown in Figure 4. During the entire process, the FESS relies on the characteristics of rapid charging and discharging to achieve rapid real-time compensation for wind power. The FESS control block diagram based on wind power frequency modulation is shown in Figure 5.    In the power outer loop control in Figure 5, the adjusted smooth wind power under ideal conditions is simulated according to the wind power output by a wind turbine, and the output power that needs to be compensated is obtained after the difference between the two. Calculate the specific value of according to = • .

• •
, and then input it to the current inner loop. At the same time, to ensure that the flywheel works within the safe range (3000-6000 rpm), the system is guaranteed to work within the safe range by limiting the input to the current inner loop.

Power Compensation Simulation
In the simulation, the actual wind power output of a wind turbine was used for compensation. The specific value of the wind power was measured every minute, and the 51min variable wind power was used as the wind power that requires compensation. In the simulation, the speed of the flywheel during the power compensation phase was 3000-6000 rpm, the moment of inertia of the flywheel rotor was 72.95 kg·m 2 , the self-dispersion rate of the system was 2% of the rated power, and the single-rated power of the FESS was 26 kW. The stand-alone energy storage capacity is 3.6 MJ, and the discharge depth is 75%. In the power outer loop control in Figure 5, the adjusted smooth wind power under ideal conditions is simulated according to the wind power output by a wind turbine, and the output power that needs to be compensated is obtained after the difference between the two. Calculate the specific value of i q according to i q = 9550·P 1.5n p ·ϕ f ·ω n , and then input it to the current inner loop. At the same time, to ensure that the flywheel works within the safe range (3000-6000 rpm), the system is guaranteed to work within the safe range by limiting the i q input to the current inner loop.

Power Compensation Simulation
In the simulation, the actual wind power output of a wind turbine was used for compensation. The specific value of the wind power was measured every minute, and the 51-min variable wind power was used as the wind power that requires compensation. In the simulation, the speed of the flywheel during the power compensation phase was 3000-6000 rpm, the moment of inertia of the flywheel rotor was 72.95 kg·m 2 , the selfdispersion rate of the system was 2% of the rated power, and the single-rated power of the FESS was 26 kW. The stand-alone energy storage capacity is 3.6 MJ, and the discharge depth is 75%. The system completes the start-up phase after 555 s, keeping the original speed unchanged, and entering the power compensation phase at 600 s, with a simulation time of 3600 s. The relevant parameters of the three-phase PMSM used are shown in Table 1. The wind power curve output by a wind turbine is shown in Figure 6, and the simulation model is shown in Figure 7 [26][27][28].        In Figure 7, the start-up phase is the maximum speed of 3000 rpm, and i q is obtained through the speed regulator (module ASR). In the power compensation stage, the module 'Wind Power' is the wind power output by the wind turbine, and the module 'Follow Power' is the adjusted power obtained under ideal conditions. The current value required to compensate the power is calculated from the rotational speedω n estimated by the module 'MRAS', and the current is limited to ensure that the work is within a safe range. i d and i q get u d and u q through the current regulator (module ACR). The module '2r-2s' converts the variables in the two-phase rotating coordinate system into the variables in the two-phase stationary coordinate system. The module 'SVPWM' is used to output the pulse signal to the power electronic converter to control the operation of the PMSM. Modules '3s-2r' are used to convert the variables in the three-phase stationary coordinate system into the variables in the two-phase rotating coordinate system. In the model, the u abc and i abc output by the PMSM are converted into u d , u q and i d , i q , which are further given to the module 'MRAS' to estimate the PMSM's rotational speed, thereby obtaining the estimated rotational speedω n and the estimated rotor positionθ.
The experiment realized compensations as high as 51 times/hour and effectively compensated the fluctuating wind power curve into a relatively smooth curve. Figure 8 shows the power curve after FESS compensation. Compared with the wind power output by a wind turbine, the wind power compensated by FESS is approximately 703 kW, and the overall power fluctuation range does not exceed 3%. Figure 9 shows a waveform diagram of the flywheel rotor speed. During the compensation process, the flywheel rotor speed was maintained at 3000-6000 rpm.
The experiment realized compensations as high as 51 times/hour and effectively compensated the fluctuating wind power curve into a relatively smooth curve. Figure 8 shows the power curve after FESS compensation. Compared with the wind power output by a wind turbine, the wind power compensated by FESS is approximately 703 kW, and the overall power fluctuation range does not exceed 3%. Figure 9 shows a waveform diagram of the flywheel rotor speed. During the compensation process, the flywheel rotor speed was maintained at 3000-6000 rpm.  The experiment realized compensations as high as 51 times/hour and effectively compensated the fluctuating wind power curve into a relatively smooth curve. Figure 8 shows the power curve after FESS compensation. Compared with the wind power output by a wind turbine, the wind power compensated by FESS is approximately 703 kW, and the overall power fluctuation range does not exceed 3%. Figure 9 shows a waveform diagram of the flywheel rotor speed. During the compensation process, the flywheel rotor speed was maintained at 3000-6000 rpm.   Figure 10 shows an enlarged part of the speed waveform. In the figure, it can be seen that the error between the estimated speed and the actual speed does not exceed 0.5%. Therefore, the results using the MRAS to estimate the speed are reliable.  Figure 10 shows an enlarged part of the speed waveform. In the figure, it can be seen that the error between the estimated speed and the actual speed does not exceed 0.5%. Therefore, the results using the MRAS to estimate the speed are reliable.

3158.7
The Actual Speed The Estimated Speed Figure 10. The enlarged part of the speed waveform. Figure 11 shows the electromagnetic torque waveform of the PMSM. It can be seen from the figure that the general trend of the electromagnetic torque is basically consistent with the wind power curve that needs adjustment.  Figure 11 shows the electromagnetic torque waveform of the PMSM. It can be seen from the figure that the general trend of the electromagnetic torque is basically consistent with the wind power curve that needs adjustment.  Figure 10. The enlarged part of the speed waveform. Figure 11 shows the electromagnetic torque waveform of the PMSM. It can be seen from the figure that the general trend of the electromagnetic torque is basically consistent with the wind power curve that needs adjustment.  Figure 12 shows the three-phase current waveform, and the change in the current size is basically the same as the change in the size of the electromagnetic torque. Figure 13 shows the enlarged part of the current waveform, and it can be seen that the current fluctuation is very small, the period is approximately 0.51 s, and the initial phase angle is 0 rad. The amplitude of the current in the start-up phase is approximately 62.72 A. In the power compensation phase, the amplitude of the three-phase alternating current changes with the torque.   Figure 12 shows the three-phase current waveform, and the change in the current size is basically the same as the change in the size of the electromagnetic torque. Figure 13 shows the enlarged part of the current waveform, and it can be seen that the current fluctuation is very small, the period is approximately 0.51 s, and the initial phase angle is 0 rad. The amplitude of the current in the start-up phase is approximately 62.72 A. In the power compensation phase, the amplitude of the three-phase alternating current changes with the torque.

Time (s)
Time (s) Figure 10. The enlarged part of the speed waveform. Figure 11 shows the electromagnetic torque waveform of the PMSM. It can be seen from the figure that the general trend of the electromagnetic torque is basically consistent with the wind power curve that needs adjustment.  Figure 12 shows the three-phase current waveform, and the change in the current size is basically the same as the change in the size of the electromagnetic torque. Figure 13 shows the enlarged part of the current waveform, and it can be seen that the current fluctuation is very small, the period is approximately 0.51 s, and the initial phase angle is 0 rad. The amplitude of the current in the start-up phase is approximately 62.72 A. In the power compensation phase, the amplitude of the three-phase alternating current changes with the torque.

Economic, Social and Ecological Benefits
Reviewing the development history of domestic wind power, Figure 14 shows the new domestic wind power hoisting situation [29]. The development of domestic wind power can be divided into two major development stages: The first stage occurred in 2010 and earlier, a typical growth period. At this stage, the domestic wind power and wind power industry showed rapid development, driven by policies in the early stage of the industry's development.

Economic, Social and Ecological Benefits
Reviewing the development history of domestic wind power, Figure 14 shows the new domestic wind power hoisting situation [29]. The development of domestic wind power can be divided into two major development stages:

Economic, Social and Ecological Benefits
Reviewing the development history of domestic wind power, Figure 14 shows the new domestic wind power hoisting situation [29]. The development of domestic wind power can be divided into two major development stages: The first stage occurred in 2010 and earlier, a typical growth period. At this stage, the domestic wind power and wind power industry showed rapid development, driven by policies in the early stage of the industry's development.
The second stage occurred from 2011 to 2020 when the demand showed obvious cyclical fluctuations. From 2011 to 2020, policies dominated the development of the country's wind power industry. Among them, favorable and unfavorable policies are intertwined, resulting in obvious cyclical fluctuations in industry demand.  Summarizing the reasons behind the growth and decline of the domestic wind power industry in the past ten years, it is mainly due to the high level of domestic wind curtailment rate and certain fluctuations. It was not until 2019 that the domestic wind curtailment rate fell to less than 5%. The high level of wind curtailment rate, on the one hand, inhibits the developers' enthusiasm for wind power investment; on the other hand, it urges regulators to strengthen supervision to solve the problem of wind curtailment. One of the main means is to control the supply of new wind power projects, which will negatively affect the newly installed capacity.
According to the national installed capacity of 400 million kW of new energy, flywheel energy storage with 3% of the installed capacity, and the cost of 2.5 yuan per watt after industrialization, the national market size will reach 30 billion yuan, the installed Summarizing the reasons behind the growth and decline of the domestic wind power industry in the past ten years, it is mainly due to the high level of domestic wind curtailment rate and certain fluctuations. It was not until 2019 that the domestic wind curtailment rate fell to less than 5%. The high level of wind curtailment rate, on the one hand, inhibits the developers' enthusiasm for wind power investment; on the other hand, it urges regulators to strengthen supervision to solve the problem of wind curtailment. One of the main means is to control the supply of new wind power projects, which will negatively affect the newly installed capacity.
According to the national installed capacity of 400 million kW of new energy, flywheel energy storage with 3% of the installed capacity, and the cost of 2.5 yuan per watt after industrialization, the national market size will reach 30 billion yuan, the installed capacity of new energy will reach 800 million kW, and the market size of flywheel energy storage will reach 60 billion yuan.
In 2020, the national curtailment of wind power was 16.61 billion kWh, and the wind power utilization rate was 96.5%, an increase of 0.5% annually, the abandonment of photovoltaic power was 5.26 billion kWh, and the photovoltaic power generation utilization rate was 98.0%, basically the same as in the previous year.
In the fourth quarter, the level of consumption and utilization of new energy across the country was relatively high, but the utilization rate of new energy declined annually. In the fourth quarter, the national abandonment of wind power was 4.97 billion kWh, an annual increase of 23.3%, and the utilization rate of wind power was 96.3%, an annual decrease of 0.3%. The utilization rate of PV power generation was 97.0%, with an annual decrease of 0.6%, and the level of new energy consumption and utilization was relatively high.
In recent years, owing to the large number of new energy installations, the reduction of the power grid has not kept up with the development of new energy, so the rate of wind curtailment has been rising. In this case, although there have been a series of measures implemented since 2018 to solve the wind curtailment problem, such as 'internal expansion and external transmission', the existence of new forms of energy storage is also necessary. From the experimental results, it can be concluded that FESS, as a new form of energy storage, will have a significant effect on reducing the wind curtailment rate.
At present, the FESS is still in the development stage, with a domestic market share of less than 1%. However, as flywheel energy storage technology matures and costs decrease in the future, its role and economic benefits will become more prominent. Future scenarios include:

1.
Promotion of large-scale price reductions of similar imported products and related technical equipment. According to the current scale of the flywheel energy storage market, this can save hundreds of millions of yuan each year. 2.
Applying the FESS array parallel technology in the power grid can improve the security and power supply quality. It is expected that the economic benefits of the FESS array on the power supply quality of the power grid during the entire life cycle can reach hundreds of millions of yuan per year. 3.
FESS can complete hundreds of thousands of charge and discharge cycles within its full life cycle, covering an area of only 50% of the traditional battery energy storage system, without discharging any toxic or harmful chemical materials, and is environmentally friendly.

Conclusions and Future Work
In this study, the FESS was used to compensate for the wind power discharge by a wind turbine in real-time, and a three-phase PMSM was used as the drive motor. The system used a double closed-loop control of speed and current in the start-up phase and the double closed-loop control of power and current in the power compensation phase. MRAS was used to estimate the speed of the flywheel rotor. The simulation effectively adjusted the wind power in real-time and verified the feasibility of using the FESS for wind power frequency modulation. The results show that applying the FESS to primary frequency modulation, such as wind and solar, has positive consequences. In the future, when demand for renewable energy power generation such as wind and solar greatly increases, using an array composed of multiple FESSs to perform primary frequency modulation on the actual output power of wind farms will prove reliable. The FESS plays a significant role in the frequency adjustment of the power grid.