A New Wind Turbine Generating System Model for Balanced and Unbalanced Distribution Systems Load Flow Analysis

: Wind turbine generating systems (WTGSs), which are conventionally connected to high voltage transmission networks, have frequently been employed as distributed generation units in today’s distribution networks. In practice, the distribution networks always have unbalanced bus voltages and line currents due to uneven distribution of single or double phase loads over three phases and asymmetry of the lines, etc. Accordingly, in this study, for the load ﬂow analysis of the distribution networks, Conventional Fixed speed Induction Generator (CFIG) based WTGS, one of the most widely used WTGS types, is modelled under unbalanced voltage conditions. The Developed model has active and reactive power expressions in terms of induction machine impedance parameters, terminal voltages and input power. The validity of the Developed model is conﬁrmed with the experimental results obtained in a test system. The results of the slip calculation based phase-domain model (SCP Model), which was previously proposed in the literature for CFIG based WTGSs under unbalanced voltages, are also given for the comparison. Finally, the Developed model and the SCP model are implemented in the load ﬂow analysis of the IEEE 34 bus test system with the CFIG based WTGSs and unbalanced loads. Thus, it is clearly pointed out that the results of the load ﬂow analysis implemented with both models are very close to each other, and the Developed model is computationally more efﬁcient than the SCP model.


Introduction
Environmental damage and declining reserves of fossil fuels, and the risks of the usage of nuclear energy, increases the interest shown in renewable energy sources day by day. Many countries have considered renewable energy sources to meet the energy needs in order to reduce the dependence on fossil fuels and foreign sources. These countries make medium/long-term plans to increase the share of renewable energy sources in energy production. Governments are applied incentives to steer the private investors on renewable energy sources by the scope of these plans, and, thus, the investments are made on domestic energy production. The European Union (EU) aims to achieve at least 20% of its energy needs from renewable energy sources until 2020, and this rate is intended to reach at least 27% for the year 2030 [1]. In the renewable electricity plan that has been published by the US Energy Efficiency and Renewable Energy Office, at least 30% of the energy needs of US are intended to be obtained from renewable energy sources by 2025 [2]. In the renewable energy action plan that was announced by the General Directorate of Renewable Energy of the Republic of Turkey, at least 30% of the energy needs are intended to be derived from renewable energy sources by 2023 [3].
Wind energy, an important renewable energy source, is a resource where countries can benefit from the advantages of their geographical location. Wind energy is converted into electrical energy by the wind turbine generator systems (WTGSs) [4]. WTGSs are generally formed with a propeller and a generator. When compared to synchronous generators, being smaller in size and having lower operation costs, asynchronous generators are widely preferred in small sized WTGSs. WTGSs can be classified under two main headings that are conventional fixed speed WTGSs and variable speed WTGSs [5]. Variable speed WTGSs have high installation and operating costs since they include high-tech power electronic components and control circuits. As a result of this, instead of variable speed induction generators, CFSIGs are preferred and are generally connected to low or medium voltage distribution networks directly.
Large-scale wind farms are typically connected to the high-voltage transmission system; however, small and medium-scale distributed generation units based on wind energy are placed in the distribution systems. Generally, distribution systems are structured and operated radially; thus, it is expected to be supplied from a single source. With the increasing usage of WTGSs, multiple feeding points are emerging on the distribution system. This can cause power flow in both directions (towards source or load), unlike the conventional system operation [6]. For example, in a case where the distribution network is dominated by wind energy resources, this distribution network can behave as an active network and can give feedback to the high voltage bus. As a distributed generation system, not only WTGSs, but also other renewable or non-renewable generation systems have an impact on the conventional distribution network's voltage profile and energy losses [6,7]. It should be acknowledged that these kinds of operations can cause a number of technical difficulties and operational issues [8].
In practice, the distribution networks always have unbalanced bus voltages and line currents due to single or double phase loads, unsymmetrical line impedances, single or two phase spurs, etc. [9]. There are balanced power flow models for directly network connected CFIGs in the literature. [10,11]. The use of these models for unbalanced distribution system power flow analysis may lead to erroneous results [12]. Therefore, a symmetrical component based model of induction generators is introduced in the study of Abdel-Akher et al. [12] and Ghorashi et al. [13]. The accurate estimation of the rotor slip value is difficult and requires computationally expensive iterative methods. In addition to this, the need for a solution for complex phase currents/voltages brings extra computational complexity for power flow analysis. Therefore, in most of the unbalanced power flow analyses in the literature, asynchronous generators are modeled as balanced and fixed Active (P)-Reactive (Q) Power sources by ignoring the effect of the imbalance (Constant P, Q Model) [14][15][16]. In some studies, this simplified model has been updated by the calculation of the reactive power demand from the grid. For this purpose, the reactive power demanded from the grid brought into a function of positive sequence voltage (V sp ) and machine impedance parameters, and, thereby, this updated model has been applied to power flow analysis (Constant P, Q(V sp ) Model) [14][15][16].
Preliminary results of this paper were presented in a conference [17]. In addition to preliminary results, this study includes load flow analysis and experimental validation of the model. In this paper, for the load flow analysis of the distribution networks, CFIG based WTGSs are modelled under unbalanced voltage conditions. The Developed model has active and reactive power expressions in terms of induction machine impedance parameters, input power and terminal voltages. In the Developed model, magnitude of the positive-sequence rotor side voltage is expressed via bi-quadratic equation, which is generally used to calculate node voltages in the load flow analysis [10,11,14,18,19]. In the study of Feijôo and Villanueva [18], a slightly modified model of the Eminoglu et al. study [10] is given and the advantages of the bi-quadratic equation based models are highlighted. In addition, such models have an iterative algorithm to find the positive-and negative-sequence rotor side powers, the sum of which is equal to the wind power. Unlike the slip calculation based models, the Developed model does not require prior knowledge on the angles of sequence voltages and currents, thus it can be implemented with a lesser number of operations.
To verify the Developed model for WTGSs, a laboratory testing was carried out. An induction generator coupled with a controllable induction motor was used. Finally, both the Developed Model and the SCP Model [12] are applied by the authors to phase domain power flow analysis in a 34 bus radial distribution test system published by IEEE [20], and the results were analyzed comparatively. In addition, the complexity analysis of these two models are compared in terms of required mathematical operations. The three-phase load flow analysis of models was performed via Matlab 2012a [21] by using open source software OpenDSS (Open Distribution System Simulator) [22] program libraries. The library files of models (DLLs (Dynamic Link Library)) were created in Delphi ® XE5.

Developed Model for WTGSs
In this study, it is aimed to model the CFSIG based WTGS for unbalanced load flow analysis of distribution power systems. It has a delta configured induction generator directly connected to the bus. The mechanical input power (P T ) of a WTGS can be calculated by regarding the air density (ρ), the area swept by the rotor (A), the power coefficient (C P ) and wind speed (u) [12]: The power coefficient is subject to the tip speed ratio (λ) expressed as follows: where ω and R denote the angular velocity and turbine rotor radius, respectively. To provide the model of WTGS, firstly, the positive-and negative-sequence equivalent circuits of the induction generator shown in Figure 1 are taken into consideration. In the same figure, V sp , V sn , V rp and V rn indicate the magnitudes of the positive-and negative-sequence voltages for the stator and referred rotor sides, respectively. For the equivalent circuits, R sc is the short-circuit equivalent resistance, which is sum of the stator winding (R s ) and referred rotor winding resistance (R r ), and X sc is the short circuit equivalent reactance, which is sum of the stator winding (X s ) and referred rotor winding reactance (X r ). R m and X m stand for the resistance and the magnetizing reactance of the core for the positive-sequence equivalent circuit. On the other hand, by looking at Ref. [13], in the negative-sequence equivalent circuit, R m and X m are neglected for the simplification. P sp and P sn stand for the positive-and negative-sequence active powers absorbed by the stator.
The sum of the positive-and negative-sequence rotor active powers, denoted as P rp and P rn , gives the mechanical input power of the WTGS: ( It is seen from the equivalent circuits that the induction generator draws the positive-and negative-sequence reactive powers (Q sp and Q sn ) from the grid. In addition, the positive-and negative-sequence reactive powers (Q rp and Q rn ), which are measured at the rotor (r) bus, should be nil. While looking at the positive-sequence equivalent circuit and the bi-quadratic equation, which was used to model the CFSIG based WTGSs for balanced systems in [9], the referred rotor positive-sequence voltage's magnitude, V rp , can be written as where Z sc stands for the induction generator's short-circuit impedance: The magnitude of the negative sequence current (I sn ), which flows from r bus to s bus in the negative-sequence equivalent circuit, and negative-sequence active power at r bus (P rn ) can be expressed in terms of the impedance parameters (R r , R sc , X sc ), the magnitudes of the rotor positive-sequence and stator negative-sequence voltages (V rp , V sn ) and positive-sequence active power at r bus (P rp ): The positive-and negative-sequence active powers at s bus can be expressed as in Equations (8) and (9): Hence, the total active power at s bus can be calculated as: On the other hand, the positive-and negative-sequence reactive powers at s bus can be expressed as follows: By regarding both reactive powers, the total reactive power at s bus can be calculated: Since P rp and P rn have unknown values, the model needs an iterative solution algorithm. The algorithm of the Developed model is summarized below:

•
Step 1: Determine initial value of P rp by assuming all of the mechanical power (P T ) delivered to the positive-sequence circuit.

•
Step 3: Find I rn and P rn by substituting the V rp and the P rp values in Equations (6) and (7). • Step 4: Calculate P s and Q s values by means of Equations (8)-(13).

•
Step 5: Find the relative difference value of P s and the relative difference value of Q s , which are calculated for the last two iterations (i + 1. and i. iterations), via Equations (14) and (15). • Step 6: If both calculated relative difference values are larger than tolerance value (ε), update the P rp value by substituting the last calculated P rn value in Equation (3), and return to the V rp calculation step. Otherwise, finalize the solution algorithm: Note that the magnitudes (V sp , V sn ) and angles (ϕ sp , ϕ sn ) of the sequence voltages at the s bus have known values, and sequence active powers (P sp , P sn ), sequence reactive powers (Q sp , Q sn ) and magnitudes (I sp , I sn ) of sequence currents at the s bus are found after completing the iterative algorithm. Additionally, the angles of the positive-and negative-sequence currents measured at the s bus, denoted as β sp and β sn , can be calculated by considering Equations (16) and (17): Finally, sequence voltages and currents of the stator side can be converted to the phase values of those quantities by means of the symmetrical components method. Then, using Equations (18) and (19), active and reactive powers can be found for m = a, b, c phases: In Equations (18) and (19), θ m denotes phase angle difference between voltage and current of phase m. If phase active and reactive powers have a minus (−) sign, it can be mentioned that these powers are injected from generator to grid; otherwise, these powers flow from grid to the generator. Note that, in the next sections, due to the fact that the injected active and demanded reactive power are given in the analysis, both of them will have a positive sign.

Load Flow Analysis Implementation of Developed Model
The details of the algorithm are provided in this section to incorporate the developed CFIG based WTGS model in load flow analyses. It is divided into two iterative blocks as given in Figures 2  and 3. The first part of the blocks is three-phase unbalanced load flow block, and the second part is provided for the Developed model. The unbalanced load flow analysis for the Developed model can be performed by a sweep based algorithm that is commonly used in distribution systems [23].  The first step at the implementation of the Developed model into the load flow analysis is the data preparation of related line and transformers' impedances, loads and generator data and initial values of the bus voltages, which are assumed to be 1 pu. In the second step, three-phase unbalanced load flow is performed. After that, generator terminal voltages, which are necessary for asynchronous machine iteration block (machine iteration block), are calculated. By converting to positive-and negative-sequence components from calculated generator terminal phase voltages, machine iteration starts and this is shown in Figure 3 after number (1). It will be stopped when Equations (14) and (15) are satisfied. Finally, the model algorithm provides the phase-domain quantities for the load flow algorithm at number (2).

Validation of Model
To show the validity of the Developed model, its results are comparatively evaluated with the experimental results obtained in a test system. In addition, the results of the SCP model [12] are also provided. The equivalent circuit parameters of the test machine are given in Table A1, and these parameters are obtained by standard no-load and blocked rotor tests. In the system, active and reactive powers of the asynchronous generator are measured under unbalanced voltages. For the same unbalanced voltages, the results of the Developed model and SCP model [12] are calculated, and they are comparatively evaluated by considering the measurement results as reference values.
In the comparative analysis, the voltage unbalance factor (VUF) [24,25] is considered to quantify the unbalance level of the test voltages. The expression of VUF can be written as: A schematic representation and a photograph of the test system, which are used for validity tests of the model, are taken and provided in Figure 4. In the first part of the validity tests, for the induction generator working with the constant rotor slip value (s = −1.2%) under the VUF values between 0% and 5%, the injected active and demanded reactive powers are plotted in Figure 5. Thus, the Developed model and SCP model [12] are examined in how they responded to VUF changes. In the second part of the tests, the asynchronous generator working with the slip values between −0.4% and −1.2% under the unbalanced terminal voltage with VUF = 5%, the injected active and demanded reactive powers are plotted in Figure 6. The VUF values have been obtained by changing the magnitudes of phase voltages while the phase angle differences remained constant at 120 degrees. In addition to that, the tolerance value (ε) of the machine iterations for validity tests has been determined as 10 −9 for both models.  The obtained and calculated reactive power demand results at VUF 0% given in Figure 5a are 0.720 pu, 0.665 pu and 0.655 pu, respectively, for Measurement, the Developed Model and the SCP Model [12]. In the same figure, reactive power demand results at VUF 5% are 0.653 pu, 0.582 pu and 0.608 pu, respectively, for Measurement, the Developed Model and the SCP Model [12]. Similarly, injected active power results at VUF 0% seen in Figure 5b are 0.405 pu, 0.366 pu and 0.358 pu, respectively, for Measurement, the Developed Model and the SCP Model [12]. In the same figure, injected active power results at VUF 5% are 0.364 pu, 0.347 pu and 0.320 pu, respectively, for Measurement, the Developed Model and the SCP Model [12]. In light of this information, as a result of neglecting Rm and Xm values which are not readily available and difficult to determine, for simplification in the negative sequence equivalent circuit, at higher VUF values, the reactive power demand results obtained by the Developed Model are moving away from the trend of the measured values. For higher VUF values, a similar situation is also observed in the results obtained by the SCP model [12].
The obtained and calculated reactive power demand results at −0.4% slip value given in Figure 6a are 0.596 pu, 0.557 pu and 0.571 pu, respectively, for Measurement, the Developed Model and the SCP Model [12]. In the same figure, reactive power demand results at −1.20% slip value are 0.653 pu, 0.582 pu and 0.608 pu, respectively, for Measurement, the Developed Model and the SCP Model [12].
Similarly, injected active power results at −0.40% slip value seen in Figure 6b are 0.087 pu, 0.103 pu and 0.080 pu, respectively, for Measurement, the Developed Model and the SCP Model [12]. In the same figure, injected active power results at −1.20% slip value are 0.364 pu, 0.347 pu and 0.320 pu, respectively, for Measurement, the Developed Model and the SCP Model [12]. Despite the change in rotor slip value, the injected active power and reactive power demand given in Figure 6a,b show that the Developed model and the SCP model give close results, and these results have the same trend as the measured ones. Actually, at BS EN 50160 [26], it is recommended that the value of the VUF is less than 2% in the distribution systems. For this reason, in order to evaluate the Developed model better, the test results of the asynchronous generator working with the slip values between −0.4% and −1.2% under the unbalanced terminal voltage with VUF = 2% are also given in Figure 7.
The obtained and calculated reactive power demand results at −0.4% slip value given in Figure 7a are 0.625 pu, 0.606 pu and 0.592 pu, respectively, for Measurement, the Developed Model and the SCP Model [12]. In the same figure, reactive power demand results at −1.20% slip value are 0.684 pu, 0.632 pu and 0.630 pu, respectively, for Measurement, the Developed Model and the SCP Model [12]. Similarly, injected active power results at −0.4% slip value seen in Figure 7b are 0.087 pu, 0.083 pu and 0.073 pu, respectively, for Measurement, the Developed Model and the SCP Model [12]. In the same figure, injected active power results at −1.20% slip value are 0.390 pu, 0.359 pu and 0.346 pu respectively for Measurement, the Developed Model and the SCP Model [12]. It is seen in Figure 7 that the results obtained by the Developed model are in close agreement with the SCP model [12] and show the same trends and similar trends as the measured ones.

Load Flow Analysis Results for the IEEE 34 Bus Test System
In this section, it is aimed to comparatively evaluate the results of the load flow analysis achieved with the Developed model and SCP model [12] for the IEEE 34 bus radial distribution test system, prepared by IEEE [20] and modified by Dugan et al. [27]. The modification of the test system was carried out by adding two generators with two service transformers to its two buses in Dugan et al. [27].
The single-line diagram of this test system is given in Figure 8. Here, it should be noted that the bus 800 refers to the main feeder in Figure 8. The wind turbines are placed at the buses 848 and 890 with the proper service transformers, and their parameters are given in Table A2. Wind turbines with induction generators have been planned as having 660 kW rated output power, and all of their parameters are included in Table A3. Figure 8. 34 bus radial test system with the placement of wind turbines (redrawed from [27]). For the test system detailed above, load flow analysis is performed using two models and thus phase voltages of each buses are obtained to compare with the phase voltages provided by Dugan et al. [27]. The tolerance value (ε) of the machine iterations at each sweep in power flow analyses has been determined as 10 −9 for both models. In addition, the tolerance value (ε) for the power flow analyses has been determined as 10 −4 .
The percentage difference of the phase voltages and angle differences for the each buses are given in Figures 9 and 10 respectively. In Figure 9, the overall differences of bus voltage do not exceed 0.23%, 0.43% and 0.28%, respectively, for phases A, B and C. Additionally, in Figure 10, the overall differences of bus voltage angle values do not exceed −0.6 • , 0.5 • and −0.5 • , respectively, for phases A, B and C.
The highest magnitude difference values of Developed Model and SCP Model [12] for phase A is 0.23% and 0.14%, for phase B is 0.43% and 0.09%, and for phase C is 0.28% and 0.08%. In addition, the highest angle difference values of the Developed Model and the SCP Model [12] for phase A is −0.6 • and −0.4 • , for phase B is 0.5 • and 0.4 • , and for phase C is −0.3 • and −0.5 • . It is seen from Figures 9 and 10 that the Developed model is in close agreement with the SCP Model [12].

Computational Performance of the Models
The iteration numbers and complexities of the studied models for the load flow of the 34 bus test system are investigated in this section. Accordingly, two terms related with iteration numbers such as general and maximum machine iterations are given in Table 1. General iteration number describes the number of iterations required in order to reach the solution of the load flow analysis and the machine iteration number represents the number of iterations required to meet the convergence criteria of the machine iteration within the load flow algorithm. The maximum machine iteration shows the maximum number of machine iterations that occurred during the load flow iterations. Table 1 shows that the load flow algorithms accompanied with the Developed Model and the SCP Model [10] converge in 11 and 12 iterations, respectively. In addition, the maximum number of machine iterations for the Developed Model and the SCP Model [12] are, respectively, two and six iterations. On the other side, numbers of operations, which are required by the models for the load flow analysis of the 34 bus test system, are shown in Figure 11. It is seen from this figure that the SCP Model [12] requires numbers of addition, subtraction, multiplication, division and square root operations as 3840, 2112, 12,672, 1272 and 744, respectively. In addition to this, the Developed model requires numbers of the respective operations as 1936, 1276, 3828, 880 and 176. Figure 11. Executed operation numbers for models in the 34 bus test system analysis.
According to the results given up to here, it can clearly be mentioned that the Developed Model is computationally more efficient than the SCP Model [12] for the load flow analysis.

Conclusions
In this study, a model to be used in balanced or unbalanced distribution system load flow analysis for CFIG based WTGSs, which is directly connected to the grid, is developed, and it is also incorporated into a sweep based load flow analysis algorithm. The accuracy of the Developed model is experimentally evaluated and extensive validation is performed against the SCP Model [12] previously presented in the literature.
Considering the experimental analysis results obtained under different voltage unbalance and rotor slip levels, it is found that the Developed Model successfully estimates the reactive power demand of the machine and provides a meaningful figure for the generated active power when compared with the measured values and the SCP Model [12]. Thus, it can be concluded that the Developed model is robust and facilitates the computation of real and reactive power outputs of the induction generator for a specified wind speed and the terminal voltage by recognizing that machine equivalent circuit parameters are available.
In addition, the Developed Model and the SCP Model [12] are used for the load flow analysis of the 34 bus unbalanced radial distribution test system. Accordingly, it is clearly pointed out that both models give results very close to the load flow results of the 34 bus test system provided by Dugan et al. [27]. Models are also compared in terms of operation number and complexity for the load flow analysis of the 34 bus test system. As a result, it is concluded that the Developed Model requires less iterations and number of operations when compared with the SCP Model [12]. Table A1. Asynchronous generator parameters used in the experimental system.

Type of Parameters Values
Power ( Table A2. Service transformer parameters used in the 34 bus test system.