Research Article  Open Access
Equivalent Modeling of DFIGBased Wind Power Plant Considering Crowbar Protection
Abstract
Crowbar conduction has an impact on the transient characteristics of a doubly fed induction generator (DFIG) in the shortcircuit fault condition. But crowbar protection is seldom considered in the aggregation method for equivalent modeling of DFIGbased wind power plants (WPPs). In this paper, the relationship between the growth of postfault rotor current and the amplitude of the terminal voltage dip is studied by analyzing the rotor current characteristics of a DFIG during the fault process. Then, a terminal voltage dip criterion which can identify crowbar conduction is proposed. Considering the different grid connection structures for single DFIG and WPP, the criterion is revised and the crowbar conduction is judged depending on the revised criterion. Furthermore, an aggregation model of the WPP is established based on the division principle of crowbar conduction. Finally, the proposed equivalent WPP is simulated on a DIgSILENT PowerFactory platform and the results are compared with those of the traditional equivalent WPPs and the detailed WPP. The simulation results show the effectiveness of the method for equivalent modeling of DFIGbased WPP when crowbar protection is also taken into account.
1. Introduction
With the rapid development of wind power generation and power electronics, the doubly fed induction generator (DFIG) has recently become the most commonly used wind turbine in wind power plant (WPP) based on its characteristics of maximum power point tracking (MPPT) [1], decoupled control of active and reactive powers [2], the use of a power convertor with a rated power of 25% of total system power [3], and so forth. To investigate the effect of the ongoing changes in power systems due to the increasing penetration of wind power, a wide range of studies, from steadystate to transient, are necessary. A detailed WPP model, in which the dynamics of each wind turbine and the internal network are fully represented, is not suitable as it could significantly increase the order of the mathematical model of the power system to be solved and thus influence the overall simulation time. As a result, the equivalent WPP model is generally recommended for reflecting the collective response of the whole WPP on large power systems.
Equivalent WPPs using aggregated wind turbines are classified as either the singlemachine representation or the multiplemachine representation [4]. In the first case, the actual WPP is modeled as a unique rescale wind turbine in the singlemachine representation. Accordingly, this is rational when all wind turbines are operating in an identical condition, whereas that might not be the case in realworld system operation. Nowadays, in order to improve the conformity of the singlemachine representation with the actual WPP, numerical identification methods for optimizing the parameters of equivalent models are applied. The optimized objects include the generator [5, 6], power convertor [7], and passive frequencydependent network [8]. The fuzzy logic system is also utilized to compute the mechanical torque compensating factor, which is integrated into a singlemachine aggregated WPP [9].
In the second case, study of the multiplemachine representation focuses on selecting reasonable clustering indices by quantifying and abstracting key features of wind turbines. The coherency method is presented in [10, 11] by clustering wind turbines with similar wind speed. In [12], the wind direction is implemented. To represent the power loss on the collective network within the WPP, a method of calculating the equivalent impedance is introduced in [13]. Moreover, research work in [14] analyzes the impact of the line impedance on the diversity in the voltage profile of each wind turbine and describes a voltageprofilebased approach to develop multiplemachine aggregated WPP. In addition, a clustering index considering the influence of fault types on coherency is depicted in [15].
Compared with conventional power resources, wind resources are stochastic and fluctuating. Equivalent models of an actual WPP might change due to different operating conditions. Therefore, the actual application should be based on historical data for the specific WPP and establish a probabilistic equivalent WPP [13, 16, 17]. This makes the equivalent WPP very practical and easy to use.
On the basis of previous equivalent methods, the research hotspot has gradually shifted to the validation and analysis of equivalent WPPs using field measurement data. The primary focus of reference [17] is to explain that the multiplemachine representation should take the diversity in protection conduction into account by investigating the discrepancies between the simulation results and the actual measurement. Perdana et al. [18] and Singh et al. [19] explore methods for validating the equivalent WPPs for SCIGbased and DFIGbased WPP, respectively. Meanwhile, reference [20] supports the fact that traditional clustering indices fail to cluster wind turbines when crowbar protection is enabled. Therefore, it is necessary to develop an equivalent method for modeling WPPs considering system protection.
This paper proposes a new equivalent method for modeling WPPs when the crowbar protection is taken into account. Based on the rotor current of the DFIG after an external threephase shortcircuit fault, the terminal voltage dip threshold is determined and from this it can be identified whether the crowbar protection is activated. Then, the threshold is revised based on the influence of the chain topology of the WPP. The revised terminal voltage dip threshold can be utilized to judge whether the crowbar is conducting. Furthermore, the equivalent WPP is established based on the division principle of crowbar conduction. Finally, the proposed equivalent WPP is simulated by the DIgSILENT PowerFactory platform and the results are compared with those of traditional equivalent WPPs and the detailed WPP. The simulation results validate the effectiveness of the equivalent method for modeling DFIGbased WPP.
2. Postfault Rotor Current Transient Analysis of DFIG
2.1. Postfault Rotor Current Calculation
Using the motor conversion, the stator and rotor voltages and fluxes in a synchronously rotating reference frame are given bywhere , , and represent voltage, current, and flux, respectively. is the resistance, and is the inductance. Subscripts and denote the stator and rotor quantities, respectively. is the mutual inductance. and are the synchronous speed and the slip speed, respectively.
One DFIG is connected to the infinite grid through a grounding transformer and a tertiary transformer, as shown in Figure 1. It is assumed that, at , an external threephase shortcircuit fault occurs, and, as a result, the terminal voltage of the DFIG abruptly changes from the steadystate voltage to the remaining voltage . The solution of (1) for can be written aswhere is the prefault rotor voltage. is the stator decay time constant. is the rotor decay time constant. and represent the rotor electrical angular speed and the integral constant, respectively, and .
When the DFIG operates in the stator voltage vector oriented control mode, , . Thus, the relationship between the output powers, voltage, and current can be written as where is the amplitude of the steadystate terminal voltage.
Based on (1) and (3), the prefault rotor voltage can be described as where
From (1), the rotor current is derived as follows:According to (2), (4), and (5), the dynamic representation for rotor current of the DFIG can be written as
2.2. Timer Action Crowbar
The crowbar circuit used in this paper is formed of a threephase rectifier, power resistor, and series IGBT switch [21]. The crowbar protection signal is triggered when the magnitude of the rotor current exceeds a threshold value set for crowbar protection insertion (2 pu) and then remains engaged for a fixed time (50 ms). When the crowbar is released, rotorside pulse width modulation (PWM) and rotor current PI control are immediately resumed. As shown in Figure 2, the DFIG wind turbine meets fault ride through requirements which refer to the capability of generation plant to remain connected and dynamically stable and offer network support throughout a serious voltage disturbance [22].
(a)
(b)
(c)
2.3. Terminal Voltage Dip Threshold
According to (6), when the DFIG is connected to the infinite grid, the postfault rotor current is jointly determined by the generator parameters, active and reactive power outputs, terminal voltage dip, steadystate terminal voltage, and generator speed. For the DFIG with the MPPT property, there is a corresponding relationship between the generator speed and the active power output. Thus, the influence of the generator speed on the rotor current can be included in the influence of the active power output. Once the generator parameters and the maximum rotor current for crowbar protection insertion are fixed, the terminal voltage dip is the function of only the power outputs and the steadystate terminal voltage. As a result, the terminal voltage dip threshold can be quantified when the crowbar protection is just activated. Crowbar conduction occurs if the terminal voltage dip is lower than the terminal voltage drop threshold. Otherwise, the crowbar protection is not activated.
We select 0.1 MW as the power step length within the scope of the DFIG rated power and 0.01 pu as the voltage step length within the amplitude scope of steadystate terminal voltage 0.97–1.07 pu [22]. The terminal voltage dip changes through regulating the external shortcircuit impedance and then finding out the terminal voltage dip threshold which makes the rotor current equal to the maximum rotor current for crowbar protection insertion. The terminal voltage dip thresholds of the DFIG whose power factor is −0.95, 1, and 0.95, respectively, are presented in Figure 3. The DFIG parameters are specified in Table 1.

(a)
(b)
(c)
3. Dynamic Equivalence of the WPP
3.1. Effect of Collector Line Impedance on Terminal Voltage Dip Threshold
In order to reduce the power loss of a lowvoltage circuit, it might be a better practice to connect multiple wind turbines to a middlevoltage bus using the daisychain structure within an actual WPP and then power infeed to the external grid from the WPP through transmission lines. In the same feeder, owing to the convergence effect of output current from other wind turbines, the voltage delta between a wind turbine and the middlevoltage bus has increased significantly, which is equivalent to increasing the electrical distance of the wind turbine to the middlevoltage bus.
Within the WPP, the differences in terminal voltages among the wind turbines are mainly caused by the different electrical distances from each wind turbine to the middlevoltage bus. Compared with the infinite grid connection structure for a single DFIG, in the WPP topology, the increasing effect of electrical distance between the wind turbine and the middlevoltage bus means the line impedance must not be ignored. At this time, the line impedance is jointed between the grounding transformer and the middlevoltage bus. Based on (6), if the impedance is incorporated into the generator stator side, the change in generator impedance parameters will have an impact on the rotor time constant, which will affect the relationship between the terminal voltage dip and rotor current. Therefore, the terminal voltage dip thresholds should be revised.
3.2. Terminal Voltage Dip Threshold Revision
In the same feeder, owing to the convergence of output currents among the wind turbines, the voltage drop has increased. In this paper, the virtual line impedance is introduced to characterize the equivalent electrical distance in the daisychain topological structure of the WPP.
As shown in Figure 4, there are n DFIGs in one feeder, and the voltage drop that is between the wind turbine number and the middlevoltage bus iswhere refers to the line impedance number ; is the conjugate value of apparent power of the wind turbine number ; is the voltage amplitude of the middlevoltage bus.
Based on the principle of equal voltage drop, the virtual line impedance which is between the wind turbine number and the middlevoltage bus isIn the infinite grid connection structure for a single DFIG, is jointed into the grounding transformer and the middlevoltage bus; then, the terminal voltage dip threshold is redetermined by means of the method in Section 2.3.
3.3. Clustering Principle of DFIGs
Based on the coherencybased method in the power system dynamic equivalence, crowbar conduction is selected as the clustering principle; for example, DFIGs are divided into two groups: (1) crowbar conduction during the fault; (2) crowbar nonconduction during the fault. Depending on the abovementioned analysis, whether the crowbar protection is activated or not can be judged by comparing the postfault terminal voltage dip of the DFIG with its terminal voltage dip threshold.
4. Simulation Results
The developed equivalent WPP is tested by comparing its transient response with the traditional equivalent WPPs and the detailed WPP.
4.1. WPP Description
The network of the WPP is presented in Figure 5. The WPP is composed of one 35/220 kV transformer, six collector systems, and 0.69 kV DFIGs that are connected to the collector system through 0.69/35 kV grounding transformers. Each collector system has 11 DFIGs, and the interval between two DFIGs within the same collector system is 400 m. Model parameters are given in Table 1.
4.2. Validation of Clustering Index
The measured wind data of the 66 DFIGs are listed in Table 2. All DFIGs operate with a unity power factor. Threephase shortcircuit faults are applied at point A, as shown in Figure 4. Faults start at 0.1 s and last for 150 ms. Two case studies are performed for different shortcircuit grounding impedance: (1) ; (2) ; therefore, the voltage of the 35 kV bus decreases to 53.7% and 55.7% of its nominal value. The postfault terminal voltage of the DFIG is compared with its terminal voltage dip threshold before revision and after revision, respectively, thus judging the DFIG whose crowbar protection is activated. Then, based on the simulation results of the detailed WPP, the correctness of the abovementioned judgment is tested, and the statistical results are listed in Table 3.


As can be seen from Table 3, when and , based on the terminal voltage dip threshold before revision, the numbers of wrong judgments of crowbar conduction are eight and seven, respectively; based on the terminal voltage dip threshold after revision, the numbers of wrong judgments of crowbar conduction are reduced to two and zero, respectively. The accuracy of the judgment of crowbar protection conduction can be significantly improved by using the terminal voltage dip threshold after revision.
4.3. Transient Response
The detailed WPP and the equivalent WPPs are modeled on the DIgSILENT PowerFactory platform. The detailed WPP consists of the 66 DFIGs, collection systems, grounding transformers, and a main transformer. Threephase shortcircuit faults are applied at point A, as shown in Figure 5. Faults start at 0.1 s and last for 150 ms.
Compared with the simulation results of the detailed WPP, error indicators are defined as follows:where and are the output variables of the complete collector system for the equivalent WPPs and the detailed WPP, respectively, and stands for the number of DFIGs within the WPP.
4.3.1. The Same Power Factor within WPP
All DFIGs operate with the same power factor, and two case studies are performed for different incoming winds: (1) 10.5 m/s, 0.5°; (2) 8.5 m/s, 174.5°. Based on wake effects, the incoming winds of each DFIG are listed in Tables 2 and 4, respectively. Three different reactivepower levels of DFIG, for example, the maximum production (0.95), unity power factor, and maximum consumption (−0.95), are used for the cases. Figures 6 and 7 illustrate the behavior of the voltages, active and reactive powers of the complete collector system for the three equivalent WPPs, and the detailed WPP following the application of a threephase shortcircuit fault.

(a)
(b)
(c)
(a)
(b)
(c)
The equivalent accuracy of the equivalent WPPs compared to the detailed WPP about the crowbar protection conduction and error indices are listed in Table 5. During the fault, whether the crowbar protection is triggered or not decides the reactive power supporting ability of the DFIG to the power grid, thereby affecting the depth of the voltage dip and the active power output of the WPP. The equivalent WPP built by this paper is more similar to the detailed WPP, due to its higher equivalent accuracy for crowbar conduction.

The crowbar conditions of the detailed WPPs are listed in Table 6. It should be noted that the DFIG might lose its capability to control the generator after crowbar conduction; thus, a DFIG whose generator speed is lower than the synchronous speed operates in the motor state. From fault elimination to crowbar protection resection, with the grid voltage recovery, a DFIG operating in the motor state will consume more active power, causing the active power at the complete collector system for the WPP to appear with a larger negative value. The equivalent WPP built by this paper can also accurately reflect this situation.

4.3.2. Different Power Factors among WPP Feeders
The incoming wind of the WPP is 10.5 m/s, 0.5°. Based on wake effects, the incoming winds of each DFIG are listed in Table 2. The power factors of the DFIGs in feeders 12, 34, and 56 are −0.95, 1, and 0.95, respectively. Threephase shortcircuit faults are applied at point A, as shown in Figure 5. Faults start at 0.1 s and last for 150 ms. Figure 8 shows the behavior of the active and reactive powers at the complete collector system for the equivalent WPPs and the detailed WPP following the application of a threephase shortcircuit fault. The crowbar condition of the detailed WPP is listed in Table 6. The equivalent accuracy of the equivalent WPPs compared to the detailed WPP about the crowbar protection conduction and error indices are listed in Table 7.

When the DFIGs operate with different power factors within the WPP, the traditional equivalent method will divide the DFIGs which have the same reactivepower level into groups. But the simulation results of the detailed WPP indicate that the crowbar conduction of DFIGs with the same power factor is not completely consistent.
5. Conclusions
After the fault, the crowbar conduction of the DFIG is related to the depth of the terminal voltage dip. This paper simplifies the judgment criterion of crowbar conduction as the terminal voltage dip threshold, which is convenient to calculate and makes it easy to identify the crowbar conduction.
Considering the influence of the connection line impedance within the same collector system on the terminal voltage dip threshold, the virtual line impedance is introduced to characterize the equivalent electrical distance from the DFIG to the mediumvoltage bus. The equivalent accuracy of crowbar conduction discrimination is improved by revising the terminal voltage dip threshold.
The equivalent WPP established based on the division principle of crowbar conduction is verified, showing that it can more accurately reflect the external characteristics of the WPP.
Appendix
Competing Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is supported by the National Natural Science Foundation of China (no. 51207039), the Science and Technology Foundation of SGCC (no. NY17201200073), and the National Energy Application Technology Research and Engineering Demonstration Projects (no. NY201104061).
References
 Z. Yu, M. E. Elbuluk, and Y. Sozer, “Stability analysis of maximum power point tracking (MPPT) method in wind power systems,” IEEE Transactions on Industry Applications, vol. 49, no. 3, pp. 1129–1136, 2013. View at: Publisher Site  Google Scholar
 M. Tazil, V. Kumar, R. C. Bansal et al., “Threephase doubly fed induction generators: an overview,” IET Electric Power Applications, vol. 4, no. 2, pp. 75–89, 2010. View at: Publisher Site  Google Scholar
 J. López, P. Sanchis, X. Roboam, and L. Marroyo, “Dynamic behavior of the doubly fed induction generator during threephase voltage dips,” IEEE Transactions on Energy Conversion, vol. 22, no. 3, pp. 709–717, 2007. View at: Publisher Site  Google Scholar
 E. Muljadi and A. Ellis, Final Project Report WECC Wind Generator Development, California Institute for Energy and Environment, 2010.
 M. A. Elizondo, L. Shuai, Z. Ning, and N. Samaan, “Model reduction, validation, and calibration of wind power plants for dynamic studies,” in Proceedings of the IEEE Power and Energy Society General Meeting, pp. 1–8, San Diego, Calif, USA, July 2011. View at: Google Scholar
 I. Erlich, F. Shewarega, C. Feltes, F. Koch, and J. Fortmann, “Determination of dynamic wind farm equivalents using heuristic optimization,” in Proceedings of the IEEE Power and Energy Society General Meeting (PES '12), pp. 1–8, IEEE, San Diego, Calif, USA, July 2012. View at: Publisher Site  Google Scholar
 Y.Q. Jin, P. Ju, and X.P. Pan, “Analysis on controller aggregation method for equivalent modeling of DFIGbased wind farm,” Automation of Electric Power Systems, vol. 38, no. 3, pp. 19–24, 2014. View at: Publisher Site  Google Scholar
 D. N. Hussein, M. Matar, and R. Iravani, “A type4 wind power plant equivalent model for the analysis of electromagnetic transients in power systems,” IEEE Transactions on Power Systems, vol. 28, no. 3, pp. 3096–3104, 2013. View at: Publisher Site  Google Scholar
 M. A. Chowdhury, W. X. Shen, N. Hosseinzadeh, and H. R. Pota, “A novel aggregated DFIG wind farm model using mechanical torque compensating factor,” Energy Conversion and Management, vol. 67, pp. 265–274, 2013. View at: Publisher Site  Google Scholar
 L. M. Fernández, C. A. García, J. R. Saenz, and F. Jurado, “Equivalent models of wind farms by using aggregated wind turbines and equivalent winds,” Energy Conversion and Management, vol. 50, no. 3, pp. 691–704, 2009. View at: Publisher Site  Google Scholar
 I. Zubia, J. X. Ostolaza, A. Susperregui, and J. J. Ugartemendia, “Multimachine transient modelling of wind farms: an essential approach to the study of fault conditions in the distribution network,” Applied Energy, vol. 89, no. 1, pp. 421–429, 2012. View at: Publisher Site  Google Scholar
 M. Ali, I.S. Ilie, J. V. Milanovic, and G. Chicco, “Wind farm model aggregation using probabilistic clustering,” IEEE Transactions on Power Systems, vol. 28, no. 1, pp. 309–316, 2013. View at: Publisher Site  Google Scholar
 E. Muljadi, C. P. Butterfield, A. Ellis et al., “Equivalencing the collector system of a large wind power plant,” in Proceedings of the IEEE Power Engineering Society General Meeting (PES '06), pp. 1–9, Montreal, Canada, June 2006. View at: Google Scholar
 Y. Cheng, M. Sahni, J. Conto, S.H. Huang, and J. Schmall, “Voltageprofilebased approach for developing collection system aggregated models for wind generation resources for grid voltage ridethrough studies,” IET Renewable Power Generation, vol. 5, no. 5, pp. 332–346, 2011. View at: Publisher Site  Google Scholar
 Q. Zhu, P. Han, M. Ding, X. Zhang, and W. Shi, “Probabilistic equivalent model for wind farms based on clusteringdiscriminant analysis,” Proceedings of the Chinese Society of Electrical Engineering, vol. 34, no. 28, pp. 4770–4780, 2014. View at: Publisher Site  Google Scholar
 E. Muljadi, Z. Mills, R. Foster, J. Conto, and A. Ellis, “Fault analysis at a wind power plant for one year observation,” in Proceedings of the IEEE Power Energy Society General Meeting, pp. 1–7, Pittsburgh, Pa, USA, July 2008. View at: Google Scholar
 Y. Zhang, E. Muljadi, D. Kosterev, and M. Singh, “Wind power plant model validation using synchrophasor measurements at the point of interconnection,” IEEE Transactions on Sustainable Energy, vol. 6, no. 3, pp. 984–992, 2015. View at: Publisher Site  Google Scholar
 A. Perdana, S. UskiJoutsenvuo, O. Carlson, and B. Lemström, “Comparison of an aggregated model of a wind farm consisting of fixedspeed wind turbines with field measurement,” Wind Energy, vol. 11, no. 1, pp. 13–27, 2008. View at: Publisher Site  Google Scholar
 M. Singh, K. Faria, S. Santoso, and E. Muljadi, “Validation and analysis of wind power plant models using shortcircuit field measurement data,” in Proceedings of the IEEE Power & Energy Society General Meeting (PES '09), pp. 1–6, Calgary, Canada, July 2009. View at: Publisher Site  Google Scholar
 J. Brochu, C. Larose, and R. Gagnon, “Validation of single and multiplemachine equivalents for modeling wind power plants,” IEEE Transactions on Energy Conversion, vol. 26, no. 2, pp. 532–541, 2011. View at: Publisher Site  Google Scholar
 G. Pannell, D. J. Atkinson, and B. Zahawi, “Minimumthreshold crowbar for a faultridethrough gridcodecompliant DFIG wind turbine,” IEEE Transactions on Energy Conversion, vol. 25, no. 3, pp. 750–759, 2010. View at: Publisher Site  Google Scholar
 State Grid Corporation of China, Technology Rule for Connecting Wind Farm Into Power Grid, State Grid Corporation of China, Beijing, China, 2009.
Copyright
Copyright © 2016 Qianlong Zhu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.