**1. Introduction**

Due to increasing environmental concerns, wind power generation has undergone rapid growth [1]. In 2013, the world wind power production capacity was 318GW, and it is forecasted to reach 712 GW, 1480GW, 2089GW, and 2672GW by 2020, 2030, 2040, and 2050, respectively, in moderate scenario [2].

In China, the figure is 123GW, 216GW, 414GW, 680GW, and 1000GW in 2015, 2020, 2030, 2040, and 2050, respectively [2, 3].

Wind energy production brings immense advantages such as sustainable, incredible domestic potential, eco-friendly, revitalizing to rural economies, low operational cost, space efficient, etc. Nevertheless, its main drawback is its stochastic over time in nature; hence, significant power fluctuations are observed in the wind farm. For a power system even with moderate wind power penetration, the fluctuations should be mitigated; otherwise, this may lead to substantial deviations in the grid frequency [1], voltage sag and flicker at the grid busses [4], steady-state voltage deviation, even equipment damage, and system collapse at large. A study in busses [4] asserts that the variable-speed pumped storage (VSPS) system is able to improve the dynamic stability and steady-state operations of the power system. Having dependable features and performances, the VSPS can play a crucial role in stability control and frequency regulator and AC voltage control during contingencies.

Previous works have focused on converter topology development [5–11], which has led to improvements in VSPS system performance and penetration of renewable energy development level [12]. Furthermore, control strategies have been intensively studied in doubly fed induction machine (DFIM)-based VSPS applications. These works mainly relate to the vector control category, including Field-Oriented Control (FOC) [9, 13, 14] and Direct Torque Control (DTC) [5, 8, 11, 13, 14]. As a result, the system performance improvement has achieved remarkable results.

In Ref. [8], a direct power control (DPC)-based frequency regulator is implemented in a full-scale converter-fed synchronous machine. DPC is also used in the application of DFIM-based VSPS systems in [11], with a focus on converter topology studies. A grid-integrated VSPS system was proposed in [5] to study the impact of faults on grid frequency stability and VSPS active power flow. In [9], a VSPS-based H-bridge cascaded multilevel converter with a stator voltage FOC strategy was proposed to suppress the effects of wind farm power fluctuations. On the other hand, the work in [1] proposed a power filter algorithm to solve the control method to adjust the frequency deviation of the power grid caused by wind power fluctuations. However, the above work does not consider the phasor model simulation technique because it is very important in the stability and control analysis of large power systems where simulation time and computational storage are critical. Since sinusoidal voltages and currents are replaced by phasors in polar form in the phasor model technique; and where electromagnetic transients are not important, dynamic simulation time is greatly reduced. Implementing a phasor model in any linear system is also an advantage, since the study of small-signal stability is based on the eigenvalue analysis of linearized power systems, which complements the dynamic simulation of nonlinear systems [15]. However, implementing phasor modeling techniques in DFIM applications, grid frequency control has always been a challenge. In the phasor model technique, the grid frequency is assumed to be a constant fundamental frequency. It is impossible to develop a phasor locked loop for synthesizing frequencies. Therefore, grid frequency control is performed through an active power control strategy, which is an open-loop control scheme for frequency, as shown in **Figure 1**. Therefore, the response of grid frequency and AC voltage lacks good dynamic performance, especially during emergencies [16].

On the other hand, grid frequency is very sensitive to load/source changes. If the load/power connected to the AC grid system is suddenly turned on/off, the grid frequency will decrease/increase, causing the power plant to generate more power or stop it. Furthermore, due to the applied frequency drop, the converter will respond to the decrease in frequency by increasing the inverter power or decreasing the rectifier power. Therefore, the frequency drop is estimated per unit increment of rectified

*Wind Power Fluctuation Compensation by Variable-Speed Pumped Storage Plant in a Grid… DOI: http://dx.doi.org/10.5772/intechopen.104938*

#### **Figure 1.**

*Open-loop frequency and AC-bus voltage control through the power control strategy.*

power fed to the closed control system. Since a change in the converter power results in a change in the DC voltage level, the frequency drop will appear as a drop in the DC voltage level in the AC/DC interface system, which in turn results in the need for power flow compensation [17]. Therefore, the way to solve this problem is that adding a synchronous machine to the grid system helps to obtain the measured and estimated frequency, because the rotor speed of the synchronous machine is absolutely the same as the scale factor of the grid frequency.

However, this chapter aims to propose a frequency droop feeding DPC for VSPS systems to accommodate grid frequency deviations due to wind power fluctuation effects. Spectral analysis is newly used to tune the stability balance and verify the dynamic performance of the proposed control system. The converter used is a threelevel back-to-back neutral-point clamped voltage source converter (NPC-VSC).

### **2. Basic description of frequency control**

The replacement system frequency is a measure of the effective power balance in the system. It is a feedback signal that allows the operator to monitor the balance of energy produced and consumed. It is a common characteristic of all location voltages, but varies over time, while system conditions change within an AC power system. Any change in power generation or load will result in a frequency deviation that can lead to system instability and minimal operation of generators and consumers. In extreme cases, system equipment may be damaged.

The concept of power system frequency stability is defined as the requirement to ensure that the system frequency remains within a given safe range during normal operation and emergency situations [1]. This is ensured by using two mechanisms, frequency control and active power balance control [18]. Frequency control is a crucial control problem in system design and operation. This is even more important today due to the complexity of modern power systems driven by ever-increasing scale, environmental constraints, emerging renewable energy sources, and uncertainty.

Based on the dynamic model of generator load based on the correlation between mismatched power (ΔPm - ΔPL) and grid frequency deviation, a simplified frequency response model can be defined for an interconnected multigenerator power system. (Δf) and can be expressed by the oscillatory differential equation [19] as

$$
\Delta P\_m(t) - \Delta P\_L(t) = 2H \frac{d\Delta f(t)}{dt} + D\Delta f(t) \tag{1}
$$

where ΔPm, ΔPL, D, and H are respectively change of the mechanical power, change of the load, the load damping coefficient, and the inertia constant.

In terms of Laplace transformation function, the expression of the change of frequency in (1) is denoted by

$$
\Delta f(s) = \frac{\Delta P\_m(s) - \Delta P\_L(s)}{2Hs + D} \tag{2}
$$

It can be seen from Eq. (2) that the frequency deviation is inversely proportional to the load damping coefficient and inertia constant.

To keep the power system frequency at the nominal value, different mechanisms can be used depending on the magnitude of the frequency deviation faced. Small frequency deviations can be adjusted by primary control and load frequency control systems. When the situation is complex and the deviation is large, emergency control schemes such as low-frequency load shedding should be used to restore the system frequency.

Frequency stability studies can range in time from seconds to minutes as both fast and slow dynamics affect it. For high-power fluctuating systems, the long-term stability issue [19] is considered so in this chapter.

### **3. The wind power**

Wind speed is random in nature and varies erratically over time. As a result, significant wind fluctuations can be observed.

For a given power factor Cp, the mechanical power (Pm) and torque (Tm) produced by the wind turbine rotor can be defined as

$$P\_m = \mathbf{0}.\mathsf{5}\rho A C\_p(\lambda, \boldsymbol{\beta}) v\_w \,^3 \tag{3}$$

$$T\_m = \frac{P\_m}{\alpha\_r} \tag{4}$$

where ρ, vw, A, and ω<sup>r</sup> are respectively the air density, wind speed, swept area, and angular speed of the generator rotor. Maximum mechanical power extraction is possibly made at maximum value Cp and is a function of the pitch angle (β) and tip-speed ratio (λ), which are related as

$$
\lambda = \frac{\alpha\_{rot} R\_T}{v\_w} \tag{5}
$$

The angular speed of the wind turbine ωrot is related to the generator rotor angular speed ω<sup>r</sup> and as a function of gear ratio Ng and given by

$$
\rho\_r = \mathcal{N}\_\mathcal{g} \rho\_{rot} \tag{6}
$$

The power coefficient expression Cp has been shown to have a unique maximum, Cp-max, for λopt = 8.1 and βopt = 0 [20].

The maximum coefficient of Cp, i.e., the Betz limit is 16/27. But maximum value of Cp is determined in this chapter. Based on the modeling turbine characteristics of [21], a generic Eq. (7) is used to model Cp(λ,β).

$$\mathbf{C}\_{p}(\boldsymbol{\lambda},\boldsymbol{\beta}) = \mathbf{C}\_{1}(\mathbf{C}\_{2}/\lambda\_{i} - \mathbf{C}\_{3}\boldsymbol{\beta} - \mathbf{C}\_{4})\mathbf{e}^{-C\_{5}/\lambda\_{i}} + \mathbf{C}\_{6}\boldsymbol{\lambda}$$

$$\Longrightarrow \frac{\mathbf{1}}{\lambda\_{i}} = \frac{\mathbf{1}}{\lambda + \mathbf{0.08\beta}} - \frac{\mathbf{0.035}}{\boldsymbol{\beta}^{3} + \mathbf{1}}\tag{7}$$

*Wind Power Fluctuation Compensation by Variable-Speed Pumped Storage Plant in a Grid… DOI: http://dx.doi.org/10.5772/intechopen.104938*

**Figure 2.**

*Wind turbine power versus speed characteristics.*

Particular optimal value of λ is defined (λopt = 8.1). Thus, for β = 0<sup>0</sup> and the coefficients of C1 = 0.5176, C2 = 116, C3 = 0.4, C4 = 5, C5 = 21, and C6 = 0.0068, the maximum value of Cp, which is Cpmax = 0.48, is achieved.

Normalizing Eq. (3) in the pu system, we get

$$P\_{m-pu} = \mathsf{C}\_{p-pu}(\lambda, \beta) \upsilon\_{w-pu}^3 \tag{8}$$

Thus, implementing the above equations, the turbine power is determined and the result is shown in **Figure 2** as power versus speed of the wind turbine.

The electromagnetic torque (Tem) is expressed as

$$T\_{em} = p\left(\Psi\_{ds}i\_{qs} - \Psi\_{qs}i\_{ds}\right) \tag{9}$$

The output active (Pg) and reactive power (Qg) of a wind turbine supplied to the grid are computed by.

$$\begin{aligned} P\_{\mathcal{g}} &= \mathbf{1.5} \left( -\nu\_{ds} i\_{ds} + \nu\_{qs} i\_{qs} - \nu\_{qr} i\_{qr} - \nu\_{dr} i\_{dr} \right) \\ Q\_{\mathcal{g}} &= \mathbf{1.5} \left( \nu\_{qs} i\_{ds} + \nu\_{ds} i\_{qs} - \nu\_{dr} i\_{qr} - \nu\_{qr} i\_{dr} \right) \end{aligned} \tag{10}$$

### **4. Modeling of the system**

#### **4.1 System description**

The structure of the proposed system is shown in **Figure 3**. The grid contains conventional synchronous generator-based thermoelectricity, represented by SG1, SG2, and SG3, and loads. VSPS stands for variable-speed pumped storage powered by a doubly fed induction machine. As shown in **Figure 3**, the VSPS induction motor is configured with the rotor connected to the rotor-side voltage source converter (VSC) and the stator connected to the grid. The grid-side VSC NPCg is connected to the grid through a coupling inductor. Details are explained in [22]. IG stands for induction generator-based wind farm.

#### **4.2 The VSPS machine modeling**

As shown in **Figure 3**, the doubly fed induction machine (DFIM)-fed VSPS unit is based on the VSC NPC converter topology. The dynamics model of the three-phase induction machine is expressed in Eqs. (11)–(15) including electrical and mechanical systems. Further details are explained in [22].

$$\begin{aligned} \frac{d\Psi\_{ds}}{dt} &= V\_{ds} - R\_{s}i\_{ds} - a\_{s}\Psi\_{qs} \\ \frac{d\Psi\_{qs}}{dt} &= V\_{q\circ} - R\_{s}i\_{q\circ} + a\_{s}\Psi\_{ds} \end{aligned} \tag{11}$$

$$\begin{aligned} \frac{d\Psi\_{dr}}{dt} &= V\_{dr} - R\_r i\_{dr} + (\alpha\_s - \alpha\_r)\Psi\_{qr} \\ \frac{d\Psi\_{qr}}{dt} &= V\_{qr} - R\_r i\_{qr} - (\alpha\_s - \alpha\_r)\Psi\_{dr} \end{aligned} \tag{12}$$

$$T\_{em} = \mathbf{1.5p} \left( \Psi\_{ds} i\_{qs} - \Psi\_{qs} i\_{ds} \right) \tag{13}$$

*Wind Power Fluctuation Compensation by Variable-Speed Pumped Storage Plant in a Grid… DOI: http://dx.doi.org/10.5772/intechopen.104938*

$$\begin{aligned} \Psi\_{ds} &= L\_i i\_{ds} + L\_m i\_{dr} \\ \Psi\_{qr} &= L\_i i\_{qr} + L\_m i\_{qr} \\ \Psi\_{dr} &= L\_r i\_{dr} + L\_m i\_{ds} \\ \Psi\_{qr} &= L\_r i\_{qr} + L\_m i\_{qr} \\ L\_s &= L\_{ls} + L\_m \\ L\_r &= L\_{lr} + L\_m \\ \frac{d\phi\_m}{dt} &= \frac{1}{2H} \left( T\_{cm} - F\phi\_m - T\_m \right) \\ L\_0 &\tag{15} \end{aligned} \tag{16}$$

$$\frac{d\theta\_m}{dt} = o\_m \tag{15}$$

where, *Ψdqs Ψdqr*, *Vqdsidqs*, *LlsLlr*, *Vqdridqr*, *Lm*, *Rs Rr*, *ωωr*, *TmTem*, *H*,*F*, *s*, and *P* are respectively dq-axis stator and rotor fluxes, dq-axis stator voltage and current, stator and rotor leakage reactance, dq-axis rotor voltage and current, magnetizing reactance, stator and rotor resistance, synchronous and rotor angular speeds, turbine shaft and electromagnetic torque, combined inertia constant, combined viscous friction coefficient, slip ratio of the induction machine, and number of pole pairs.

#### **4.3 The VSC-NPC converter modeling for deploying VSPS**

The VSC is a self-commutating converter that is fast and controllable for AC/DC interface applications. One type of VSC converter is a three-level NPC power converter. The tertiary block contains three arms, each with its own four switch assemblies with antiparallel diodes and two NPCs.

The VSC self-commutated converter is a fast and controllable converter for AC/DC interface applications. One type of VSC converter is a three-level NPC power converter. The three-level block contains three arms with their own four switch assemblies with antiparallel diodes and two NPCs. Nonetheless, since this chapter is devoted to verifying the performance of the proposed system according to the phasor model technique, since the proposed system is large and should exhibit wind energy compensation and frequency stability not possible to see in the detailed model technique, assumptions are adopted. It is assumed that power losses in the converter are ignored. The switching dynamics can also be ignored because the PWM frequency in the threelevel VSC NPC is much larger than the system frequency [23]. Therefore, the model of VSC NPC was developed based on (16).

$$
\mu\_{dqr} = \frac{m\_{dqr}V\_{dc}}{2} \tag{16}
$$

**Figure 4.** *DC-link phasor model.*

where udqr, udqg and mdqr, mdqg are converter-injected voltages and modulated signals respectively in the rotor and grid-side converters.

The converter dynamics can be given by.

$$\frac{dV\_{dc}}{dt} = \frac{1}{C\_{eq}} \dot{\iota}\_c = \frac{1}{C\_{eq}V\_{dc}} \left(P\_g - P\_r\right) \tag{17}$$

Or in Laplace transform function

$$V\_{dc}(\mathbf{s}) = \frac{1}{\mathfrak{s}\mathbf{C}\_{eq}} I\_{\mathfrak{c}}(\mathbf{s}) = \frac{1}{\mathfrak{s}\mathbf{C}\_{eq}V\_{dc}(\mathbf{s})} \left[P\_{\mathfrak{g}}(\mathbf{s}) - P\_r(\mathbf{s})\right] \tag{18}$$

where Ceq and Vdc are the equivalent capacitance and DC voltage, respectively, and ic is the DC link capacitor charging current. Pg and Pr are active power flow in the grid and rotor side of the converters, respectively. Hence, Eq. (18) can be illustrated by **Figure 4** [22].

#### **4.4 Control strategies for the VSC NPC system**

#### *4.4.1 Active power control*

Using the DPC strategy, the response speed is fast, and the control system directly and effectively compensates for wind power fluctuations. Even if this is a local control scheme placed in the VSPS unit, as long as the VSPS operates within specified limits, an important output of an adjustable and almost constant power flow can be seen from the network of the power system. The power command Popt of the VSPS is determined by the capacity of the VSPS and its efficiency.

In this regard, an external setpoint Pset should be assumed for the VSPS, which is the input to the power control system. The set point depends on the VSPS capacity and grid conditions. The value of this setpoint is specifically used to optimize the system energy balance and is updated at a slower rate. However, the main focus is on controlling wind power fluctuations and regulating VSPS power to mitigate the impact on the grid.

The active power Ps and reactive power Qs equations provide the basis for the control strategy modeling, and given as.

$$\begin{aligned} P\_s &= \mathbf{1.5} (\upsilon\_{qs} i\_{qs} + \upsilon\_{ds} i\_{ds}) \\ Q\_s &= \mathbf{1.5} (\upsilon\_{qs} i\_{ds} - \upsilon\_{ds} i\_{qs}) \end{aligned} \tag{19}$$

Based on the detail works in [22], we have.

$$\begin{aligned} P\_s &= \mathbf{1.5V}\_s (L\_m/L\_s) i\_{qr} \\ Q\_s &= \mathbf{1.5V}\_s \{ (\Psi\_s/L\_s) - (L\_m/L\_s) i\_{dr} \} \end{aligned} \tag{20}$$

which implies the independent control of the real power and reactive power in the application of DFIM in VSPS.

Starting from the voltage equations of (11)–(14), assuming that the voltage drop across the stator resistance is very small compared with the grid voltage, and that the magnitude of the stator flux is fairly constant, we have the mathematical equation of (21). *Wind Power Fluctuation Compensation by Variable-Speed Pumped Storage Plant in a Grid… DOI: http://dx.doi.org/10.5772/intechopen.104938*

$$V\_{qr} = R\_r i\_{qr} + L\_r^\* \frac{di\_{qr}}{dt} + (\alpha\_s - \alpha\_r)(L\_r i\_{dr} + L\_m i\_{dr})$$

$$V\_{dr} = R\_r i\_{dr} + L\_r^\* \frac{di\_{dr}}{dt} - (\alpha\_s - \alpha\_r) \left(L\_r i\_{qr} + L\_m i\_{qs}\right)$$

where

$$L\_r^\* = L\_r - \left(L\_m^2 / L\_s\right) \tag{21}$$

From (21), the rotor dynamics model is determined and defined as

$$\frac{I\_{dqr}(\mathfrak{s})}{I\_{dqr-r\mathfrak{e}f}(\mathfrak{s})} = \mathbf{G}\_{\mathfrak{i}}(\mathfrak{s}) = \frac{\mathbf{1}}{\frac{L^{\ast}}{\mathbb{R}\_{r}}\mathfrak{s} + \mathbf{1}} \tag{22}$$

Thus, the current control (inner loops) structure of **Figure 5** is supposed to be developed. The PI control of *Ki(s)* is defined by

$$K\_i(\mathbf{s}) = \frac{P\_i \mathbf{s} + I\_i}{\mathbf{s}} \tag{23}$$

where *Pi* and *Ii* are proportional and integral gains. *Pi* and *Ii* are determined and tuned based on the control stability theory.

From Eq. (23), we obtain

$$I\_{qr}(\mathbf{s}) = G\_i(\mathbf{s}) I\_{qr - r \neq}(\mathbf{s}) \tag{24}$$

Multiplying both sides of (24) by *1.5VsLm/Ls*, we get

$$\mathbf{1.5V\_{s}L\_{m}/L\_{s}I\_{qr}(s) = G\_{i}(s)\mathbf{1.5V\_{s}L\_{m}/L\_{s}I\_{qr-r\not\!f}(s)}\tag{25}$$

Therefore, based on Eq. (20), from (25), we can deduce (26).

$$P\_s(\mathfrak{s}) = G\_i(\mathfrak{s}) P\_{s-r\!f}(\mathfrak{s}) \tag{26}$$

That is the dynamic model for active power control of rotor-side converter VSC NPCr is the same as the current control dynamic model. The PI control of Kp(s) of **Figure 6** is defined in (30)

$$K\_p(\mathfrak{s}) = \frac{P\_p \mathfrak{s} + I\_p}{\mathfrak{s}} \tag{27}$$

**Figure 5.** *The current control structure of the VSC NPC1.*

**Figure 6.** *The active power control structure of the VSC NPC1.*

Considering the power command ΔPwind to the VSPS power control unit, it is computed by

$$
\Delta P\_{wind} = P\_{wind} - P\_{set} \tag{28}
$$

where Pwind is the wind farm power measured, and Pset an assumed forecasted power output. This Pwind is added to the power command Popt to feed the fluctuated wind power to the control loop for regulating.

#### *4.4.2 Frequency droop control for the improvement of active power control performance*

Controlling the VSPS unit in primary frequency control is very important in power systems. By implementing a frequency droop control scheme, the response of the power system to emergencies is improved in the case of fluctuations in the dominant wind farm power generation. Power transfer within the DC link can be modulated by frequency droop control. It is specially designed for converters to provide frequency support to the grid. The control structure shown in **Figure 7** provides a governor-like droop behavior through active power flow. Droop-type control is built on the concept of power synchronization control, where grid synchronization is achieved regardless of dedicated synchronization units [24]. This droop control system provides a DC link with frequency droop characteristics. Load sharing is possible for other generator sets in the grid system.

The introduction of VSPS can support the improvement of system performance to cope with all emergencies of the power system. Therefore, the additional power

**Figure 7.** *Typical wind speed characteristics.*

*Wind Power Fluctuation Compensation by Variable-Speed Pumped Storage Plant in a Grid… DOI: http://dx.doi.org/10.5772/intechopen.104938*

**Figure 8.** *A frequency droop-fed DPC control strategy for VSPS system.*

command ΔPf due to the frequency variation of the VSPS active power control system (29) is determined.

$$
\Delta P\_f = -K\_{droop} \left( f\_{ref} - f\_{grid} \right) \tag{29}
$$

where Kdroop is the droop constant, fref is the 50 Hz reference frequency used in this chapter, and fgrid is the grid frequency. Through the active power closed-loop control system, frequency droop control can be used as a support for grid frequency control. Therefore, in order to obtain more possible controllability in the operation of the power system, the frequency droop control of the VSPS system must be introduced significantly. The above strategies are combined for the active power control scheme in the VSPS unit in the proposed system, as shown in **Figure 8**.

The models and details of the reactive power control in both rotor-side and grid-side converters and DC voltage control in the grid-side converter are presented in [22].

#### **4.5 The wind power fluctuation and its spectrum characteristics**

The distribution of wind speed and direction varies by region and season overtime. The spectral characteristics of wind power fluctuations recorded over time for a typical wind farm can be obtained using the Fast Fourier Transform (FFT) algorithm. **Figure 7** shows wind speed fluctuations; the corresponding wind power generation and its spectral characteristics are shown in **Figure 9a** and **b**, respectively.

Spectral analysis is the process of estimating the power spectral density (PSD), which characterizes the frequency content or random processes of a signal, from its time domain representation. Spectrum automatically decomposes a signal or random process into different frequencies and identifies periodicity. FFT-based nonparametric methods make no assumptions about the input data, can be used for any type of signal, produce more accurate spectral estimates, and allow us to convert flow timedomain signals to frequency domain and vice versa. The frequency domain representation of a signal reveals important signal characteristics that are difficult to analyze in the time domain.

**Figure 9.** *The wind-turbine output power; (a) in time domain, (b) the frequency spectrum characteristics.*

One of the definitions of power spectral density (PSD) function [25] is given by

$$\phi(\alpha) = \lim\_{N \to \infty} E\left\{ \frac{1}{N} \left| \sum\_{t=1}^{N} \mathcal{y}(t) e^{-i\alpha t} \right|^2 \right\} \tag{30}$$

Certainly, the PSD *ϕ*(*ω)* measures the power at frequency ω in the signal's autocovariance sequence *r*(*k*), i.e., *r*(*k*) = *r\**(�*k*). Since *ϕ*(*ω*) is a power density, it should be real valued and nonnegative, i.e., *ϕ*(*ω*) ≥ *0* for all *ω*. Thus, according to [25], the other equivalent definition of the PSD is given by

$$\phi(\alpha) = r(0) + 2\sum\_{k=1}^{\infty} r(k)\cos\left(\alpha k\right) \tag{31}$$

where *r*(*0*) is the averaging value of the function over the sampled time, and its inverse transformation is

$$r(k) = \frac{1}{2\pi} \int\_{-\pi}^{\pi} \phi(\alpha) e^{i\alpha k} d\alpha \tag{32}$$

Therefore, the spectral characteristic of the wind power fluctuation is well described and analyzed by the above PSD equations, and the simulating result of **Figure 9a** is shown in **Figure 9b**.

As shown in **Figure 9b**, the wind power fluctuation is insignificant for the spectral frequency greater than about 1.2 Hz. Hence, the VSPS is engaged to compensate the

*Wind Power Fluctuation Compensation by Variable-Speed Pumped Storage Plant in a Grid… DOI: http://dx.doi.org/10.5772/intechopen.104938*

significant wind power fluctuation components and the spectrum of the VSPS control response is determined accordingly.

According to statistics, the typical frequency of wind power fluctuations is recorded from a few minutes to a few days [9]. For slow fluctuations that are not an obvious problem, such as daily, the grid has adequate means of control. However, it is difficult to adjust minute-to-minute fluctuations only by traditional means such as thermal power. Therefore, to solve such problems caused by the randomness of wind power, VSPS is a possible option for stabilizing and smoothing grid frequency changes.

### **4.6 The VSPS system spectrum analysis and its implication for wind power fluctuation compensation**

The power control in the DFIM-based wind farm VSPS system is essentially the control of the VSC-fed VSPS system. As shown in **Figure 10**, the frequency-droop-fed DPC-controlled VSPS system integrated in the grid tracks the power fluctuation command well. Small distortions are observed when the command signal changes from one state to another.

We know that the dynamic characteristics of the generator set and pump set of VSPS are closely related to the control method and system structure, but it is difficult to obtain a specific value by using the analytical method. Therefore, this chapter uses the spectral analysis method to verify the results obtained by the time domain analysis, as shown in **Figure 11**. The spectral characteristics are determined using the driving mathematically equivalent total active power control system equation of (33). Therefore, the simulation results of the spectral characteristics of the controlled active power response of the VSPS are shown in **Figure 12**. It is observed that the bandwidth of the response of the VSPS control system is large enough to cover the bandwidth power fluctuation of the spectral results of the wind as shown in **Figure 10b**. This means that the frequency-fed DPC control strategy of the VSPS system can use the derivation function of (33) to

**Figure 10.** *The tracking of active power instruction in the droop fed DPC-driven VSPS system.*

#### **Figure 12.**

*Single-line diagram of the VSPS-wind farm-grid integrated system simulating setup.*

quickly and flexibly adjust the power fluctuations caused by the wind farm energy in the grid integrated system.

$$G(s) \approx \frac{\frac{P\_p}{I\_p}s + \mathbf{1}}{\frac{\left(\frac{1+P\_i}{I\_i}\right)k\_{drag}}{I\_p}s^2 + \frac{\left(P\_p+1\right)}{I\_p}s + \mathbf{1}}}\tag{33}$$

*Wind Power Fluctuation Compensation by Variable-Speed Pumped Storage Plant in a Grid… DOI: http://dx.doi.org/10.5772/intechopen.104938*

### **5. Simulation of models**

In this chapter, a case study is presented involving a 300 MW DFIM-based VSPS and a grid system integrated with a wind farm consisting of seven identical 15 MW wind turbine induction generators. To study the effect of wind fluctuations on grid frequency, a model of a medium-sized power system consisting of two 200MVA hydropower plants and one 15MVA diesel generator set and conventional synchronous generators was established. Phasor modeling techniques in MATLAB/Simulink are applied. For conventional synchronous machines, the primary voltage regulation method in automatic voltage regulators based on the standard IEEE type I and the primary frequency regulation method in turbine speed governors are used. The simulation setup is shown in **Figure 12**.
