# Patent

## Title
Supercapacitor-based grid fault ride-through system

## Abstract

The supercapacitor-based grid fault ride-through system provides a dynamic model of a wind generation system including a Voltage Source Converter (VSC), which functions as a Static Compensator (STATCOM). The power control capability of the STATCOM is extended by incorporating energy stored in a supercapacitor. The system implements a vector control technique based on the decoupling of real and reactive power. Simulation results show that a fixed speed induction generator is capable of withstanding a significant grid fault when aided by the supercapacitor-based grid fault ride-through system. Moreover, the induction generator regains its pre-fault status immediately after the fault is cleared.

## Background
BACKGROUND OF THE INVENTION
1. Field of the Invention
The present invention relates generally to electrical power conversion, and particularly to a supercapacitor-based grid fault ride-through system that enhances grid fault ride-through capability of a wind power generation system.
2. Description of the Related Art
The increased level of wind penetration into the power system has resulted in the revision of Grid Codes for wind generators in many countries. Grid codes issued during the last few years invariably demand that wind farms must withstand voltage dips to a certain percentage of the nominal voltage (down to 0% in some cases) and for a specified duration. Such requirements are known as Fault Ride Through (FRT) or Low Voltage Ride Through (LVRT), and are described by a voltage vs. time characteristic curve.
The FRT requirements also include fast active and reactive power restoration to the pre-fault values after the system voltage returns to normal operation levels. Some codes impose increased reactive power generation by the wind turbines during the disturbance in order to provide voltage support. Nevertheless, the wind parks are required to behave, in principle, like other generators, e.g., conventional generating units.
Induction generators (IGs) are widely used as wind generators since they are relatively inexpensive, rigid, and require low maintenance. Although variable-speed wind turbines (WTs) are getting popularity, statistics show that a significant percentage of the total installations are of fixed-speed WTs with cage-rotor IG. But this kind of generator consumes a large amount of reactive power during normal operating conditions. This consumption jumps sharply upon the occurrence of grid fault due to a large increase in slip. Instability in induction generator terminal voltage emerges from the deficiencies in the supply of demanded reactive power during the transient period.
Solutions that have been suggested to prevent voltage collapse after the faults include dynamic reactive power sources, i.e., a static reactive volt-ampere compensator (SVC) or a static synchronous compensator (STATCOM). But a common STATCOM does not have the capability of active power compensation because the DC capacitor with it is not a bulk storage device. It can only affect the active power flow in the power system indirectly by regulating the voltage at the common point of connection with the transmission line. A STATCOM with a braking resistor or a battery energy storage system (BESS) have also been proposed for the stability enhancement of a large wind farm. To have short-term active and reactive power exchange ability during disturbance, power electronics-based energy capacitor system (ECS) has also been proposed recently. The two most promising short-term energy storage devices suitable for wind energy applications are flywheels and supercapacitors.
Thus, a supercapacitor-based grid fault ride-through system solving the aforementioned problems is desired.
SUMMARY OF THE INVENTION
The supercapacitor-based grid fault ride-through system provides a dynamic model of a wind generation system, including a Static Synchronous Compensator (STATCOM). The power control capability of the STATCOM is extended by incorporating energy stored in a supercapacitor. The system implements a vector control technique based on the decoupling of real and reactive power. Simulation results show that a fixed speed induction generator is capable of withstanding a significant grid fault when aided by the supercapacitor-based grid fault ride-through system. Moreover, the induction generator regains its pre-fault status immediately after the fault is cleared. Supercapacitors have high power density, present good efficiency, and have a much higher life cycle than batteries. In addition, their manufacturing cost is decreasing rapidly.
We consider a two-mass model of a wind-generator's drive-train system, and develop a complete model of a STATCOM/supercapacitor-based (SESS-based) grid fault ride-through system. The real power (P)-reactive power (Q) control of the SESS-based grid fault ride-through system is decoupled. Due to the STATCOM's fast reactive power control capability and the supercapacitor's quick charge and discharge properties, the control strategy of the supercapacitor-based grid fault ride-through system demonstrates very good potential to make a fixed speed wind generator grid code compliant.
These and other features of the present invention will become readily apparent upon further review of the following specification and drawings.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a schematic diagram of a wind power generation system incorporating a supercapacitor-based grid fault ride-through system according to the present invention.
FIG. 2 is a schematic diagram showing a load line model of the supercapacitor-based grid fault ride-through system according to the present invention.
FIG. 3 is a block diagram of a decoupled controller used in the supercapacitor-based grid fault ride-through system according to the present invention.
FIG. 4 is a block diagram of a bidirectional DC-DC control system of the supercapacitor-based grid fault ride-through system according to the present invention.
FIG. 5 is a plot of the low voltage ride through standard set by FERC.
FIG. 6 is a plot showing variations in slip deviation.
FIG. 7 is a plot showing terminal voltage variation.
FIG. 8 is a plot showing torque variation.
FIG. 9 is a plot showing variation in injected reactive power by the controller.
FIG. 10 is a plot showing variation in real power delivered to a grid.
FIG. 11 is a plot showing variation in injected real power by the controller.
FIG. 12 is a plot showing variation in supercapacitor current.
FIG. 13 is a plot showing variation in supercapacitor voltage.



Similar reference characters denote corresponding features consistently throughout the attached drawings.
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
It will be understood that the diagrams in the Figures depicting the supercapacitor-based grid fault ride-through system are exemplary only. The control processes and circuitry defined herein may be embodied in a dedicated electronic device having a microprocessor, microcontroller, digital signal processor, application specific integrated circuit, field programmable gate array, any combination of the aforementioned devices, or other analog or digital device that combines the functionality of the supercapacitor-based grid fault ride-through system on discrete components, a single chip or multiple chips configured to carry out the control processes described herein.
Moreover, the control processes may be embodied in a general purpose computer having the appropriate peripherals attached thereto and software stored on non-transitory computer readable media, such as hard drives, programmable memory chips, floppy disks, USB drives, and the like, that can be loaded into main memory and executed by a processing unit to carry out the functionality and steps of the control processes.
The supercapacitor-based grid fault ride-through system provides a dynamic model of a wind power generation system including a Static Synchronous Shunt Compensator (STATCOM). The power control capability of the STATCOM is extended by incorporating energy stored in a supercapacitor. The system implements a vector control technique based on the decoupling of real and reactive power. Simulation results show that a fixed speed induction generator is capable of withstanding a significant grid fault when aided by the supercapacitor-based grid fault ride-through system. Moreover, the induction generator regains its pre-fault status immediately after the fault is cleared.
FIG. 1 shows the fixed speed wind turbine cage-generator configuration. The configuration includes a horizontal axis wind-generator 80 connected to a power grid through a step-up transformer (Tr1) and a double circuit transmission line 60. A capacitor bank is used at the generator terminal (bus 1) for reactive power compensation of the induction generator at steady state, A voltage source converter (VSC)-DC link capacitor combination 14 (STATCOM) is connected at the generator terminal through a step-down transformer (Tr2). The supercapacitor-based grid fault ride-through system 10 includes a DC-DC converter 12 and supercapacitor Csc. The system 10 is connected across the dc-link capacitor Cdc of the STATCOM 14 to make up the STATCOM/SESS 11. The models for the different components of the wind-generator system are given in the following:
In the wind turbine model, the mechanical power output of a wind turbine is related to the wind speed Vω by:




  
    
      
        
          P
          m
        
        =
        
          
            1
            2
          
          ⁢
          ρ
          ⁢
          
            
          
          ⁢
          A
          ⁢
          
            
          
          ⁢
          
            
              C
              p
            
            ⁡
            
              (
              
                λ
                ,
                β
              
              )
            
          
          ⁢
          
            
              V
              w
              3
            
            .
          
        
      
    
    
      
        (
        1
        )
      
    
  




Here, ρ is the air density and A is the area swept by the turbine blades. The power coefficient Cp(λ,β) depends on both the blade pitch angle β and the tip speed ratio defined as:





  
    
      
        λ
        =
        
          
            
              Ω
              ⁢
              
                
              
              ⁢
              R
            
            
              V
              w
            
          
          .
        
      
    
    
      
        (
        2
        )
      
    
  




where R is wind turbine rotor radius and Ω is the mechanical angular velocity. An expression for Cp is given by:





  
    
      
        
          
            
              C
              p
            
            ⁡
            
              (
              
                λ
                ,
                β
              
              )
            
          
          =
          
            
              0.5176
              ⁢
              
                (
                
                  
                    116
                    
                      λ
                      i
                    
                  
                  -
                  
                    0.4
                    ⁢
                    β
                  
                  -
                  5
                
                )
              
              ⁢
              
                e
                
                  
                    -
                    21
                  
                  
                    λ
                    i
                  
                
              
            
            +
            
              0.00068
              ⁢
              
                
              
              ⁢
              λ
            
          
        
        ⁢
        
          

        
        ⁢
        where
        ⁢
        
          

        
        ⁢
        
          
            1
            
              λ
              i
            
          
          =
          
            
              1
              
                λ
                +
                
                  0.08
                  ⁢
                  β
                
              
            
            -
            
              
                0.035
                
                  
                    β
                    3
                  
                  +
                  1
                
              
              .
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




In the Induction Generator (IG) and Drive-train model the induction generator is modeled as an equivalent voltage source, E′=e′d+je′q, behind the transient impedance, Z′=Rs+jx′. For drive-train, two-mass model is adopted here as much higher inertia wind turbine rotor is connected to the low inertia IG rotor with a relatively soft shaft. Essentially, the whole system is represented through a set of differential-algebraic equations of the form:

{dot over (x)}=f(x,z)

0=g(x,z),  (4)

where x=[e′d, e′q, ωr, ωt, θs] and z=[ids, iqs, νds, νqs, Pm].

The stator voltages and currents are related through:





  
    
      
        
          
            
              
                
                  
                    v
                    ds
                  
                  =
                  
                    
                      
                        -
                        
                          R
                          s
                        
                      
                      ⁢
                      
                        i
                        ds
                      
                    
                    +
                    
                      
                        x
                        ′
                      
                      ⁢
                      
                        i
                        qs
                      
                    
                    +
                    
                      e
                      d
                      ′
                    
                  
                
              
            
            
              
                
                  
                    v
                    qs
                  
                  =
                  
                    
                      
                        -
                        
                          R
                          s
                        
                      
                      ⁢
                      
                        i
                        qs
                      
                    
                    -
                    
                      
                        x
                        ′
                      
                      ⁢
                      
                        i
                        ds
                      
                    
                    +
                    
                      e
                      q
                      ′
                    
                  
                
              
            
          
          }
        
        .
      
    
    
      
        (
        5
        )
      
    
  




The electromagnetic torque is computed as:

Te=e′dids+e′qiqs  (6)

The subscripts d and q stand for direct and quadrature axis values normalized in per unit (p.u.), respectively.

In the transmission line and load model shown in FIG. 2, Vs, and VB represent terminal voltage and infinite bus voltage, respectively, of the wind induction generator IG 80. Z=(R+jX) is the impedance 59 of transmission line 60 of FIG. 1, including the transformer Tr1. The induced grid fault 70 occurs at end of the transmission line distal from the wind induction generator IG 80. The STATCOM current 9 completes the load line model of the system 10. Compensating capacitor and the local load are lumped together and are represented by the admittance term Y11=g11+jb11. Induction generator terminal voltage is given by:

Vs=VB+(R+jX)(Is+Ist−Vs(g11+jb11))  (7)

where Ist is the current injected by the static compensator-supercapacitor energy storage system (SESS-based grid fault ride-through system 11).

Writing currents and voltages of (7) in terms of their d-q components, equations (5) and (7) can be solved for ids and iqs as:




  
    
      
        
          [
          
            
              
                
                  i
                  ds
                
              
            
            
              
                
                  i
                  qs
                
              
            
          
          ]
        
        =
        
          
            
              [
              
                K
                1
              
              ]
            
            ⁡
            
              [
              
                
                  
                    
                      e
                      d
                      ′
                    
                  
                
                
                  
                    
                      e
                      q
                      ′
                    
                  
                
              
              ]
            
          
          +
          
            
              [
              
                K
                2
              
              ]
            
            ⁡
            
              [
              
                
                  
                    
                      i
                      std
                    
                  
                
                
                  
                    
                      i
                      stq
                    
                  
                
              
              ]
            
          
          +
          
            
              [
              
                K
                3
              
              ]
            
            ⁡
            
              [
              
                
                  
                    
                      V
                      Bd
                    
                  
                
                
                  
                    
                      V
                      Bq
                    
                  
                
              
              ]
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




where, K1, K2 and K3 are matrices whose elements are the system's impedance and admittance parameters, respectively.

In the Supercapacitor and STATCOM model 10, the supercapacitor is modeled as an ideal capacitance and an equivalent series resistance (ESR). The ESR is shown in FIG. 1 as Rsc, and accounts for resistive losses in the dielectric, plate material, and electrolytic solution.
The STATCOM 14 is modeled as a controllable voltage source Vst. The dynamic voltage-current relationship of the STATCOM 14 including the converter transformer Tr2 is obtained from:




  
    
      
        
          
            
              L
              st
            
            ⁢
            
              
                d
                
                  I
                  st
                
              
              
                d
                t
              
            
          
          +
          
            
              R
              st
            
            ⁢
            
              I
              st
            
          
        
        =
        
          
            V
            
              st
              ⁢
              
                
              
            
          
          -
          
            V
            s
          
        
      
    
    
      
        (
        9
        )
      
    
  




Here, Rst and Lst are the resistance and inductance of STATCOM 14, including the converter transformer. Equation (9) can be broken-up to a synchronously rotating d-q reference frame and, including the dynamics of the de link voltage Vdc, the STATCOM—Supercapacitor Energy Storage System 11 can be modeled as:




  
    
      
        
          
            d
            
              i
              std
            
          
          
            d
            t
          
        
        =
        
          
            ω
            b
          
          ⁡
          
            (
            
              
                
                  -
                  
                    
                      R
                      st
                    
                    
                      L
                      st
                    
                  
                
                ⁢
                
                  i
                  std
                
              
              +
              
                
                  ω
                  
                    ω
                    b
                  
                
                ⁢
                
                  i
                  stq
                
              
              +
              
                
                  
                    m
                    ⁢
                    
                      
                    
                    ⁢
                    
                      V
                      
                        d
                        ⁢
                        
                          
                        
                        ⁢
                        c
                      
                    
                  
                  
                    L
                    
                      st
                      ⁢
                      
                        
                      
                    
                  
                
                ⁢
                
                  cos
                  ⁡
                  
                    (
                    
                      ψ
                      +
                      
                        θ
                        s
                      
                    
                    )
                  
                
              
              -
              
                
                  
                    V
                    s
                  
                  
                    L
                    st
                  
                
                ⁢
                cos
                ⁢
                
                  
                
                ⁢
                
                  θ
                  s
                
              
            
            )
          
        
      
    
    
      
        (
        10
        )
      
    
  
  
    
      
        
          
            d
            
              i
              stq
            
          
          
            d
            t
          
        
        =
        
          
            ω
            b
          
          ⁡
          
            (
            
              
                
                  -
                  
                    ω
                    
                      ω
                      b
                    
                  
                
                ⁢
                
                  i
                  std
                
              
              -
              
                
                  
                    R
                    st
                  
                  
                    L
                    st
                  
                
                ⁢
                
                  i
                  stq
                
              
              +
              
                
                  
                    m
                    ⁢
                    
                      
                    
                    ⁢
                    
                      V
                      
                        d
                        ⁢
                        
                          
                        
                        ⁢
                        c
                      
                    
                  
                  
                    L
                    st
                  
                
                ⁢
                
                  sin
                  ⁡
                  
                    (
                    
                      ψ
                      +
                      
                        θ
                        s
                      
                    
                    )
                  
                
              
              -
              
                
                  
                    V
                    s
                  
                  
                    L
                    st
                  
                
                ⁢
                sin
                ⁢
                
                  
                
                ⁢
                
                  θ
                  s
                
              
            
            )
          
        
      
    
    
      
        (
        11
        )
      
    
  
  
    
      
        
          
            d
            
              V
              
                d
                ⁢
                
                  
                
                ⁢
                c
              
            
          
          
            d
            t
          
        
        =
        
          
            
              -
              
                m
                
                  C
                  
                    d
                    ⁢
                    
                      
                    
                    ⁢
                    c
                  
                
              
            
            ⁢
            
              (
              
                
                  
                    i
                    std
                  
                  ⁢
                  
                    cos
                    ⁡
                    
                      (
                      
                        ψ
                        +
                        
                          θ
                          s
                        
                      
                      )
                    
                  
                
                +
                
                  
                    i
                    stq
                  
                  ⁢
                  
                    sin
                    ⁡
                    
                      (
                      
                        ψ
                        +
                        
                          θ
                          s
                        
                      
                      )
                    
                  
                
              
              )
            
          
          -
          
            
              
                D
                r
              
              ⁢
              
                I
                sc
              
            
            
              C
              
                d
                ⁢
                
                  
                
                ⁢
                c
              
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




where m and Ψ are the modulation index and phase angle defined by VSC pulse width modulation (PWM), respectively, and Isc is the supercapacitor current, and Dr is the dc-dc converter duty ratio. The supercapacitor voltage Vsc is given by:





  
    
      
        
          
            V
            
              s
              ⁢
              
                
              
              ⁢
              c
            
          
          =
          
            
              E
              
                s
                ⁢
                
                  
                
                ⁢
                c
              
            
            +
            
              
                R
                
                  s
                  ⁢
                  
                    
                  
                  ⁢
                  c
                
              
              ⁢
              
                I
                sc
              
            
          
        
        ;
        
          
            with
            ⁢
            
              
            
            ⁢
            
              I
              sc
            
          
          =
          
            
              C
              sc
            
            ⁢
            
              
                d
                
                  E
                  sc
                
              
              
                d
                t
              
            
          
        
      
    
    
      
        (
        13
        )
      
    
  




The algebraic equations and differential equations from (1) to (12) represent a complete model of a wind generator integrated with a SESS-based grid fault ride-through system.
Regarding control of the STATCOM 14, more specifically, the decoupled P-Q control, the complex power injected by the STATCOM, S=(Vsejθ)Ist*, gives:

Pst=Vs cos θsistd+Vs sin θsistq 

Qst=Vs sin θsistd−Vs cos θsistq  (14)

The decoupled P-Q control strategy is obtained by defining a new set of STATCOM currents:

Istnew=Iste−jθs  (15)

This makes:

Pst=Vsistdnew; Qst=−Vsistqnew  (16)

Pst and Qst can then be controlled independently by istdnew and Istqnew. Equations (10)-(11) can be rewritten as:




  
    
      
        
          
            
              d
              
                d
                t
              
            
            ⁡
            
              [
              
                
                  
                    
                      i
                      std
                      new
                    
                  
                
                
                  
                    
                      i
                      stq
                      new
                    
                  
                
              
              ]
            
          
          =
          
            
              
                
                  ω
                  b
                
                ⁡
                
                  [
                  
                    
                      
                        
                          -
                          
                            
                              R
                              st
                            
                            
                              L
                              st
                            
                          
                        
                      
                      
                        0
                      
                    
                    
                      
                        0
                      
                      
                        
                          -
                          
                            
                              R
                              st
                            
                            
                              L
                              st
                            
                          
                        
                      
                    
                  
                  ]
                
              
              ⁡
              
                [
                
                  
                    
                      
                        i
                        std
                        new
                      
                    
                  
                  
                    
                      
                        i
                        stq
                        new
                      
                    
                  
                
                ]
              
            
            +
            
              [
              
                
                  
                    
                      x
                      1
                    
                  
                
                
                  
                    
                      x
                      2
                    
                  
                
              
              ]
            
          
        
        ⁢
        
          

        
        ⁢
        where
      
    
    
      
        (
        17
        )
      
    
  
  
    
      
        
          [
          
            
              
                
                  x
                  1
                
              
            
            
              
                
                  x
                  2
                
              
            
          
          ]
        
        =
        
          
            ω
            b
          
          ⁡
          
            [
            
              
                
                  
                    
                      
                        ω
                        
                          ω
                          b
                        
                      
                      ⁢
                      
                        i
                        stq
                        new
                      
                    
                    +
                    
                      
                        
                          m
                          ⁢
                          
                            
                          
                          ⁢
                          
                            V
                            
                              d
                              ⁢
                              
                                
                              
                              ⁢
                              c
                            
                          
                        
                        
                          L
                          st
                        
                      
                      ⁢
                      cos
                      ⁢
                      
                        
                      
                      ⁢
                      ψ
                    
                    -
                    
                      
                        V
                        s
                      
                      
                        L
                        st
                      
                    
                  
                
              
              
                
                  
                    
                      
                        -
                        
                          ω
                          
                            ω
                            b
                          
                        
                      
                      ⁢
                      
                        i
                        std
                        new
                      
                    
                    +
                    
                      
                        
                          m
                          ⁢
                          
                            
                          
                          ⁢
                          
                            V
                            
                              d
                              ⁢
                              
                                
                              
                              ⁢
                              c
                            
                          
                        
                        
                          L
                          st
                        
                      
                      ⁢
                      sin
                      ⁢
                      
                        
                      
                      ⁢
                      ψ
                    
                  
                
              
            
            ]
          
        
      
    
    
      
        (
        18
        )
      
    
  




Defining exd=mVdc cos Ψ and exq=mVdc sin Ψ, equation (18) can be rewritten as:




  
    
      
        
          
            e
            
              xd
              ⁢
              
                
              
            
          
          =
          
            
              
                
                  L
                  st
                
                
                  ω
                  b
                
              
              ⁢
              
                (
                
                  
                    x
                    1
                  
                  -
                  
                    ω
                    ⁢
                    
                      
                    
                    ⁢
                    
                      i
                      istq
                      new
                    
                  
                
                )
              
            
            +
            
              V
              m
            
          
        
        ;
        
          
            and
            ⁢
            
              
            
            ⁢
            
              e
              xq
            
          
          =
          
            
              
                L
                st
              
              
                ω
                b
              
            
            ⁢
            
              (
              
                
                  x
                  2
                
                +
                
                  ω
                  ⁢
                  
                    
                  
                  ⁢
                  
                    i
                    istd
                    new
                  
                
              
              )
            
          
        
      
    
    
      
        (
        19
        )
      
    
  
  
    
      
        m
        =
        
          
            
              
                
                  
                    e
                    xd
                    2
                  
                  +
                  
                    e
                    xq
                    2
                  
                
              
              
                V
                
                  d
                  ⁢
                  
                    
                  
                  ⁢
                  c
                
              
            
            ⁢
            
              
            
            ⁢
            and
            ⁢
            
              
            
            ⁢
            Ψ
          
          =
          
            
              tan
              
                -
                1
              
            
            ⁡
            
              (
              
                
                  e
                  xq
                
                
                  e
                  xd
                
              
              )
            
          
        
      
    
    
      
        (
        20
        )
      
    
  




Equation (17) shows that istd and istq(new) respond to x1 and x2, respectively, through a simple first-order transfer function, proving decoupled control of istd and istq. The control rule is thus completed by defining the feedback loops and PI compensation as follows:




  
    
      
        
          x
          1
        
        =
        
          
            
              (
              
                
                  k
                  
                    p
                    ⁢
                    
                      
                    
                    ⁢
                    1
                  
                
                +
                
                  
                    k
                    
                      i
                      ⁢
                      
                        
                      
                      ⁢
                      1
                    
                  
                  s
                
              
              )
            
            ⁢
            
              (
              
                
                  i
                  std
                  ref
                
                -
                
                  i
                  std
                
              
              )
            
            ⁢
            
              
            
            ⁢
            and
            ⁢
            
              
            
            ⁢
            
              x
              2
            
          
          =
          
            
              (
              
                
                  k
                  
                    p
                    ⁢
                    
                      
                    
                    ⁢
                    2
                  
                
                +
                
                  
                    k
                    
                      i
                      ⁢
                      
                        
                      
                      ⁢
                      2
                    
                  
                  s
                
              
              )
            
            ⁢
            
              (
              
                
                  i
                  stq
                  ref
                
                -
                
                  i
                  stq
                
              
              )
            
          
        
      
    
    
      
        (
        21
        )
      
    
  




The inverter voltage vector is controlled through the m and Ψ defined in equation (20). A block diagram of the control scheme 300 is presented in FIG. 3.

With respect to generation of Pref and Qref for the STATCOM, the control objectives include: (i) regulation of voltage at bus 1 by injecting/absorbing required reactive power, and (ii) minimizing real power variations in Pm2. Qref is generated by sensing the deviation of voltage magnitude at bus 1, i.e., (Vsref−Vs) where Vsref is the pre-disturbance voltage magnitude. Similarly, Pref is generated by sensing the deviation of bus voltage angle, i.e., (θsref−θs), where θsref is the pre-disturbance bus voltage angle. PI controllers are used to ensure zero steady-state error. FIG. 3 shows the generation of Pref and Qref. Corresponding istdref and istqref are generated using equation (16).
In the DC-DC converter-control system, a bidirectional DC-DC converter 12 (most clearly shown in FIG. 4) includes switching components S1 and S2 for controlling the charge and discharge of the supercapacitors Csc. Exemplary bidirectional DC-DC converter 12, shown in FIG. 4, is a buck-boost converter and operates in ‘Buck Mode’ to recharge the supercapacitors Csc, whereas ‘Boost Mode’ transfers the stored energy to the DC link, maintaining the DC link voltage Vdc constant during real power delivery. FIG. 4 shows the main circuit of the bidirectional DC-DC converter 12, which includes three energy storage components, i.e., the DC-link capacitor (Cdc), the choke (Lsc) and the supercapacitors (Csc). Among the objectives of the control are: (i) to regulate the DC-link voltage with energy stored in the supercapacitor, and (ii) to keep the supercapacitor voltage Vsc in its desired operating range. These are achieved through the use of DC-link voltage controller 402, and SC current controller 404. The measured DC-link voltage is compared to the reference Vdcref, and the error is taken into the PI voltage controller, which outputs the reference currents Iscref. The voltage reference for the converter Vpwmref is obtained as the sum of SC voltage Vsc and the inductor voltage reference VLref obtained from the current loop PI controller 404. The response of current controller 404 must be fast enough, with good reference tracking capability. Following design procedure known in the art, the gain values of the current controller 404 are calculated.
Simulations were carried out to satisfy Low Voltage Ride Through (LVRT) requirements issued by U.S. Federal Energy Regulatory Commission (FERC). This code stipulates that the wind generator has to remain connected to the power system if the voltage remains at a level greater than 15% of the nominal voltage for a period that does not exceed 0.625 seconds, as shown in plot 500 of FIG. 5. For comparison, no control, Q control only, and both P & Q control using the SESS-based grid fault ride-through system 10 are considered.
A bolted three-phase fault at the grid connection point is simulated. The fault is cleared after 625 ms. The transient responses are shown in plots 600, 700, and 800 of FIGS. 6-8, respectively, for a nominal loading of 75%. They display the changes in generator slip (Δs), stator voltage variations, and variations in electromagnetic torque developed, respectively. The uncontrolled system response is unstable. During the fault, the induction generator rotor accelerates, and the terminal voltage and flux collapse. When the fault is cleared, it requires extra reactive power support for quick recovery of its flux and terminal voltage. In the absence of such support, it fails to recover its terminal voltage, as can be seen in FIG. 7. The slip continues to increase as the decelerating torque becomes negligible. In this situation, a protective system (either over-speed protection or under-voltage protection) temporarily isolates the wind generator from the grid.
Incorporation of the SESS-based grid fault ride-through system 10 changes the situation dramatically. The STATCOM's Q controller quickly brings back the generator terminal voltage by promptly injecting required reactive power (Qst), as can be seen from plot 900 of FIG. 9. This also keeps the rotor slip, which returns to the nominal value after a brief oscillation, under control. The induction generator 80 can continue to supply power to the grid without requiring any shutdown-and-reconnection exercise.
The real merit of the SESS-based grid fault ride-through system 10 can be observed when it is equipped with both P and Q controllers. The P controller cooperates with the Q controller by additionally suppressing the variation of real power. With their combined action, variation of slip becomes very small, as can be seen in plot 600 of FIG. 6. Oscillations in induction generator output power and terminal voltage also are reduced significantly. More importantly, smoothed power is delivered to the grid immediately after the fault is cleared, as can be seen in plot 1000 of FIG. 10. The P controller of SESS-based grid fault ride-through system 10 injects/absorbs required real power (Pst) in opposite phase of real power variation, thereby canceling it. Plot 1100 of FIG. 11 shows the injected real power from the SESS-based grid fault ride-through system. During a fault, the STATCOM controllers 402 and 404 facilitate absorption of real power. Thus, the STATCOM control acts as a decelerating torque for the induction generator 80. Because of this absorption of real power, corresponding variation of supercapacitor current (Isc) and voltage (Vsc) are shown in plots 1200 and 1300 of FIGS. 12 and 13, respectively.
Due to lack of voltage and frequency control capability, a conventional fixed speed induction generator (FSIG) suffers from meeting some of the important grid integration requirements, such as reactive power compensation and fault ride through. The application of the supercapacitor supported STATCOM system 10 achieves a grid compliant FSIG. The system 10 implements a decoupled P-Q control strategy for STATCOM 14. It is observed that wind-driven FSIG, when supported by supercapacitor-based grid fault ride-through system 10 for STATCOM 14, can achieve significant grid fault tolerance capability. The system 10 provides fault ride-through with very small rotor speed variation, constant terminal voltage, and resumption of smooth power delivery to the grid immediately after the fault is cleared. This demonstrates the excellent performance of the supercapacitor-based P-Q control scheme.
It is to be understood that the present invention is not limited to the embodiments described above, but encompasses any and all embodiments within the scope of the following claims.

## Claims

1. A supercapacitor-based grid fault ride-through system, comprising:
a voltage source converter (VSC) including a DC-link input electrical circuit for receiving input electrical power from a power source, a voltage pulse width modulation (PWM) modulation index control input, a phase control input, and an output electrical circuit adapted for connection to a power grid for providing output electrical power to the power grid;
a voltage controller connected to the voltage PWM modulation index control input;
a current controller connected to the phase control input, said voltage controller and said current controller each being a proportional-integral (PI) controller;
means for detecting a fault from the power grid, the means for fault detecting being connected to the voltage and current controllers, wherein said means for detecting a fault from said power grid comprises means for sensing deviation of voltage magnitude at a local bus connected to the power grid and means for sensing deviation of voltage phase angle at the local bus connected to the power grid;
a supercapacitor having an electrical connection to the DC-link input electrical circuit of the voltage source converter, said supercapacitor providing storage of electrical energy for real and reactive power exchange;
a bidirectional DC-DC converter including a pair of switches connected between the electrical connection of the supercapacitor and the input electrical circuit of the voltage source converter, the DC-DC bidirectional converter, via the switch pair, selectively switching the electrical connection of the supercapacitor to the voltage source converter, wherein when a fault is detected by the means for detection, responsively, the voltage and current controllers regulate the voltage source converter during a time period of the fault using the electrical energy stored in the supercapacitor, thereby enabling the voltage source converter to independently provide compensating real power and reactive power to the power grid during the fault,
wherein each said proportional-integral (PI) controller comprises:
means for generating a summation of a VSC output quadrature reference current and an actual VSC output quadrature current, the VSC output quadrature reference current being based on a summation of a pre-fault reference voltage and a during fault actual voltage at the local bus;
means for generating a summation of a VSC output direct reference current and an actual VSC output direct current at the local bus, the VSC output direct reference current being based on a summation of a pre-fault reference phase angle and a during fault actual phase angle at the local bus;
means for generating a quadrature error signal exq and a direct error signal exd based on the quadrature VSC output current summations, the direct VSC output current summations, and an inductive reactance of the VSC; and
means for generating voltage PWM modulation index and phase angle control signals based on the quadrature error signal exq, the direct error signal exd, and a voltage Vdc at the DC-link input electrical circuitry, the PWM modulation index and phase control signals being characterized by relations:





  
    m
    =
    
      
        
          
            
              
                e
                xd
                2
              
              +
              
                e
                xq
                2
              
            
          
          
            V
            
              d
              ⁢
              
                
              
              ⁢
              c
            
          
        
        ⁢
        
          
        
        ⁢
        and
        ⁢
        
          
        
        ⁢
        Ψ
      
      =
      
        
          tan
          
            -
            1
          
        
        ⁡
        
          (
          
            
              e
              xq
            
            
              e
              xd
            
          
          )
        
      
    
  
  ,




where m is the voltage PWM modulation control signal and Ψ is the phase angle control signal; and
means for routing the m and Ψ control signals to the VSC, thereby controlling modulation index and phase angle defined by pulse width modulation of the VSC.

2. The supercapacitor-based grid fault ride-through system according to claim 1, wherein said bi-directional DC-DC converter is a buck-boost converter, said DC-DC converter being operable in buck mode for charging said supercapacitor and in boost mode for transferring stored energy of said supercapacitor to said voltage source converter.

3. The supercapacitor-based grid fault ride-through system according to claim 1, further comprising a wind turbine generator connected to said supercapacitor-based grid fault ride-through system and to the power grid.

4. The supercapacitor-based grid fault ride-through system according to claim 1, further comprising:
a DC-link capacitor connected to said supercapacitor; and
a choke inductor connected to said supercapacitor, the DC-link capacitor and the choke inductor providing energy storage capacity additional to energy storage capacity of said supercapacitor.