Power flow simulation results with wind farm and SVC.

## 1. Introduction

Even though the use of wind generators for converting wind energy into electricity is beneficial from the environmental standpoint, their consideration in the active power dispatch makes the already complex task of achieving power system controllability even more demanding. Consequently, the quantification of the effects that large-scale integration of wind generation will cause on the network is a very important matter that requires special attention when planning and operating an electrical power system. Arguably, power flow analysis is the most popular computational calculation performed in a power system’s planning and operation, and this study has been selected to quantify the electrical response of wind generators when they are integrated in Flexible AC Transmission Systems.

Mathematical models of several types of wind generators have been developed in which their active and reactive power outputs are obtained based on the steady-state equivalent representation of the induction machine. The power injection method is then used to include these models into the power flow formulation, which is solved by using a sequential approach to obtain an operating point of the power system. In this approach, only the network’s state variables are calculated through a conventional power flow algorithm, while a subproblem is formulated for updating the state variables of wind generators as well as their power injections at the end of each power flow’s iteration.

Instead of using the power injection concept, another way of representing a wind generator is by means of an equivalent variable impedance expressed in terms of the slip of the generator and its rotor and stator winding parameters. This impedance is included in the system’s admittance matrix, and the network nodal voltages are computed through the power flow analysis. Based on these voltages, the air-gap power of the wind generator is calculated and used to iteratively compute the value of the slip of the induction generator that produces the match between the air-gap power and the mechanical power extracted from the wind.

In general terms, all the methods discussed above share the characteristic of using a sequential approach to calculate the state variables of the wind generators, and none of them considers the integration of FACTS controllers in the network where the WECSs are embedded.

A fundamentally different approach for the modelling of WECS, within the context of the power flow problem, is a method that simultaneously combines the state variables associated with the wind generators, the FACTS controllers and the transmission network in a single frame-of-reference for a unified iterative solution through a Newton-Raphson (NR) technique. From the convergence standpoint, the unified method is superior to the sequential one because the interaction between the network, FACTS controllers and wind generators is better represented during the iterative solution. Furthermore, it arrives at the solution with a quadratic convergence regardless of the network size. Hence, the key contribution of this work is to provide a comprehensive and general approach for the analysis of power flows in Flexible AC Transmission Systems containing wind generators in a unified single-frame of reference.

## 2. Power flow including FACTS controllers and wind generators

The unified approach suggested in [1] is extended to compute the power flow solution of a power system containing FACTS controllers and WECS furthermore, the approach is represented by a single set of nonlinear power flow mismatch equations *f(X* _{AC }, *X* _{F }, *X* _{WG }*)*=0, where *X* _{AC }is a vector of all nodal voltage magnitudes and angles, *X* _{F }stands for the state variables of the FACTS controllers and *X* _{WG }is a vector of all state variables associated with the wind generators. The linearised power flow mismatch equations corresponding to the wind farms are then combined with those associated with the FACTS controllers and the rest of the network, as given by (1), which are solved iteratively by the NR method:

where Δ*R* _{F }and Δ*R* _{WG }represent the mismatch equations of the FACTS controllers and wind generators, respectively. The NR method starts from an initial guess for all the state variables and updates the solution at each iteration *i* until a predefined tolerance is fulfilled. In this unified solution, all the state variables are adjusted simultaneously in order to compute the steady-state operating condition of the power system. Hence, this method has strong convergence characteristics.

## 3. Modelling of FACTS devices

Among all FACTS devices used to improve the steady-state performance of power systems [1], the Static Var Compensator (SVC) and Thyristor-Controlled Series Capacitor (TCSC) are the controllers considered.

### 3.1. Static VAR compensator

An integrated SVC and step-down transformer model is obtained by combining the admittances of both components *Y* _{T-SVC }*=* *Y* _{T }+ *Y* _{SVC }as proposed in [1]. The linearised power flow equations are given by (2) considering the firing angle *α* _{SVC }of the SVC as the state variable within the NR method

where *G* _{T-SVC} and *B* _{T-SVC }are functions dependent on *α* _{SVC }:

### 3.2. Thyristor-controlled series compensator

The TCSC firing angle power flow model is represented as an equivalent series reactance which is associated with the firing angle *α* _{TCSC. }This reactance can be expressed as [1] (see Appendix):

Assuming the TCSC controls the active power flowing from bus *k* to bus *m* to a specified value of

where the power flow mismatch for the TCSC module is defined as*B* _{km }=1/*X* _{TCSC }.

## 4. Modelling of wind generators

Wind generators are categorized according to how they operate when they are connected to the grid. The Fixed-Speed Wind Generators (FSWG) are thus named so because their speed is mainly set according to the system’s frequency [2]. In this category are the Stall-Regulated Fixed-Speed Wind Generators (SR-FSWG) and the Pitch-Regulated Fixed-Speed Wind Generators (PR-FSWG). A variant of the previous models is the semi-variable speed wind generator (SSWG), which uses a wound-rotor induction generator with an external resistor added to the rotor circuit in order to achieve a power regulation when wind speeds are above the rated one [3]. Also, variable-speed wind generators are being employed worldwide with the doubly-fed induction generator being the most used. However, another emergent topology that is being widely accepted is the wind generator based on a Permanent Magnet Synchronous Generator (PMSG) with a full-scale converter in which the gearbox can be omitted [4].

Mathematical modelling of FSWG, SSWG and PMSG-based wind generators for power flow studies is addressed below. In reference [5], the models of fixed- and semi-variable speed wind generators are suitably derived for power flow analysis and can be readily integrated in the formulation presented herein. For this reason, only a brief description of these models is given next.

### 4.1. Fixed-speed wind generators

This generator is directly connected to the network through a step-up transformer, and its final operating point depends upon the electrical frequency as well as the nodal voltage at the generator’s terminals. The generated reactive and active powers are determined by Equations (5) and (6), respectively, and the stator and rotor currents of the induction generator can be expressed according to Equations (7) and (8) [5] :

where *s* is the machine’s slip, *V* is the terminal voltage, and the constants from *A* to *W* are as given in the Appendix. Furthermore, the power converted from mechanical to electrical form (*P* _{conv }) can be computed by using Equation (9), where *R* _{2 }represents the rotor resistance

#### 4.1.1. Stall-regulated fixed-speed wind generator

The mechanical power *P* _{m }[W] extracted from the wind by this generator is given by [6]

where

*ρ*is the air density [kg/m^{3}],*A*is the swept area of the blades [m^{2}],*V*_{w }is the wind speed [m/s],*R*is the radius of the rotor [m],*n*_{gb}is the gearbox ratio,*ω*_{s }is the angular synchronous speed [rad/s],*β*is the pitch angle [degrees],*ω*_{T }is the angular speed of the turbine [rad/s],and the constants

*c*_{1 }to*c*_{9 }are the parameters of the wind turbine’s design.

Thus, assuming that the SR-FSWG is connected at bus *k*, the power mismatches equations are (11)-(13), and the set of linearised equations that has to be assembled and combined with the Jacobian matrix and the power mismatch vector of the entire network is shown in Equation (14) [5]:

where *P* _{g }(*V,s*) and *Q* _{g }(*V,s*) are given by (6) and (5), respectively, *P* _{LK }and *Q* _{LK }represent the active and reactive powers drawn by the load at bus *k*, respectively, and *P* _{k }^{cal }and *Q* _{k }^{cal }are active and reactive power injections given by

#### 4.1.2. Pitch-regulated fixed-speed wind generator

Since this wind generator has a blade pitch angle mechanism which actuates to limit the power extracted from the wind [7], the generated active power *P* _{g,pr }can be obtained from its power curve and is considered constant at a value *Q* _{g,pr }=*Q* _{g }needs to be calculated [5]. Therefore, the internal power equilibrium point in the wind generator has to be computed by Equation (17):

where *P* _{losses,s }and *P* _{losses,r }are the three-phase stator and rotor power losses, respectively, and the core losses in the induction machine are neglected. As mentioned previously, the set of linearised power flow mismatch equations regarding the PR-FSWG is given by Equations (18)-(21) when the generator is connected at bus *k*:

### 4.2. Semi-variable speed wind generator

In this type of generator, the slip of the induction machine cannot be regarded as the state variable because of the external resistance *R* _{ext }added in the rotor circuit [5]. Hence, a total rotor circuit resistance, *R* _{x }= (*R* _{2 }+*R* _{ext }) */s*, is considered as a single-state variable associated with the rotor circuit, which is adjusted to satisfy the power mismatch equations during the NR power flow calculation [5,8]. Hence, the reactive and active powers, the stator and rotor currents as well as the power converted from mechanical to electrical form will be dependent functions on *R* _{x }and can be expressed as

where the constants *A*´, *C*´, *E*´, *H*´, *K*´, *M*´ and *T*´ are given in the Appendix. Since the generated active power *P* _{m,ss }can be estimated as follows:

Therefore, when the SSWG is connected at bus *k*, its set of mismatch power flow equations is

where the internal energy balance in the induction machine Δ*P* _{WT3 }= -(*P* _{m,ss }- *P* _{conv }) is derived by using Equations (26)-(27).

### 4.3. PMSG-based wind generator

This type of wind generator possesses a PMSG and a full-rated converter to connect the generator to the network, resulting in complete speed and reactive power control [9]. Hence, all the generated power is supplied to the power system through a machine-side converter and grid-side converter. The schematic diagram of this topology is shown in figure 1(a). Reactive power support is one of the characteristics that make this machine attractive for wind power production. In this case, the inclusion of the explicit representation of the wind generator step-up transformer is considered, which allows for direct voltage magnitude control at the high-voltage side of the transformer. The proposed model for power flow studies is shown in the figure 1(b) in which *P* _{g,pmsg }represents the output power set by the wind generator power curve for a given wind speed, *V* _{msc }and *V* _{gsc }are the voltage at the machine-side converter and grid-side converter terminal, respectively, and *Z* _{st }is the step-up transformer impedance.

The power flow equations for the PMSG-based wind generator are derived assuming the following voltage at the grid-side converter terminal: V_{ gsc }= *V* _{gsc }*e* ^{jδgsc }. Based on Fig. 1(b), the active and reactive powers flowing from the grid-side converter terminal to the *k-th* bus are

For the active and reactive powers at bus *k*, the subscripts *gsc* and *k* are exchanged in Equations (32) and (33). Therefore, the NR-based power flow formulation is given by

Note that Equation (36) represents the power constraint in the AC/DC/AC converter in which active power losses are neglected.

## 5. Case studies with wind farms and FACTS devices

This section shows how the proposed approach performs when considering a power system having FACTS devices and wind generators.

### 5.1. Five-bus test system with FSWGs, SSWGs and a SVC

The typical five-bus test system is used to provide an example with the inclusion of a wind farm consisting of ten SR-FSWGs, ten PR-FSWGs and ten SSWGs operating at a wind speed of 16 m/s with which all wind generators are injecting their maximum power. Additionally, a SVC is placed at bus five in order maintain its terminal voltage magnitude at 1 pu. The conventional generators are set to control voltage magnitudes at 1 pu. Parameters of the wind farm are given in the Appendix.

In order to show the effect of wind farms and a SVC in the operation of the power system, the following three scenarios are presented: (a) the base case where the wind farm and SVC are not considered, (b) the case where only the wind farm is running and (c) the case where the SVC is connected at the bus at which the wind farm is connected in order to provide voltage support. The results regarding each case are reported in Table I.

Results | Scenario | ||

(a) | (b) | (c) | |

V1 | 1.000 | 1.000 | 1.000 |

V2 | 1.000 | 1.000 | 1.000 |

V3 | 0.971 | 0.971 | 0.977 |

V4 | 0.971 | 0.971 | 0.979 |

V5 | 0.967 | 0.966 | 1.000 |

V6 | --- | 0.945 | 0.981 |

Pg,WF | --- | 25.349 | 25.341 |

Qg,WF | --- | -10.802 | -9.857 |

Qg,SVC | --- | --- | 36.897 |

The simulated wind farm is a reactive power consumer since it lacks a reactive power control as seen from Table 1. When no SVC is considered, its reactive power absorption exceeds 10 MVAr, resulting in a low-voltage magnitude at node six. On the other hand, when the SVC is set in operation, not only the low-voltage side of the wind farm transformer is boosted, but also the system voltage profile.

Since all wind generators are operating at the same wind speed, clearly the state variables calculated for each group of wind generators corresponding to each model will result in the same value. Furthermore, if the system voltage profile changes, as occurs with the inclusion of the SVC, another operating point is found at each wind generator as shown in Table 2.

Scenario | SR-FSWG | PR-FSWG | SSWG | SVC |

s | s | Rx | αSVC (deg) | |

(b) | -0.00506 | -0.00679 | -67.69437 | --- |

(c) | -0.00467 | -0.00613 | -75.16952 | 136.35024 |

### 5.2. Five-bus test system with PMSG-based wind generators and a TCSC

In this case, a PMSG-based wind farm is located at bus 4 with 30 wind generators operating at a rated wind speed, i.e. 15 m/s. Be aware that each PMSG-based wind generator provides reactive power support by controlling its terminal voltage magnitude. Also, a TCSC is placed for controlling the active power flowing through the transmission line connected between nodes 4 and 5 at *P* _{4-5 }= 20 MW, as shown in figure 2. The next scenarios are analyzed: (a) the base case where the wind farm and TCSC are not considered, (b) the power system including only the PMSG-based wind farm and (c) the network having simultaneously the TCSC and the wind farm. The results are reported in Table 3.

The system voltage profile is improved when the PMSG-based wind farm is integrated to the grid as seen from Table 3. This is mainly due to two reasons; there is a redistribution of power flows in the network, e.g. the load connected at node four is being supplied locally by the wind farm. On the other hand, the PMSG-based wind farm is providing voltage support, resulting in a voltage magnitude of 1 pu at the low-voltage side of the wind farm transformer even when the TCSC is set in operation to increase the power transfer to 20 MW in the line connecting nodes 4 and 5. The wind generators and TCSC state variables estimated by the NR algorithm for each scenario are the following: scenario (b) *V* _{gsc }= 1.017 pu, *δ* _{gsc }= 9.071^{o}, and scenario (c) *V* _{gsc }= 1.016 pu, *δ* _{gsc }= 8.862^{o}, *α* _{TCSC }= 143.499^{o}.

Results | Scenario | ||

(a) | (b) | (c) | |

V1 | 1.000 | 1.000 | 1.000 |

V2 | 1.000 | 1.000 | 1.000 |

V3 | 0.971 | 0.984 | 0.984 |

V4 | 0.971 | 0.987 | 0.988 |

V5 | 0.967 | 0.973 | 0.971 |

V6 | --- | 1.000 | 1.000 |

Pg,WF | --- | 60.000 | 60.000 |

Qg,WF | --- | 2.566 | 2.218 |

P1-2 | 89.683 | 46.792 | 45.883 |

P1-3 | 40.639 | 20.252 | 21.206 |

P2-3 | 24.821 | 11.565 | 13.114 |

P2-4 | 28.076 | 9.279 | 11.253 |

P2-5 | 55.032 | 45.465 | 41.053 |

P4-3 | -18.686 | 13.628 | 11.171 |

P4-5 | 6.258 | 15.590 | 20.000 |

## 6. Discussion

This chapter has put forward the NR-based power flow algorithm which is capable of computing the steady-state operating point of electric networks containing WECS and FACTS devices. The solution problem is formulated in a single-frame of reference, resulting in an efficient iterative solution. Additionally, a PMSG-based wind generator model for power flow studies is presented, which allows for direct voltage magnitude control at the high-voltage side of the wind generator transformer. Numerical examples have shown that FACTS controllers are a practical alternative to integrate WECS into power systems without degrading their operative performance.

## Appendix

- TCSC parameters:

- FSWG parameters:

- SSWG parameters:

The data for each wind farm is (on a base power of 100 MVA): wind farm step-up transformer impedance is 0.2 pu, and the impedance of each wind generator transformer is 4.1667 pu. Also, the data of each wind generator model are given in Table 4

SR-FSWG | PR-FSWG | SSWG | |

Z1 | 0.0027 + j0.025 | 0.0 + j0.09985 | 0.00269 + j0.072605 |

Z2 | 0.0022 + j0.046 | 0.00373 + j0.10906 | 0.002199 + j0.04599 |

Zm | j1.38 | j3.54708 | j1.37997 |

Vnom | 690 | 690 | 690 |

Pnom | 900 | 600 | 1000 |

where *Z* _{1 }is the stator impedance [Ω], *Z* _{2 }is the rotor impedance [Ω], *Z* _{m }is the magnetizing impedance [Ω], *V* _{nom }is the rated voltage [V] and *P* _{nom }is the rated power of the wind generator [kW].

The coefficients of Equation (10) are as follows: *c* _{1 }= 0.5; *c* _{2 }= 116; *c* _{3 }= 0.4; *c* _{4 }= 0.0; *c* _{5 }= 0; *c* _{6 }= 5; *c* _{7 }=21; *c* _{8 }= 0.08; and *c* _{9 }= 0.035; *β =* 0. Also, for the PMSG-based wind generator, its rated voltage is 690 V, and its rated power is 2000 kW.

## Acknowledgements

The authors gratefully acknowledge the financial support granted to MSc. Luis M. Castro by the Consejo Nacional de Ciencia y Tecnología (CONACYT) México, and the University of Michoacán (U.M.S.N.H) for allowing him to undertake PhD studies. The authors gratefully acknowledge the financial support granted to Fuerte-Esquivel CR and Angeles-Camacho C by the FI and the II at the UNAM under the research project 2102.