LiMingguang
                     (Mingguang Li)
                     1
                     YangDejian
                     (Dejian Yang)
                     1
                     강용철
                     (Yong-Cheol Kang)
                     3†
                     홍준희
                     (Junhee Hong)
                     2
               
                  - 
                           
                        (Dept. of Electrical Engineering, Chonbuk National University, Korea)
                        
 
                  - 
                           
                        (Dept. of Energy IT, Gachon University, Korea)
                        
 
                  - 
                           
                        (Dept. of Energy IT, Gachon University, Korea)
                        
 
               
             
            
            
            Copyright © The Korean Institute of Electrical Engineers(KIEE)
            
            
            
            
            
               
                  
Key words
               
               Variable-speed wind turbine, Power smoothing, Frequency mitigation, Over-frequency section, Control gain
             
            
          
         
            
                  1. Introduction
               For an electric power system with a high penetration level of wind power, the significant
                  variation of wind speed results in difficulties in keeping the system frequency within
                  a narrow range(1). This is because variable-speed wind turbines (VSWTs) such as doubly-fed induction
                  generators (DFIGs) and fully-rated converter-based VSWTs perform maximum power point
                  tracking (MPPT) operation(2). To mitigate these problems, some countries including Korea have specified the requirements
                  on the ramp rates of the output power of a wind power plant (WPP)(3,4).
               
               To mitigate the frequency fluctuation, if the system frequency exceeds the nominal
                  frequency (over-frequency section, OFS), the output power generated from a VSWT should
                  be reduced; conversely, in the under-frequency section (UFS) the output power from
                  a VSWT should be increased. A number of methods have been proposed that can mitigate
                  the frequency fluctuation by smoothing the output power of a VSWT(5-11); these solutions can be roughly divided into two types(11): those with or without energy storage systems (ESSs). In(5-7), ESSs such as supercapacitors, flywheels, or batteries were suggested to mitigate
                  the output power fluctuation of a VSWT. These additional devices effectively mitigate
                  the power fluctuation, however, an extra cost for installation and maintenance is
                  required, particularly for a large-scale VSWT.
               
               To avoid this, the power-smoothing schemes were suggested that utilizes the rotating
                  masses of a VSWT as an ESS. These schemes employ an additional control loop operating
                  in conjunction with the MPPT control loop(8-10). They can mitigate the output power fluctuation by releasing the stored energy in
                  the rotating masses of a VSWT to the grid in the UFS or storing energy into the rotating
                  masses in the OFS. This means that the heavy rotating masses of a VSWT, which has
                  the comparable inertia constant to the conventional synchronous generators, can be
                  effectively utilized as an ESS. Additional control loops were suggested based on the
                  measured frequency: the rate of change of frequency loop and/or frequency deviation
                  ($\Delta f$) loop(8-10). In(8,9), a fixed gain is used for the additional control loops. The use of a large gain can
                  improve the power- smoothing capability of a VSWT, but it might result in over-deceleration
                  of the rotor speed in the low-rotor-speed region. In this case, the output power of
                  a VSWT should be abruptly reduced to recover the rotor speed; consequently, this causes
                  a significant frequency fluctuation in the power system. To avoid this, the use of
                  a small gain is inevitable, thereby it provides a limited contribution to mitigating
                  the frequency fluctuation. To improve the power-smoothing capability while avoiding
                  over-deceleration of the rotor speed, a power- smoothing scheme was proposed that
                  employs a variable gain of the rotor speed(10). In this scheme, the control gain is set to be proportional to the rotor speed; thus,
                  it can effectively release the stored energy in a VSWT to smooth power in the UFS.
                  However, this scheme uses the same control gain for both the OFS and UFS. Thus, even
                  though in the low-rotor-speed region a VSWT has a large potential for storing energy
                  into the rotating masses, this scheme provides a limited contribution to storing energy
                  in the OFS because of a small gain.
               
               This paper proposes a power-smoothing scheme of a VSWT that improves the frequency-mitigating
                  capability in the OFS. The proposed scheme employs an additional control loop based
                  on the frequency deviation operating in conjunction with the MPPT control loop. To
                  improve the frequency-smoothing capability, the control gain of the additional control
                  loop for the OFS is proposed that is set to be inversely proportional to the rotor
                  speed. In contrast, the control gain in the UFS is set to be proportional to the rotor
                  speed. The performance of the proposed scheme is investigated under various wind penetration
                  levels using an EMTP-RV simulator.
               
             
            
                  2. Power-Smoothing Schemes of a DFIG
               This section describes the overall features of the two conventional power-smoothing
                  schemes suggested in[9] and[10], and the proposed scheme, which use the $\Delta f$ loop for power smoothing. In this
                  paper, the scheme in[9] is denoted as Scheme #1 and the scheme in[10] as Scheme #2.
               
               The power reference ($P_{ref}$) in Scheme #1, Scheme #2, and the proposed scheme consists
                  of the output for the MPPT control loop ($P_{MPPT}$) and the output of the $\Delta
                  f$ loop ($\Delta P$) as in:
               
               
               where $\omega_r$ is the rotor speed; and $k_g$ is constant and is set to 0.512.
               This paper aims to improve the power-smoothing capability of a DFIG, particularly
                  in the OFS. Fig. 1 shows the operational features of the proposed scheme. While Scheme #2 uses the same
                  control gain for both the UFS and OFS, the proposed scheme employs the different control
                  gains separately defined in the UFS and OFS, as in:
               
               
               
               where $K_{UFS}$($\omega_r$) and $K_{OFS}$($\omega_r$) are the control gains in the
                  UFS and OFS, respectively; $\omega_{min}$ and $\omega_{max}$ are set to 0.7 p.u. and
                  1.25 p.u., respectively.
               
               Thus, $\Delta P$ in the UFS and OFS can be obtained as:
               
               $K_{UFS}$($\omega_r$), which represents the lower straight line between ($\omega_{min}$,
                  0) and ($\omega_{max}$, 300) in Fig. 1(b), is deliberately set to be proportional to $\omega_r$ while setting $K_{UFS}$($\omega_{min}$)
                  to be zero so that the power-smoothing capability of the proposed scheme is similar
                  to Scheme #2. The reason for this is as follows. If $K_{UFS}$($\omega_{max}$) is set
                  to be larger than 300, the proposed scheme can provide better power - smoothing capability
                  even in the UFS than Scheme #2 by releasing a larger amount of energy. However, this
                  paper aims to improve the power- smoothing capability by employing the different control
                  gain in the OFS while keeping the power-smoothing capability in the UFS similar to
                  that of Scheme #2. As shown in Fig. 1(b), $K_{UFS}$($\omega_r$) is slightly larger than that in Scheme #2.
               
               Conversely, in the OFS a DFIG should store energy into the rotating masses to mitigate
                  the output power fluctuation of a DFIG, thereby reducing the frequency deviation.
                  In the low-rotor-speed region, a DFIG has a larger potential for storing energy into
                  the rotating masses than in the high-rotor-speed region and thus the control gain
                  in the low-rotor-speed region should be set to be larger. This means that the control
                  gain for the OFS should be monotonously decreasing with $\omega_r$; in addition, the
                  control gain in the OFS should be larger than that in Scheme #2. Thus, $K_{OFS}$($\omega_r$),
                  which represents the upper straight line between ($\omega_{min}$, 600) and ($\omega_{max}$,
                  300) in Fig. 1(b), is set to be inversely proportional to $\omega_r$. The reason for this is as follows.
                  Larger gains than 600 at $\omega_{min}$ and 300 at $\omega_{max}$ can be set to improve
                  the power-smoothing capability, but attention should be paid on more activations of
                  the pitch-angle controller, thereby causing energy losses. In contrast, smaller gains
                  than 600 at $\omega_{min}$ and 300 at $\omega_{max}$ might cause less power- smoothing
                  capability. In this paper, these values are obtained based on a large number of simulation
                  results considering the activations of the pitch-angle controller and power- smoothing
                  capability; however, $K_{OFS}$($\omega_{min}$) and $K_{OFS}$($\omega_{max}$) can be
                  set to other values depending on the design purposes.
               
               As shown in Fig. 1(b), the control gain for Scheme #1 is irrespective of $\omega_r$ while the control gain
                  in Scheme #2 is monotonously increasing with $\omega_r$. In addition, for Scheme #1
                  and Scheme #2 the same control gains are used in both the UFS and OFS. In contrast,
                  the different control gains for the UFS and OFS are defined. $K_{UFS}$($\omega_r$)
                  is set to be proportional to $\omega_r$, which is similar to that of Scheme #2. Thus,
                  ΔP in the proposed scheme in the UFS is similar to that in Scheme #2. Conversely,
                  $K_{OFS}$($\omega_r$) is set to be much larger than that of Scheme #2. Thus, ΔP in
                  the proposed scheme in the OFS is much larger than that in Scheme #2. Therefore, the
                  power-smoothing capability of the proposed scheme can be improved in the OFS while
                  keeping the power-smoothing capability in the UFS similar to that of Scheme #2.
               
               
                     Fig. 1. Operational features of the proposed power-smoothing scheme
 
             
            
                  3. Model System Layout
               Simulations were carried out to validate the efficacy of the proposed power-smoothing
                  scheme based on an EMTP-RV simulator. As shown in Fig. 2, the model system consists of six synchronous generators, a static load of 220 MW
                  and 79 MVAr, a 330-MW asynchronous motor, and an aggregated DFIG-based WPP. To model
                  the power system that has a low ramping capability, all synchronous generators are
                  modeled as steam turbine generators; the steam turbine governor model is the IEEEG1
                  steam governor model[12]. The droop settings for synchronous generators are set to 5%, which is the typical
                  setting of synchronous generators used in Korea’s power system.
               
               
             
            
                  4. Simulation Results
               This section compares the performance of the power- smoothing schemes under the scenarios
                  by varying wind speeds with the wind penetration levels of 18% and 36%, in which the
                  installed capacity for the WPP are 100 MW for Case 1 and 200 MW for Case 2, respectively.
                  In this paper, the wind penetration level is defined as the installed capacity of
                  a WPP divided by the total load[13].
               
               Fig. 3shows the wind profile of the WPP with the average wind speed of 7.5 m/s ranging from
                  4.7 m/s to 11.3 m/s. The performance of the proposed scheme is investigated compared
                  to those of Scheme #2, Scheme #1 with the fixed gain of 25, and MPPT operation.
               
               
                     Fig. 3. Wind profile of the WPP
 
               
                     4.1 Case 1: Wind Penetration Level = 18%
                  As shown in Fig. 4(a), the frequency fluctuation of the proposed scheme is significantly mitigated more
                     than those of Scheme #1 and Scheme #2 because the proposed scheme significantly smooths
                     the output power fluctuation of the WPP by adjusting $K_{UFS}$($\omega_r$) and $K_{OFS}$($\omega_r$)
                     in the UFS and OFS, respectively (see Fig. 4(b)).
                  
                  As shown in Table. 1, the root mean square (RMS) values of $\Delta f$ in Scheme #1, Scheme #2, and the
                     proposed scheme are 0.039 Hz, 0.029 Hz, and 0.024 Hz, respectively. The maximum frequency
                     deviation ($\Delta f_{max}$) in the proposed scheme is 0.055 Hz, which is less than
                     that of Scheme #1 by 0.045 Hz and less than that of Scheme #2 by 0.008 Hz. In addition,
                     the operating range of $\omega_r$ in the proposed scheme is wider than those in the
                     conventional schemes because the proposed control gain in the OFS is larger than those
                     of the conventional schemes. As a result, the proposed scheme effectively utilizes
                     the rotating masses of a DFIG while performing power smoothing.
                  
                  
                        Fig. 4. Simulation results for Case 1
 
                  
                        
                        
Table 1. Results for Case 1
                     
                     
                        						
                        
                              
                                 
                                    | 
                                       
                                    										
                                     
                                    									
                                   | 
                                 
                                       
                                    										
                                     MPPT operation 
                                    									
                                  | 
                                 
                                       
                                    										
                                     Scheme #1 
                                    									
                                  | 
                                 
                                       
                                    										
                                     Scheme #2 
                                    									
                                  | 
                                 
                                       
                                    										
                                     Proposed scheme 
                                    									
                                  | 
                              
                           
                           
                                 
                                    | 
                                       
                                    										
                                     RMS $\left\{\Delta f\right\}$ (Hz) 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.044 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.039 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.029 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.024 
                                    									
                                  | 
                              
                              
                                    | 
                                       
                                    										
                                     $\Delta f_{max}$ (Hz) 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.107 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.100 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.063 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.055 
                                    									
                                  | 
                              
                              
                                    | 
                                       
                                    										
                                     Operating Range of $\omega_r$ for WPP (p.u.) 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.269 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.300 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.359 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.406 
                                    									
                                  | 
                              
                           
                        
                        					
                     
                     
                   
                
               
                     4.2 Case 2: Wind Penetration Level = 36%
                  As the wind penetration level increases, the system frequency fluctuation becomes
                     more significant; therefore, this subsection investigates the performance of the proposed
                     scheme for a higher wind penetration level.
                  
                  Fig. 5shows the results for Case 2, which is identical to Case 1 except for a higher wind
                     penetration level. As in Case 1, the frequency fluctuation is lessened compared to
                     those of the conventional schemes, as shown in Fig. 5(a).
                  
                  As shown in Table. 2, the RMS value of $\Delta f$ in the proposed scheme is 0.042 Hz, which is less than
                     that in Scheme #1 by 0.050 Hz and less than that in Scheme #2 by 0.015 Hz. $\Delta
                     f_{max}$ in the proposed scheme is 0.075 Hz, which is less than that in Scheme #1
                     by 0.160 Hz and less than that in Scheme #2 by 0.020 Hz. Note that the RMS value of
                     $\Delta f$, in an MPPT operation in Case 2 are approximately three times those in
                     Case 1. However, the RMS value of $\Delta f$ and $\Delta f_{max}$ in the proposed
                     scheme in Case 2 are less than twice those in Case 1. Thus, the proposed scheme can
                     improve the frequency-mitigating capability compared to Scheme #2, Scheme #1, and
                     an MPPT operation for a higher wind penetration level. Further, as in Case 1, the
                     operating range of $\omega_r$ in the proposed scheme is 0.431 p.u., which is the largest
                     among the schemes. As a result, the proposed scheme effectively utilizes the rotating
                     masses while performing power smoothing.
                  
                  
                        Fig. 5. Simulation results for Case 2
 
                  
                        
                        
Table 2. Results for Case 2
                     
                     
                        						
                        
                              
                                 
                                    | 
                                       
                                    										
                                     
                                    									
                                   | 
                                 
                                       
                                    										
                                     MPPT operation 
                                    									
                                  | 
                                 
                                       
                                    										
                                     Scheme #1 
                                    									
                                  | 
                                 
                                       
                                    										
                                     Scheme #2 
                                    									
                                  | 
                                 
                                       
                                    										
                                     Proposed scheme 
                                    									
                                  | 
                              
                           
                           
                                 
                                    | 
                                       
                                    										
                                     RMS $\left\{\Delta f\right\}$ (Hz) 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.122 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.092 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.057 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.042 
                                    									
                                  | 
                              
                              
                                    | 
                                       
                                    										
                                     $\Delta f_{max}$ (Hz) 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.309 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.235 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.095 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.075 
                                    									
                                  | 
                              
                              
                                    | 
                                       
                                    										
                                     Operating Range of $\omega_r$ for WPP (p.u.) 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.259 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.342 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.404 
                                    									
                                  | 
                                 
                                       
                                    										
                                     0.431 
                                    									
                                  | 
                              
                           
                        
                        					
                     
                     
                   
                
             
            
                  5. Conclusions
               This paper proposes a power-smoothing scheme of a DFIG that significantly improves
                  the frequency-mitigating capability in the OFS. The proposed scheme employs an additional
                  control loop based on the frequency deviation operating in conjunction with the MPPT
                  control loop. To improve the frequency-mitigating capability in the OFS, the different
                  control gain is proposed that is set to be inversely proportional to the rotor speed.
               
               The results clearly indicate that the proposed scheme effectively smooths the output
                  power fluctuation of a WPP by employing the different control gain for the OFS, thereby
                  mitigating the frequency fluctuation under various wind penetration levels. In addition,
                  the proposed scheme effectively utilizes the rotating masses of a VSWT while performing
                  power smoothing. Hence, it will help reduce the size of ESSs required to mitigate
                  the system frequency fluctuation caused by varying wind speeds.
               
             
          
         
            
                  Acknowledgements
               This work was supported by the Korea Institute of Energy Technology Evaluation and
                  Planning(KETEP) and the Ministry of Trade, Industry & Energy(MOTIE) of the Republic
                  of Korea (No. 20172010000190).
               
             
            
                  
                     References
                  
                     
                        
                        Bevrani Hassan, 2014, Robust Power System Frequency Control, 2nd ed. New York, Springer,
                           Boston, MA

 
                      
                     
                        
                        Lalor G., Mullane A., O’Malley M., 2005, Frequency control and wind turbine technologies,
                           IEEE Trans. Power Syst., Vol. 20, No. 4, pp. 1905-1913

 
                      
                     
                        
                        Fagan E. , Grimes S. , McArdle J. , Smith P. , Stronge M. , 2005, Grid code provisions
                           for wind generators in Ireland, IEEE PES General Meeting,, Vol. 2, No. , pp. 1241-1247

 
                      
                     
                        
                         EirGrid, 2008, Controllable wind farm power station grid code provisions, EirGrid
                           (Ireland) Grid Code, ver. 3.1

 
                      
                     
                        
                        Francisco Díaz-González, Fernando  D. Bianchi, Andreas  Sumper, Oriol  Gomis-Bellmunt,
                           Mar 2014, Control of a flywheel energy storage system for power smoothing in wind
                           power plants, IEEE Trans. Energy Convers., Vol. 29, No. 1, pp. 204-214

 
                      
                     
                        
                        J.  Pegueroles-Queralt, F.  D. Bianchi, O.  Gomis- Bellmunt, 2015, A power smoothing
                           system based on supercapacitors for renewable distributed generation, IEEE Trans.
                           Ind. Electron., Vol. 62, No. 1, pp. 343-350

 
                      
                     
                        
                        Li X. , Hui  D., Lai X. , 2013, Battery Energy Storage Station (BESS)-Based Smoothing
                           Control of Photovoltaic (PV) and Wind Power Generation Fluctuations, IEEE Trans. Sustain.
                           Energy, Vol. 4, No. 2, pp. 464-473

 
                      
                     
                        
                        I.  D. Margaris, S.  A. Papathanassiou, N.  D. Hatziargyriou, A.  D. Hansen, P.  Sorensen,
                           2012, Frequency control in autonomous power systems with high wind power penetration,
                           IEEE Trans. Sustain. Energy, Vol. 3, No. 2, pp. 189-199

 
                      
                     
                        
                        M.  Wang-hansen, R.  Josefsson, H.  Mehmendovic, 2013, Frequency controlling wind
                           power modeling of control strategies, IEEE Trans. Sustain. Energy, Vol. 4, No. 4,
                           pp. 954-959

 
                      
                     
                        
                        Kim Y., Kang M. , Muljadi E. , Park J. W., Kang Y. C. , 2017, Power smoothing of a
                           variable-speed wind turbine generator in association with the rotor speed–dependent
                           gain, IEEE Trans. Sustain. Energy, Vol. 8, No. 3, pp. 990-999

 
                      
                     
                        
                        A.   M. Howlader, N.  Urasaki, A.  Yona, T.  Senjyu, A.  Y. Saber, 2013, A review
                           of output power smoothing methods for wind energy conversion systems, Renew. Sustain.
                           Energy Rev., Vol. 26, No. , pp. 135-146

 
                      
                     
                        
                        Byerly R. T. , Aanstad O. , Berry D. H. , Dunlop R. D. , Ewart D. N. , Fox B. M. ,
                           Johnson L. H. , Tschappat D. W. , 1973, Dynamic models for steam and hydro turbines
                           in power system studies, IEEE Trans. Power App. Syst., Vol. PAS-92, No. 6, pp. 1904-1915

 
                      
                     
                        
                        T.  Ackermann, 2012, Overview of Integration Studies–Methodologies and Results, in
                           Wind Power in Power Systems, 2nd ed., England: John Wiley & Sons, Ltd.

 
                      
                   
                
             
            저자소개
            
            
            
            
                  Mingguang Li (Mingguang Li)
 
            
            
            		His B.S. degree in electrical engineering and its automation from Northeast Electric
            Power University, China, in 2014. 
            
            		He is currently pursuing his M.S. degree from the Department of Electrical Engineering,
            Chonbuk National University, Korea. 
            
            		His research interest is in the frequency regulation of wind power plants. 
            
            		Tel : 063-270-2391 
            
            		Fax : 063-270-2394 
            
            		E-mail : lmg9156@gmail.com
            
            
            
                  Dejian Yang (Dejian Yang)
 
            
            
            		His B.S. and M.S. degrees from Mudanjiang Normal University, China, in 2013 and
            Chonbuk National University, Jeonju, Korea, in 2016, respectively. 
            
            		He is currently pursuing his Ph.D. degree from the Department of Electrical Engineering,
            Chonbuk National University. 
            
            		He was an assistant researcher at the Wind Energy Grid-Adaptive Technology (WeGAT)
            Research Center, which is supported by the Ministry of Science, ICT and future Planning
            (MSIP), Korea. 
            
            		His research interest includes the frequency support of wind power plants. 
            
            		Tel : 063-270-2391 
            
            		Fax : 063-270-2394 
            
            		E-mail : dejian@jbnu.ac.kr
            
            
            
            
            
            		His B.S., M.S., and Ph.D. degrees in electrical engineering from Seoul National
            University, Korea, in 1991, 1993, and 1997, respectively. 
            
            		From 1999 to 2017, he was a professor with the Department of Electrical Engineering,
            Chonbuk National University, Jeonju, Korea. 
            
            		He was the director of the WeGAT Research Center supported by the MSIP, Korea. 
            
            		Since 2017, he has joined Gachon University, Kyunggi-do, Korea. 
            
            		His research interests include the development of control and protection techniques
            for wind power plants. 
            
            		Tel : 031-750-8560 
            
            		Fax : 031-750-8571 
            
            		E-mail : augustinekang33@gmail.com
            
            
            
            
            
            		His B.S., M.S., and Ph.D. degrees in electrical engineering from Seoul National
            University, Korea, in 1987, 1989, and 1995, respectively. 
            
            		He has been a professor of the Department of Energy IT at Gachon University, Kyunggi-do,
            Korea, since 1995. 
            
            		His research interests include the smart grid, supergrid, renewable energy, digitalization
            of electric power, and energy policy. 
            
            		Tel : 031-750-5350 
            
            		Fax : 031-750-8571 
            
            		E-mail : hongpa@gachon.ac.kr