Patent ID: 11962150
Assignee: NORTH CHINA ELECTRIC POWER UNIVERSITY
Field: Electrical machinery, apparatus, energy (Electrical engineering)
Classification: CPC H | IPC H

Claim 0:
1. A method for protecting a power system having inverter-interfaced renewable energy sources (IIRESs), wherein the power system includes an inverter, an inductor-capacitor-inductor (LCL) limiter, a control system, and at least one renewable energy source; and the control system includes a current controller including a saturation limiter and a proportional and integral (PI) controller, a phase-locked system, and a low-voltage ride-through (LVRT)control unit, the method comprising:
by using a Park transformation matrix T2r/3s, determining an output voltage u of the inverter in a dq-coordinate system, as, {
     
      
       
        
         
          u
          d
         
         =
         
          
           
            -
            L
           
           ⁢
           
            
             d
             ⁢
             
              i
              d
             
            
            
             d
             ⁢
             t
            
           
          
          -
          
           R
           ⁢
           
            i
            d
           
          
          +
          
           ω
           ⁢
           
            Li
            q
           
          
          +
          
           e
           d
          
         
        
       
      
      
       
        
         
          u
          q
         
         =
         
          
           
            -
            L
           
           ⁢
           
            
             d
             ⁢
             
              i
              q
             
            
            
             d
             ⁢
             t
            
           
          
          -
          
           R
           ⁢
           
            i
            q
           
          
          -
          
           ω
           ⁢
           
            Li
            d
           
          
          +
          
           e
           q
          
         
        
       
      
     
    
   
   
    
     (
     1
     ), wherein:
the ω is an angular frequency detected by a phase-locked loop (PLL) of the phase-locked system in a synchronous reference frame (SRF), uφ is an inverter output voltage, eφ is a terminal voltage, and iφ is a grid current;
φ∈{d, q} indicates variables in a synchronous reference frame, and R and L are a total resistance and inductance of the power filter;
determining a modulated voltage u* of the output voltage u of the inverter in the dq-coordinate system, as, {
     
      
       
        
         
          u
          d
          *
         
         =
         
          
           -
           
            
             m
             ˆ
            
            d
           
          
          +
          
           
            ω
            
             c
             ⁢
             o
             ⁢
             m
            
           
           ⁢
           L
           ⁢
           
            i
            q
           
          
          +
          
           e
           d
          
         
        
       
      
      
       
        
         
          u
          q
          *
         
         =
         
          
           -
           
            
             m
             ˆ
            
            q
           
          
          -
          
           
            ω
            
             c
             ⁢
             o
             ⁢
             m
            
           
           ⁢
           L
           ⁢
           
            i
            d
           
          
          +
          
           e
           q
          
         
        
       
      
     
    
   
   
    
     (
     2
     ), wherein:
the {circumflex over (m)}φ is an output of the saturation limiter of the current controller, ωcom is an angular frequency of the compensation link of the current controller, and uφ* is a modulated voltage;
upon detecting a grid fault, obtaining current references I* in the dq-coordinate system by the LVRT control unit, as, {
     
      
       
        
         
          
           I
           d
           *
          
          ≤
          
           min
           ⁢
           
            {
            
             
              
               
                I
                max
                2
               
               -
               
                I
                q
                
                 *
                 2
                
               
              
             
             ,
              
             
              
               
                2
                ⁢
                
                 P
                 0
                
               
               
                3
                ⁢
                
                 e
                 d
                
               
              
              ⁢
              
               I
               rated
              
             
            
            }
           
          
         
        
       
       
        
         
          
           I
           d
           *
          
          ≥
          
           
            1
            .
            5
           
           ⁢
           
            (
            
             
              0
              .
              9
             
             -
             
              E
              
               p
               ⁢
               o
               ⁢
               s
              
             
            
            )
           
           ⁢
           
            I
            rated
           
          
         
        
       
      
      ;
     
    
   
   
    
     (
     3
     ), wherein:
the Iφ* is a current reference value of a current controller, Imax is a maximal value of the fault current, P0 is a pre-fault active power of the inverter, Epos is a per-unit value of the positive sequence component of the terminal voltage; and
Imax=1.2Irated such that Irated indicates the rated current;
determining the {circumflex over (m)}φ which is the output of the saturation limiter of the current controller, as, m
        ˆ
       
       φ
      
      (
      t
      )
     
     =
     
      {
      
       
        
         
          
           
            
             k
             p
            
            (
            
             
              I
              φ
              *
             
             -
             
              i
              φ
             
            
            )
           
           +
           
            
             k
             i
            
            ⁢
            
             ∫
             
              
               (
               
                
                 I
                 φ
                 *
                
                -
                
                 i
                 φ
                
               
               )
              
              ⁢
              dt
               
             
            
           
          
         
         
          
           ,
           
            
             
              ❘
              "\[LeftBracketingBar]"
             
             
              m
              φ
             
             
              ❘
              "\[RightBracketingBar]"
             
            
            <
            
             M
             0
            
           
          
         
        
        
         
          
           
            
             +
             
              M
              0
             
            
            ⁢
               
            or
           
            
           -
           
            M
            0
           
            
          
         
         
          
           ,
           
            
             
              ❘
              "\[LeftBracketingBar]"
             
             
              m
              φ
             
             
              ❘
              "\[RightBracketingBar]"
             
            
            ≥
            
             M
             0
            
           
          
         
        
       
       ;
      
     
    
   
   
    
     (
     4
     ), wherein the kp and ki are a proportional gain and an integral gain of a current controller, M0 is a saturation limit of a saturation limiter of the current controller, and mφ is an input of the saturation limiter of the current controller;
determining the ω in a time domain which is the angular frequency detected by the SRF-PLL, as, ω
       ⁢
       
        (
        t
        )
       
      
      =
      
       
        ω
        0
       
       +
       
        Δ
        ⁢
        
         θ
         [
         
          
           
            -
            
             
              ρ
              1
              2
             
             
              2
              ⁢
              λ
             
            
           
           ⁢
           
            e
            
             
              -
              
               ρ
               1
              
             
             ⁢
             t
            
           
          
          +
          
           
            
             ρ
             2
             2
            
            
             2
             ⁢
             λ
            
           
           ⁢
           
            e
            
             -
             
              ρ
              
               2
               t
              
             
            
           
          
         
         ]
        
       
      
     
     ;
    
   
   
    
     (
     5
     ), wherein the Δθ is a phase jump value of the terminal voltage, ω0 is a power angular frequency, and constant coefficients λ, ρ1 and ρ2 are determined by PI parameters of a PLL and a terminal voltage;
determining, by combing the (1), (2) and (5), a fault current in a first stage of a transient phase of the grid fault, as, {
     
      
       
        
         
          
           
            i
            d
           
           (
           t
           )
          
          =
          
           
            
             C
             1
            
            ⁢
            
             e
             
              
               -
               
                1
                
                 τ
                 1
                
               
              
              ⁢
              t
             
            
            ⁢
            
             cos
             ⁡
             (
             
              Δω
              ⁢
              t
             
             )
            
           
           -
           
            
             C
             2
            
            ⁢
            
             e
             
              
               -
               
                1
                
                 τ
                 1
                
               
              
              ⁢
              t
             
            
            ⁢
            
             sin
             ⁡
             (
             
              Δω
              ⁢
              t
             
             )
            
           
           +
           
            C
            3
           
          
         
        
       
       
        
         
          
           
            i
            q
           
           (
           t
           )
          
          =
          
           
            
             C
             1
            
            ⁢
            
             e
             
              
               -
               
                1
                
                 τ
                 1
                
               
              
              ⁢
              t
             
            
            ⁢
            
             sin
             ⁡
             (
             
              Δω
              ⁢
              t
             
             )
            
           
           +
           
            
             C
             2
            
            ⁢
            
             e
             
              
               -
               
                1
                
                 τ
                 1
                
               
              
              ⁢
              t
             
            
            ⁢
            
             cos
             ⁡
             (
             
              Δω
              ⁢
              t
             
             )
            
           
           +
           
            C
            4
           
          
         
        
       
      
      ,
      
       
        
         t
         0
        
        ≤
         
        t
        ≤
        
         t
         d
        
       
       ;
      
     
    
   
   
    
     (
     7
     ), wherein the C1-C4 are constant coefficients derived from initial current values and relevant derivative values at a fault time t0, C3 and C4 are related to variant factors: M0, R, L, and Δω=ωcom−ω(t), an attenuation time constant τ1 satisfies 1/τ1=R/L, and td is an end time of a first stage of a transient phase of a grid fault;
determining, by combing the (1), (2) and (5), a fault current in a second stage of the transient phase of the grid fault, as, {
     
      
       
        
         
          
           
            i
            d
           
           (
           t
           )
          
          =
          
           
            
             C
             5
            
            ⁢
            
             e
             
              
               -
               
                1
                
                 τ
                 2
                
               
              
              ⁢
              t
             
            
           
           +
           
            
             C
             6
            
            ⁢
            
             e
             
              
               -
               
                1
                
                 τ
                 3
                
               
              
              ⁢
              t
             
            
           
           +
           
            C
            7
           
          
         
        
       
       
        
         
          
           
            i
            q
           
           (
           t
           )
          
          =
          
           
            
             C
             8
            
            ⁢
            
             e
             
              
               -
               
                1
                
                 τ
                 2
                
               
              
              ⁢
              t
             
            
           
           +
           
            
             C
             9
            
            ⁢
            
             e
             
              
               -
               
                1
                
                 τ
                 3
                
               
              
              ⁢
              t
             
            
           
           +
           
            C
            
             1
             ⁢
             0
            
           
          
         
        
       
      
      ,
      
       
        
         t
         d
        
        ≤
        t
        ≤
        
         t
         q
        
       
       ;
      
     
    
   
   
    
     (
     9
     ), wherein the C5-C10 are constant coefficients derived from a q-axis current and its derivatives at moment td, attenuation time constants τ2 and τ3 satisfy 1/τ2={2R+kp+[kp2−(2ΔωL)2]1/2}/(2L) and 1/τ3={2R+kp−[kp2−(2ΔωL)2]1/2}/(2L), and tq is an end time of a second stage of a transient phase of the grid fault;
determining, by combing the (1), (2) and (5), a fault current in a third stage of the transient phase of the grid fault, as, {
     
      
       
        
         
          
           
            i
            d
           
           (
           t
           )
          
          =
          
           
            
             e
             
              
               α
               ⁢
               t
              
              +
              
               λ
               1
              
             
            
            ⁢
            
             cos
             ⁡
             (
             
              
               β
               ⁢
               t
              
              +
              
               μ
               1
              
             
             )
            
           
           +
           
            
             e
             
              
               γ
               ⁢
               t
              
              +
              
               λ
               2
              
             
            
            ⁢
            
             cos
             ⁡
             (
             
              
               δ
               ⁢
               t
              
              +
              
               μ
               2
              
             
             )
            
           
           +
           
            C
            
             1
             ⁢
             1
            
           
          
         
        
       
       
        
         
          
           
            i
            q
           
           (
           t
           )
          
          =
          
           
            
             e
             
              
               α
               ⁢
               t
              
              +
              
               λ
               1
              
             
            
            ⁢
            
             sin
             ⁡
             (
             
              
               β
               ⁢
               t
              
              +
              
               μ
               1
              
             
             )
            
           
           +
           
            
             e
             
              
               γ
               ⁢
               t
              
              +
              
               λ
               2
              
             
            
            ⁢
            
             sin
             ⁡
             (
             
              
               δ
               ⁢
               t
              
              +
              
               μ
               2
              
             
             )
            
           
           +
           
            C
            
             1
             ⁢
             2
            
           
          
         
        
       
      
      ,
      
       
        t
        ≥
        
         t
         q
        
       
       ;
      
     
    
   
   
    
     (
     14
     ), wherein α, γ, λ1, λ2, β, δ, μ1, μ2, C11, and C12 are constant coefficients and related to variant factors: R, L, kp, ki, Δω, Id* and Iq*;
switching the control system to a fault control mode by tracking the fault currents in the first stage, second stage and third stage of the transient phase of the grid fault, to the current references obtained by LVRT unit; and
evaluating an effectiveness of a power frequency line protection of the power system, and designing control loops of the IIRESs to improve performance of power frequency line protection.