# Patent

## Title
System and a method for controlling at least one voltage converter having a plurality of cells in series

## Abstract

A method of controlling at least one voltage converter having a plurality of cells in series, comprising an AC part and a DC part, characterized in that the AC input voltage (Vei) of each cell is determined directly by the use of a high speed current control loop relating to the AC part and a lower speed voltage control loop relating to the cells, the method including choosing the following voltage control law:




  
    Ψ
    i
  
  =
  
    
      
        C
        i
      
      2
    
    ⁢
    
      (
      
        
          
            2
            
              
                R
                pi
              
              ⁢
              
                C
                i
              
            
          
          ⁢
          
            Z
            i
          
        
        +
        
          
            d
            
              d
              t
            
          
          ⁢
          
            Z
            
              i
              ⁢
              
                
              
              ⁢
              _
              ⁢
              
                
              
              ⁢
              ref
            
          
        
        -
        
          K
          
            1
            ⁢
            Zi
          
        
        -
        
          
            K
            
              2
              ⁢
              Zi
            
          
          ⁢
          
            sign
            ⁡
            
              (
              
                E
                Zi
              
              )
            
          
        
      
      )
    
  




where:

    
    
        K1zi and K2zi are positive adjustment gains;
        Ci is the continuous capacitance of the capacitor Ci of each cell;
        Rpi is the losses associated with each cell;
        Zi—ref is the reference value of Zi=(UDCi)2, UDCi; being the direct voltage across the capacitor Ci; and
        EZi is such that EZi=Zi−Zi—ref.

## Background
CROSS REFERENCE TO RELATED APPLICATIONS OR PRIORITY CLAIM
This application is a National Phase of PCT/EP2009/063642, filed Oct. 19,2009, entitled, “SYSTEM AND METHOD FOR CONTROLLING AT LEAST ONE VOLTAGE CONVERTER WITH SEVERAL CELLS IN SERIES”, and which claims priority of French Patent Application No. 08 57130, filed Oct. 21, 2008.
TECHNICAL FIELD
This invention relates to a system and to a method for controlling at least one voltage converter having a plurality of cells in series, which may be a converter of the type known as a chain link static synchronous compensator (or “chain link STATCOM”).
In order to simplify the description given below, it relates by way of example to a voltage converter having a plurality of cells and being of the chain link STATCOM type.
PRESENT STATE OF THE PRIOR ART
In a very high voltage power distribution network, the electricity is generated by alternators in the form of three-phase alternating current (AC), and the voltage is then increased to very high voltage by step-up transformers before being transported on the distribution network. The very high voltage enables the power to be transported over long distances, while lightening structures, and minimizing losses by heating. It does however remain limited by problems of insulation of the equipment, and by electromagnetic radiation effects. The level of the very high voltage that offers a good compromise extends from 400 kilovolts (kV) to 800 kV.
To enable the power to go between a source and a receiver, it is necessary for the voltage at the source to be at a phase offset of angle θ relative to the voltage at the receiver. This angle θ is called the internal angle of the power line, or transmission angle.
If Vs is the voltage on the source end, Vr the voltage on the receiver end, and X1 the purely inductive impedance of the line, the active power P and reactive power Q supplied by the source are given by the following expressions respectively:



  
    P
    =
    
      
        
          
            V
            s
          
          ⁢
          
            V
            r
          
        
        
          X
          ⁢
          
            
          
          ⁢
          1
        
      
      ⁢
      sin
      ⁢
      
        
      
      ⁢
      θ
    
  




  
    Q
    =
    
      
        
          V
          s
          2
        
        -
        
          
            V
            s
          
          ⁢
          
            V
            r
          
          ⁢
          cos
          ⁢
          
            
          
          ⁢
          θ
        
      
      
        X
        ⁢
        
          
        
        ⁢
        1
      
    
  




  
    
      P
      
        ma
        ⁢
        
          
        
        ⁢
        x
      
    
    =
    
      
        
          V
          s
        
        ⁢
        
          V
          r
        
      
      
        X
        ⁢
        
          
        
        ⁢
        1
      
    
  



These expressions show that the active and reactive powers transmitted on an inductive line are a function of the voltages Vs and Vr, the impedance X1, and the transmission angle θ.
In order to augment the power that can be transmitted over the line, it is accordingly possible to proceed in the following three ways:

    
    
        To augment the voltages Vs and Vr. The insulating distance of the installation, and the space called for in dimensioning the installation, at once impose a limitation on this approach. The radiated electromagnetic field is greater. There is therefore an environmental impact to be taken into account. In addition, the equipment is more expensive and is costly to maintain.
        To act on the transmission angle θ. This angle is a function of the active power supplied by the generators. The maximum angle corresponding to Pmax is θ=π/2. For larger angles, the process enters the descending part of the P=f(θ) curve, which is an unstable zone. To work with angles θ that are too large is to risk losing control of the power system, especially in the event of a transient fault—for example grounding of the phases—on the system where return to normal operation gives rise to a transient increase in the transmission angle in order to evacuate the energy that was produced while the fault condition subsisted, and that was unable to be used by the load and was therefore stored as kinetic energy in the rotors of the generators. It is thus imperative that this angle does not exceed the limit of stability.
        To act on the value of the impedance X1, which is caused to diminish by connecting a capacitor in series with the line, thereby compensating for the reactive power generated by the power transmission line itself. For a given transmission angle, the transmissible power rises as the value of the impedance X1 falls.
    
    


The use of flexible alternating current transmission systems (FACTS), which are electronic power devices with a high speed of response, for controlling the active power, opens up new perspectives for better use of power transmission systems. They permit continuous and rapid action on the various parameters of the system such as phase offset, voltage, and impedance. Power transfers are therefore better controlled, and voltages better maintained, thereby enabling margins of stability to be widened, while also moving towards the power transmission lines being used to transfer the largest currents possible, at the limit of the thermal capacity of these high and very high voltage lines.
FACTS can be grouped in two families, namely parallel FACTS and series FACTS, which are as follows:

    
    
        The parallel FACTS consist in particular of: Mechanical Switched Capacitor (MSC) systems; Static Var Compensator (SVC) systems; and Static Synchronous Compensator (STATCOM) systems.
        The series FACTS consist in particular of: Fixed Capacitor (FC) systems; Thyristor Switched Series Capacitor (TCSC) systems; and Static Synchronous Series Compensator (SSSC) systems.
    
    


It is the STATCOM systems that have the reputation of being reactive power controllers of the future; they can be used in the same way as a synchronous compensator (synchronous machine). The STATCOM can provide dynamic compensation of the voltage, increase transient stability, and improve the damping of a power transmission or distribution network.
A STATCOM is an electronic power apparatus connected in parallel, which injects into the power transmission system an alternating current that is controlled in phase and in amplitude. The amplitude of the voltages modulated by the converter (in phase with the system voltages), enables reactive power to be absorbed or supplied. It thus provides dynamic adjustment of the production or consumption of reactive power. It is also capable of damping transient power oscillations.
The main function of a STATCOM is to provide support for the voltage at critical points in a power transmission and distribution network, by injecting a reactive current into the system. In general, the reactive current injected by a STATCOM is proportional to the difference in voltage between the STATCOM and the system itself. In consequence the control of the reactive current depends on control of the voltage by the converter.
A specific STATCOM topology is based on the assembly of a plurality of standard circuits consisting of converters, and referred to herein as cells; this is known as a Chain Link STATCOM. Besides the control of the reactive current supplied by such a converter, the main problem to be resolved is that of distribution of the total voltage over each cell of the converter.
The document referenced [1] at the end of this description gives a model for a chain link STATCOM, together with some simplified control laws for it. However:

    
    
        firstly, having regard to the approximations of modeling and control, the existing control laws do not have a very high degree of robustness (there are parametric variations, harmonics, and so on), nor is dynamic performance good; and
        secondly, there is no general method of the analytic control of a number N of cells connected in series for this type of converter (association of N cells connected as an H bridge, controlled by pulse width modulation (PWM)).
    
    


The object of the invention is to propose a system and a method that make use of a novel control law to give better distribution of the total voltage to each cell of a chain link STATCOM converter.
SUMMARY OF THE INVENTION
The invention provides a control system for at least one voltage converter having a plurality of cells in series, wherein the module and phase components of an input voltage Vei relative to a frame of reference (d, q) are delivered to each cell, comprising an AC part and a DC part, the system being characterized in that it includes a high speed current control loop relating to the AC part that delivers a control voltage Ve, and a lower speed voltage control loop relating to the cells that delivers an active power ψi, such that:




  
    Ψ
    i
  
  =
  
    
      
        C
        i
      
      2
    
    ⁢
    
      (
      
        
          
            2
            
              
                R
                pi
              
              ⁢
              
                C
                i
              
            
          
          ⁢
          
            Z
            i
          
        
        +
        
          
            d
            
              d
              t
            
          
          ⁢
          
            Z
            
              i
              ⁢
              
                
              
              ⁢
              _
              ⁢
              
                
              
              ⁢
              ref
            
          
        
        -
        
          
            K
            
              1
              ⁢
              Z
              ⁢
              
                
              
              ⁢
              i
            
          
          ⁢
          
            E
            Zi
          
        
        -
        
          
            K
            
              2
              ⁢
              Zi
            
          
          ⁢
          
            sign
            ⁡
            
              (
              
                E
                Zi
              
              )
            
          
        
      
      )
    
  




where:

    
    
        K1zi and K2zi are positive adjustment gains;
        Ci is the continuous capacitance of the capacitor Ci of each cell;
        Rpi is the losses associated with each cell;
        Zi—ref is the referenced value of Zi=(UDCi)2, UDCi being the direct voltage across the capacitor Ci; and
        EZi is such that EZi=Zi−Zi—ref,
        in that the link between these two control loops is obtained via the consumption of active power in the cells representing the output of the lower speed control loop, and
        in that it comprises means for obtaining a current reference ie—refd for the high speed control loop by effecting the summation Σψi of the active power consumptions in the N cells with
    
    






  
    
      i
      
        e
        ⁢
        
          
        
        ⁢
        _
        ⁢
        
          
        
        ⁢
        ref
      
      d
    
    =
    
      
        
          
            ∑
            
              i
              =
              1
            
            N
          
          ⁢
          
            Ψ
            i
          
        
        -
        
          
            V
            e
            q
          
          ⁢
          
            i
            
              e
              ⁢
              
                
              
              ⁢
              _
              ⁢
              
                
              
              ⁢
              ref
            
            q
          
        
      
      
        V
        e
        d
      
    
  
  ,




    
    
         and by making use of a phase locked loop on the control voltage Ve.
    
    


The system of the invention may with advantage further include a phase locked loop on the total output voltage component of the AC part. It may with advantage further include a module for error correction in the phase locked loop, such that the total output voltage component of the AC part is given by the following equation:




  
    V
    e
  
  =
  
    
      V
      ePLL
    
    +
    
      
        K
        PLL
      
      ⁢
      
        
          ∑
          
            i
            =
            1
          
          N
        
        ⁢
        
          (
          
            
              ∫
              0
              1
            
            ⁢
            
              
                E
                Zi
              
              ⁢
              
                d
                t
              
            
          
          )
        
      
    
  




where:





  
    K
    PLL
  
  ⁢
  
    
      ∑
      
        i
        =
        1
      
      N
    
    ⁢
    
      (
      
        
          ∫
          0
          1
        
        ⁢
        
          
            E
            Zi
          
          ⁢
          
            d
            t
          
        
      
      )
    
  




    
    
         is a correction term;
        KPLL is the adjustment gain; and
        EZi is the error in tracking the square of the DC voltage.
    
    


It may with advantage include a module for delivering to the cells module and phase components relative to a frame of reference (d, q). It may with advantage comprise a chain link STATCOM converter on each phase of a three-phase power distribution network.
The invention also provides a method of controlling at least one voltage converter having a plurality of cells in series, wherein the module and phase components of an input voltage Vei relative to a frame of reference (d, q) are delivered to each cell, comprising an AC part and a DC part, the method being characterized in that the AC (Vei) input voltage of each cell is determined directly by the use of a high speed current control loop relating to the AC part that delivers a control voltage Ve, and a lower speed voltage control loop relating to the cells that delivers an active power ψi, such that:




  
    Ψ
    i
  
  =
  
    
      
        C
        i
      
      2
    
    ⁢
    
      (
      
        
          
            2
            
              
                R
                pi
              
              ⁢
              
                C
                i
              
            
          
          ⁢
          
            Z
            i
          
        
        +
        
          
            d
            
              d
              t
            
          
          ⁢
          
            Z
            
              i
              ⁢
              
                
              
              ⁢
              _
              ⁢
              
                
              
              ⁢
              ref
            
          
        
        -
        
          
            K
            
              1
              ⁢
              Zi
            
          
          ⁢
          
            E
            Zi
          
        
        -
        
          
            K
            
              2
              ⁢
              Zi
            
          
          ⁢
          
            sign
            ⁡
            
              (
              
                E
                Zi
              
              )
            
          
        
      
      )
    
  




where:

    
    
        K1zi and K2zi are positive adjustment gains;
        Ci is the continuous capacitance of the capacitor Ci of each cell;
        Rpi is the losses associated with each cell;
        Zi—ref is the referenced value of Zi=(UDCi)2, UDCi being the direct voltage across the capacitor Ci; and
        EZi is such that EZi=Zi−Zi—ref,
        in that the link between these two control loops is obtained via the consumption of active power in the cells representing the output of the lower speed control loop, and
        in that a current reference ie—refd is derived for the high speed control loop by effecting the summation Σψi of the active power consumptions in the N cells with
    
    






  
    
      i
      
        e
        ⁢
        
          
        
        ⁢
        _
        ⁢
        
          
        
        ⁢
        ref
      
      d
    
    =
    
      
        
          
            ∑
            
              i
              =
              1
            
            N
          
          ⁢
          
            Ψ
            i
          
        
        -
        
          
            V
            e
            q
          
          ⁢
          
            i
            
              e
              ⁢
              
                
              
              ⁢
              _
              ⁢
              
                
              
              ⁢
              ref
            
            q
          
        
      
      
        V
        e
        d
      
    
  
  ,




    
    
         and by making use of a phase locked loop on the control voltage Ve.
    
    


Advantageously, the method may further include the use of a phase locked loop on the total output voltage component of the AC part.
Use may with advantage also be made of error correction in the phase locked loop, such that the total output voltage component of the AC part is given by the following equation:




  
    V
    e
  
  =
  
    
      V
      ePLL
    
    +
    
      
        K
        PLL
      
      ⁢
      
        
          ∑
          
            i
            =
            1
          
          N
        
        ⁢
        
          (
          
            
              ∫
              0
              1
            
            ⁢
            
              
                E
                Zi
              
              ⁢
              
                d
                t
              
            
          
          )
        
      
    
  




where:





  
    K
    PLL
  
  ⁢
  
    
      ∑
      
        i
        =
        1
      
      N
    
    ⁢
    
      (
      
        
          ∫
          0
          1
        
        ⁢
        
          
            E
            Zi
          
          ⁢
          
            d
            t
          
        
      
      )
    
  




    
    
         is a correction term;
        KPLL is the adjustment gain; and
        EZi is the error in tracking the square of the DC voltage.
    
    


Advantageously, module and phase components relative to a frame of reference (d, q) are delivered to the cells.
Advantageously, a chain link STATCOM converter is used on each phase of a three-phase power distribution network.
The proposed invention rests on the following features:

    
    
        introduction of a dynamic model of the converter, extended to a number N of cells connected in series;
        application of control techniques, on the dynamic model so developed, that are non-linear, robust and generalized to N cells; and
        improvement of the phase locked loop (PLL) proposed so as to improve resolution of the measured voltage and current vectors that are necessary for the control law into their components in a frame of reference (d, q).
    
    


The system and method of the invention therefore have the following advantages:

    
    
        modeling of the converter generalized for N cells connected in series;
        formulation of a general control law that is valid regardless of the number of cells connected in series;
        validity of the control law regardless of the type of electronic power circuit making use of pulse width modulation;
        robustness of the control law against any parametric variations in a chain link STATCOM converter; and
        good dynamic response.
    
    





BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 shows a configuration of a single-phase chain link STATCOM system of the prior art but making use of the method of this invention.
FIG. 2 is a diagram of one of the DC cells that are connected in series in the same configuration of the prior art.
FIG. 3 shows a simplified single-phase chain link STATCOM system making use of the method of this invention.
FIG. 4 shows the voltage and current control system of the invention.
FIG. 5 shows the generation of distributed voltage with the same vector modules (infinite number of solutions) in accordance with the prior art.
FIG. 6 shows one embodiment of the system of the invention.
FIG. 7 shows a three-phase configuration that makes use of the method of this invention.
FIGS. 8A to 18C show various results that demonstrate the performance of the embodiment shown in FIG. 6.



DETAILED DESCRIPTION OF PARTICULAR EMBODIMENTS
A STATCOM converter with a plurality of cells in series, such as that described in the document referenced [2], comprises, for each phase, N cells 12 connected in series. The DC voltage of each of these cells is monitored by controlling the active component of the current. In addition the reactive component of the current is controlled in order to adjust the control of voltage Vres in the power line of the transmission system.
The diagram of such a single phase chain link STATCOM system, with an AC (alternating current) part 13 and a DC (direct current) part 14, is shown in FIG. 1. It includes a high frequency (HF) filter 10 consisting of a resistor R1, an inductor L1 and a capacitor Cf, together with a transformer 11, the impedance of which is represented by the resistance R2 and inductance L2. The capacitance of the capacitor of each cell 12 is given by Ci (i=1, . . . N), while the losses associated with each cell are represented by the impedance Rpi (i=1, . . . N).
Considering the AC part 13 of the system, associated with the HF filter 10 and transformer 11, the following equation gives the voltage on the distribution network:




  
    
      
        
          
            V
            ref
          
          =
          
            
              
                R
                1
              
              ⁢
              i
            
            +
            
              
                L
                1
              
              ⁢
              
                d
                
                  d
                  t
                
              
              ⁢
              i
            
            +
            
              V
              cf
            
          
        
        ⁢
        
          

        
        ⁢
        and
      
    
    
      
        (
        1
        )
      
    
  
  
    
      
        
          i
          e
        
        =
        
          i
          -
          
            i
            cf
          
        
      
    
    
      
        (
        2
        )
      
    
  




The voltage across the capacitor Cf is given by the equation:




  
    
      
        
          
            d
            
              d
              t
            
          
          ⁢
          
            V
            cf
          
        
        =
        
          
            1
            
              C
              f
            
          
          ⁢
          
            i
            cf
          
        
      
    
    
      
        (
        3
        )
      
    
  




then, by substituting equation (2) into equation (3) we get:





  
    
      
        
          
            d
            
              d
              t
            
          
          ⁢
          
            V
            cf
          
        
        =
        
          
            
              1
              
                C
                f
              
            
            ⁢
            i
          
          -
          
            
              1
              
                C
                f
              
            
            ⁢
            
              i
              e
            
          
        
      
    
    
      
        (
        4
        )
      
    
  




It is also possible to get the following relationship:




  
    
      
        
          V
          cf
        
        =
        
          
            
              R
              2
            
            ⁢
            
              i
              e
            
          
          +
          
            
              L
              2
            
            ⁢
            
              d
              
                d
                t
              
            
            ⁢
            
              i
              e
            
          
          +
          
            V
            e
          
        
      
    
    
      
        (
        5
        )
      
    
  




By combining equations (1), (4) and (5), the behavior of the circuit of the AC part 13 in FIG. 1 is given by the following set of differential equations:




  
    
      
        {
        
          
            
              
                
                  
                    d
                    
                      d
                      t
                    
                  
                  ⁢
                  i
                
                =
                
                  
                    
                      -
                      
                        
                          R
                          1
                        
                        
                          L
                          1
                        
                      
                    
                    ⁢
                    i
                  
                  -
                  
                    
                      1
                      
                        L
                        1
                      
                    
                    ⁢
                    
                      V
                      cf
                    
                  
                  +
                  
                    
                      1
                      
                        L
                        1
                      
                    
                    ⁢
                    
                      V
                      res
                    
                  
                
              
            
          
          
            
              
                
                  
                    
                      d
                      
                        d
                        t
                      
                    
                    ⁢
                    
                      V
                      cf
                    
                  
                  =
                  
                    
                      
                        1
                        
                          C
                          f
                        
                      
                      ⁢
                      i
                    
                    -
                    
                      
                        1
                        
                          C
                          f
                        
                      
                      ⁢
                      
                        i
                        e
                      
                    
                  
                
                ⁢
                
                  
                
              
            
          
          
            
              
                
                  
                    d
                    
                      d
                      t
                    
                  
                  ⁢
                  
                    i
                    e
                  
                
                =
                
                  
                    
                      -
                      
                        
                          R
                          2
                        
                        
                          L
                          2
                        
                      
                    
                    ⁢
                    
                      i
                      e
                    
                  
                  +
                  
                    
                      1
                      
                        L
                        2
                      
                    
                    ⁢
                    
                      V
                      cf
                    
                  
                  -
                  
                    
                      1
                      
                        L
                        2
                      
                    
                    ⁢
                    
                      V
                      e
                    
                  
                
              
            
          
        
      
    
    
      
        (
        6
        )
      
    
  




If we now consider the circuit shown in FIG. 2, we get:




  
    
      
        {
        
          
            
              
                
                  
                    
                      d
                      
                        d
                        t
                      
                    
                    ⁢
                    
                      U
                      
                        D
                        ⁢
                        
                          
                        
                        ⁢
                        
                          C
                          i
                        
                      
                    
                  
                  =
                  
                    
                      1
                      
                        C
                        i
                      
                    
                    ⁢
                    
                      i
                      
                        C
                        i
                      
                    
                  
                
                ⁢
                
                  
                
              
            
          
          
            
              
                
                  i
                  
                    C
                    i
                  
                
                =
                
                  
                    
                      i
                      es
                    
                    -
                    
                      i
                      si
                    
                  
                  =
                  
                    
                      i
                      ep
                    
                    -
                    
                      i
                      
                        R
                        i
                      
                    
                    -
                    
                      i
                      si
                    
                  
                
              
            
          
          
            
              
                
                  
                    i
                    
                      R
                      i
                    
                  
                  =
                  
                    
                      
                        U
                        
                          R
                          i
                        
                      
                      
                        R
                        
                          p
                          i
                        
                      
                    
                    =
                    
                      
                        U
                        ei
                      
                      
                        R
                        
                          p
                          i
                        
                      
                    
                  
                
                ⁢
                
                  
                
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




We then get the following expression for the DC voltage across the capacitor Ci of one cell 12, as shown in FIG. 2:




  
    
      
        
          
            d
            
              d
              t
            
          
          ⁢
          
            U
            
              D
              ⁢
              
                
              
              ⁢
              
                C
                i
              
            
          
        
        =
        
          
            
              1
              
                C
                i
              
            
            ⁢
            
              i
              ep
            
          
          -
          
            
              1
              
                C
                i
              
            
            ⁢
            
              i
              si
            
          
          -
          
            
              1
              
                
                  R
                  pi
                
                ⁢
                
                  C
                  i
                
              
            
            ⁢
            
              U
              ei
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




Since no load is connected in the DC part of the system:

isi=0  (9)

So equation (8) becomes:




  
    
      
        
          
            d
            
              d
              t
            
          
          ⁢
          
            U
            
              D
              ⁢
              
                
              
              ⁢
              
                C
                
                  i
                  ⁢
                  
                    
                  
                
              
            
          
        
        =
        
          
            
              1
              
                C
                i
              
            
            ⁢
            
              i
              ep
            
          
          -
          
            
              1
              
                
                  R
                  pi
                
                ⁢
                
                  C
                  i
                
              
            
            ⁢
            
              U
              ei
            
          
        
      
    
    
      
        (
        10
        )
      
    
  




In addition, we have:

UDCi=URi=Uei  (11)

So equation (10) be put in the following form:




  
    
      
        
          
            d
            
              d
              t
            
          
          ⁢
          
            U
            
              D
              ⁢
              
                
              
              ⁢
              
                C
                
                  i
                  ⁢
                  
                    
                  
                
              
            
          
        
        =
        
          
            
              1
              
                C
                i
              
            
            ⁢
            
              i
              ep
            
          
          -
          
            
              1
              
                
                  
                    R
                    pi
                  
                  ⁢
                  
                    C
                    i
                  
                
                ⁢
                
                  
                
              
            
            ⁢
            
              U
              
                D
                ⁢
                
                  
                
                ⁢
                
                  C
                  i
                
              
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




The relationship between the inputs and outputs of AC and DC voltages and currents is given by the following equation:




  
    
      
        
          [
          
            
              
                
                  i
                  e
                
              
            
            
              
                
                  V
                  ei
                
              
            
          
          ]
        
        =
        
          
            [
            
              
                
                  
                    1
                    
                      
                        W
                        Ri
                      
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                  
                
                
                  0
                
              
              
                
                  0
                
                
                  
                    
                      W
                      Ri
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
              
            
            ]
          
          ⁡
          
            [
            
              
                
                  
                    i
                    ep
                  
                
              
              
                
                  
                    U
                    ei
                  
                
              
            
            ]
          
        
      
    
    
      
        (
        13
        )
      
    
  




We are now obtaining a mathematical model that describes the behavior of the electrical diagram seen in FIG. 2, by combining equations (12) and (13) as follows:




  
    
      
        {
        
          
            
              
                
                  
                    d
                    
                      d
                      t
                    
                  
                  ⁢
                  
                    U
                    
                      D
                      ⁢
                      
                        
                      
                      ⁢
                      
                        C
                        i
                      
                    
                  
                
                =
                
                  
                    
                      -
                      
                        1
                        
                          
                            R
                            
                              p
                              i
                            
                          
                          ⁢
                          
                            C
                            i
                          
                        
                      
                    
                    ⁢
                    
                      U
                      
                        D
                        ⁢
                        
                          
                        
                        ⁢
                        
                          C
                          i
                        
                      
                    
                  
                  +
                  
                    
                      1
                      
                        C
                        
                          i
                          ⁢
                          
                            
                          
                        
                      
                    
                    ⁢
                    
                      
                        W
                        Ri
                      
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                    ⁢
                    
                      i
                      e
                    
                  
                
              
            
          
          
            
              
                
                  
                    V
                    ei
                  
                  =
                  
                    
                      
                        W
                        Ri
                      
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                    ⁢
                    
                      U
                      
                        D
                        ⁢
                        
                          
                        
                        ⁢
                        
                          C
                          i
                        
                      
                    
                  
                
                ⁢
                
                  
                
              
            
          
        
      
    
    
      
        (
        14
        )
      
    
  




The model of the system is then given by associating equations (6) and (14) with the following output filter voltage:




  
    
      
        
          V
          e
        
        =
        
          
            
              V
              
                e
                ⁢
                
                  
                
                ⁢
                1
              
            
            +
            
              V
              
                e
                ⁢
                
                  
                
                ⁢
                2
              
            
            +
            ...
            +
            
              V
              eN
            
          
          =
          
            
              ∑
              
                i
                =
                1
              
              N
            
            ⁢
            
              V
              ei
            
          
        
      
    
    
      
        (
        16
        )
      
    
  




So we get:




  
    
      
        {
        
          
            
              
                
                ⁢
                
                  
                    
                      d
                      
                        d
                        t
                      
                    
                    ⁢
                    i
                  
                  =
                  
                    
                      
                        -
                        
                          
                            R
                            1
                          
                          
                            L
                            1
                          
                        
                      
                      ⁢
                      i
                    
                    -
                    
                      
                        1
                        
                          L
                          1
                        
                      
                      ⁢
                      
                        V
                        cf
                      
                    
                    +
                    
                      
                        1
                        
                          L
                          1
                        
                      
                      ⁢
                      
                        V
                        res
                      
                    
                  
                
              
            
          
          
            
              
                
                ⁢
                
                  
                    
                      d
                      
                        d
                        t
                      
                    
                    ⁢
                    
                      V
                      cf
                    
                  
                  =
                  
                    
                      
                        1
                        
                          C
                          f
                        
                      
                      ⁢
                      i
                    
                    -
                    
                      
                        1
                        
                          C
                          f
                        
                      
                      ⁢
                      
                        i
                        e
                      
                    
                  
                
              
            
          
          
            
              
                
                ⁢
                
                  
                    
                      d
                      
                        d
                        t
                      
                    
                    ⁢
                    
                      i
                      e
                    
                  
                  =
                  
                    
                      
                        -
                        
                          
                            R
                            2
                          
                          
                            L
                            2
                          
                        
                      
                      ⁢
                      
                        i
                        e
                      
                    
                    +
                    
                      
                        1
                        
                          L
                          2
                        
                      
                      ⁢
                      
                        V
                        cf
                      
                    
                    -
                    
                      
                        1
                        
                          L
                          2
                        
                      
                      ⁢
                      
                        V
                        e
                      
                    
                  
                
              
            
          
          
            
              
                
                ⁢
                
                  
                    V
                    e
                  
                  =
                  
                    
                      V
                      
                        e
                        ⁢
                        
                          
                        
                        ⁢
                        1
                      
                    
                    +
                    
                      V
                      
                        e
                        ⁢
                        
                          
                        
                        ⁢
                        2
                      
                    
                    +
                    
                      V
                      
                        e
                        ⁢
                        
                          
                        
                        ⁢
                        3
                      
                    
                    +
                    
                      V
                      
                        e
                        ⁢
                        
                          
                        
                        ⁢
                        4
                      
                    
                    +
                    
                      
                        ∑
                        
                          i
                          =
                          1
                        
                        4
                      
                      ⁢
                      
                        V
                        ei
                      
                    
                  
                
              
            
          
          
            
              
                
                ⁢
                
                  
                    
                      V
                      ei
                    
                    =
                    
                      
                        
                          W
                          Ri
                        
                        ⁡
                        
                          (
                          t
                          )
                        
                      
                      ⁢
                      
                        U
                        
                          D
                          ⁢
                          
                            
                          
                          ⁢
                          
                            C
                            i
                          
                        
                      
                    
                  
                  ,
                  
                    i
                    =
                    1
                  
                  ,
                  ...
                  ⁢
                  
                    
                  
                  ,
                  N
                
              
            
          
          
            
              
                
                  
                    d
                    
                      d
                      t
                    
                  
                  ⁢
                  
                    U
                    
                      D
                      ⁢
                      
                        
                      
                      ⁢
                      
                        C
                        i
                      
                    
                  
                
                =
                
                  
                    
                      -
                      
                        1
                        
                          
                            R
                            
                              p
                              i
                            
                          
                          ⁢
                          
                            C
                            i
                          
                        
                      
                    
                    ⁢
                    
                      U
                      
                        D
                        ⁢
                        
                          
                        
                        ⁢
                        
                          C
                          i
                        
                      
                    
                  
                  +
                  
                    
                      1
                      
                        C
                        i
                      
                    
                    ⁢
                    
                      
                        W
                        Ri
                      
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                    ⁢
                    
                      i
                      e
                    
                  
                
              
            
          
        
      
    
    
      
        (
        17
        )
      
    
  




The above model is not very useful if it is desired to control the transfer of active power between the AC part and the DC part. It is therefore proposed to control the square of the value of the voltages, (UDCi)2. We then have the following definition:

Zi=(UDCi)2  (18)

The time derivative of this is:




  
    
      
        
          
            d
            
              d
              t
            
          
          ⁢
          
            Z
            i
          
        
        =
        
          
            
              -
              
                2
                
                  
                    R
                    pi
                  
                  ⁢
                  
                    C
                    i
                  
                
              
            
            ⁢
            
              Z
              i
            
          
          +
          
            
              2
              
                C
                i
              
            
            ⁢
            
              U
              
                D
                ⁢
                
                  
                
                ⁢
                
                  C
                  i
                
              
            
            ⁢
            
              
                W
                Ri
              
              ⁡
              
                (
                t
                )
              
            
            ⁢
            
              i
              e
            
          
        
      
    
    
      
        (
        19
        )
      
    
  




Using the second equation of (14), we get:




  
    
      
        
          
            d
            
              d
              t
            
          
          ⁢
          
            Z
            i
          
        
        =
        
          
            
              -
              
                2
                
                  
                    R
                    pi
                  
                  ⁢
                  
                    C
                    i
                  
                
              
            
            ⁢
            
              Z
              i
            
          
          +
          
            
              2
              
                C
                i
              
            
            ⁢
            
              V
              ei
            
            ⁢
            
              i
              e
            
          
        
      
    
    
      
        (
        20
        )
      
    
  




The active power transfer in the system is represented by:




  
    
      
        {
        
          
            
              
                P
                =
                
                  
                    ∑
                    
                      i
                      =
                      1
                    
                    N
                  
                  ⁢
                  
                    P
                    i
                  
                
              
            
          
          
            
              
                
                  P
                  i
                
                =
                
                  
                    V
                    ei
                  
                  ⁢
                  
                    i
                    e
                  
                
              
            
          
        
      
    
    
      
        (
        21
        )
      
    
  




The system of equations (17) may be rewritten as follows:




  
    
      
        
      
    
    
      
        (
        22
        )
      
    
  
  
    
      
        {
        
          
            
              
                
                ⁢
                
                  
                    
                      d
                      
                        d
                        t
                      
                    
                    ⁢
                    i
                  
                  =
                  
                    
                      
                        -
                        
                          
                            R
                            1
                          
                          
                            L
                            1
                          
                        
                      
                      ⁢
                      i
                    
                    -
                    
                      
                        1
                        
                          L
                          1
                        
                      
                      ⁢
                      
                        V
                        cf
                      
                    
                    +
                    
                      
                        1
                        
                          L
                          1
                        
                      
                      ⁢
                      
                        V
                        res
                      
                    
                  
                
              
            
            
              
                (
                e
                )
              
            
          
          
            
              
                
                ⁢
                
                  
                    
                      d
                      
                        d
                        t
                      
                    
                    ⁢
                    
                      V
                      cf
                    
                  
                  =
                  
                    
                      
                        1
                        
                          C
                          f
                        
                      
                      ⁢
                      i
                    
                    -
                    
                      
                        1
                        
                          C
                          f
                        
                      
                      ⁢
                      
                        i
                        e
                      
                    
                  
                
              
            
            
              
                (
                f
                )
              
            
          
          
            
              
                
                ⁢
                
                  
                    
                      d
                      
                        d
                        t
                      
                    
                    ⁢
                    
                      i
                      e
                    
                  
                  =
                  
                    
                      
                        -
                        
                          
                            R
                            2
                          
                          
                            L
                            2
                          
                        
                      
                      ⁢
                      
                        i
                        e
                      
                    
                    +
                    
                      
                        1
                        
                          L
                          2
                        
                      
                      ⁢
                      
                        V
                        cf
                      
                    
                    -
                    
                      
                        1
                        
                          L
                          2
                        
                      
                      ⁢
                      
                        V
                        e
                      
                    
                  
                
              
            
            
              
                (
                a
                )
              
            
          
          
            
              
                
                ⁢
                
                  
                    V
                    e
                    d
                  
                  =
                  
                    
                      ∑
                      
                        i
                        =
                        1
                      
                      N
                    
                    ⁢
                    
                      V
                      ei
                      d
                    
                  
                
              
            
            
              
                (
                b
                )
              
            
          
          
            
              
                
                ⁢
                
                  
                    V
                    e
                    q
                  
                  =
                  
                    
                      ∑
                      
                        i
                        =
                        1
                      
                      N
                    
                    ⁢
                    
                      V
                      ei
                      q
                    
                  
                
              
            
            
              
                (
                c
                )
              
            
          
          
            
              
                
                  
                    
                      
                        d
                        
                          d
                          t
                        
                      
                      ⁢
                      
                        Z
                        i
                      
                    
                    =
                    
                      
                        
                          
                            
                              -
                              
                                2
                                
                                  
                                    R
                                    
                                      p
                                      i
                                    
                                  
                                  ⁢
                                  
                                    C
                                    i
                                  
                                
                              
                            
                            ⁢
                            
                              Z
                              i
                            
                          
                          +
                          
                            
                              2
                              
                                C
                                i
                              
                            
                            ⁢
                            
                              (
                              
                                
                                  
                                    V
                                    ei
                                    d
                                  
                                  ⁢
                                  
                                    i
                                    e
                                    d
                                  
                                
                                +
                                
                                  
                                    V
                                    ei
                                    q
                                  
                                  ⁢
                                  
                                    i
                                    e
                                    q
                                  
                                
                              
                              )
                            
                          
                        
                        →
                        i
                      
                      =
                      1
                    
                  
                  ,
                  ...
                  ⁢
                  
                    
                  
                  ,
                  N
                
                ⁢
                
                  
                
              
            
            
              
                (
                d
                )
              
            
          
        
      
    
    
      
        
      
    
  




Equations (22e), (22f) and 22(a) represent the AC part of the system of the invention, while equations (22b) and (22c) represent the voltage distribution link between the total voltage Ve at the output of the AC part and the AC input on each cell 12. The behavior of each cell 12 is finally given by equation (22d), in which the term (Veidied+Veiqieq) is the active power exchanged.
If we neglect the high frequency filter in the range of working frequencies, the electrical circuit to be controlled, and shown in the diagram of FIG. 1, may be simplified as shown in FIG. 3.
The model for the state of this system is given by the following equations:




  
    
      
        
      
    
    
      
        (
        23
        )
      
    
  
  
    
      
        {
        
          
            
              
                
                ⁢
                
                  
                    
                      d
                      
                        d
                        t
                      
                    
                    ⁢
                    
                      i
                      e
                    
                  
                  =
                  
                    
                      
                        -
                        
                          
                            R
                            
                              2
                              ⁢
                              
                                
                              
                              ⁢
                              T
                            
                          
                          
                            L
                            
                              2
                              ⁢
                              
                                
                              
                              ⁢
                              T
                            
                          
                        
                      
                      ⁢
                      
                        i
                        e
                      
                    
                    +
                    
                      
                        1
                        
                          L
                          
                            2
                            ⁢
                            
                              
                            
                            ⁢
                            T
                          
                        
                      
                      ⁢
                      
                        V
                        cf
                      
                    
                    -
                    
                      
                        1
                        
                          L
                          
                            2
                            ⁢
                            
                              
                            
                            ⁢
                            T
                          
                        
                      
                      ⁢
                      
                        V
                        e
                      
                    
                  
                
              
            
            
              
                (
                a
                )
              
            
          
          
            
              
                
                ⁢
                
                  
                    V
                    e
                    d
                  
                  =
                  
                    
                      ∑
                      
                        i
                        =
                        1
                      
                      N
                    
                    ⁢
                    
                      V
                      ei
                      d
                    
                  
                
              
            
            
              
                (
                b
                )
              
            
          
          
            
              
                
                ⁢
                
                  
                    V
                    e
                    q
                  
                  =
                  
                    
                      ∑
                      
                        i
                        =
                        1
                      
                      N
                    
                    ⁢
                    
                      V
                      ei
                      q
                    
                  
                
              
            
            
              
                (
                c
                )
              
            
          
          
            
              
                
                  
                    
                      
                        d
                        
                          d
                          t
                        
                      
                      ⁢
                      
                        Z
                        i
                      
                    
                    =
                    
                      
                        
                          
                            
                              -
                              
                                2
                                
                                  
                                    R
                                    
                                      p
                                      i
                                    
                                  
                                  ⁢
                                  
                                    C
                                    i
                                  
                                
                              
                            
                            ⁢
                            
                              Z
                              i
                            
                          
                          +
                          
                            
                              2
                              
                                C
                                i
                              
                            
                            ⁢
                            
                              (
                              
                                
                                  
                                    V
                                    ei
                                    d
                                  
                                  ⁢
                                  
                                    i
                                    e
                                    d
                                  
                                
                                +
                                
                                  
                                    V
                                    ei
                                    q
                                  
                                  ⁢
                                  
                                    i
                                    e
                                    q
                                  
                                
                              
                              )
                            
                          
                        
                        →
                        i
                      
                      =
                      1
                    
                  
                  ,
                  ...
                  ⁢
                  
                    
                  
                  ,
                  N
                
                ⁢
                
                  
                
              
            
            
              
                (
                d
                )
              
            
          
        
      
    
    
      
        
      
    
  




The system of equation (23) is used for the whole control synthesis. In the control diagram, it is supposed that the system represents a connection based on the time of two control loops. The first control loop (20), having a rapid dynamic (high speed) given by the equation (23a), relates to the AC part. The second control loop (21), having a slower dynamic (lower speed), and relating to the cells, is given by equation (23d). The link between these two loops 20 and 21 is the voltage as resolved into components in the frame of reference (d, q) as illustrated by equations (23b) and (23c).
A high speed current control circuit 22, and a means 23 for generating a sinusoidal signal, are shown in the high speed control loop 20. A lower speed voltage control circuit 24 is shown in the lower speed control loop 21 that relates to the cells, this circuit delivering the active power ψi.
The link between these two control loops 20 and 21 is obtained via the consumption of active power in the cells representing the output of the lower speed control loop 21.
The current reference ie—refd is derived for the high speed control loop 20 by effecting the summation Σψi (see reference numeral 27 in FIG. 4) of the active power consumptions in the N cells. This is possible by means of the use of two voltage control loops, one high speed and one lower speed, the current reference ie—refd delivered by the lower speed loop being seen as a variable with slow variation (i.e. a constant variable) derived from the faster loop. The interconnected control algorithm is obtained in FIG. 4 by making use of a phase locked loop (PLL) 25 on the control voltage Ve.
If equation (23d) is now considered, the lower speed voltage control circuit 24 is differentiated in a differentiator module 26, using a method in accordance with the Lyapunov theory. If the following quantity is considered:

Ψi=Veidied+Veiqieq  (24)

Equation (23d) can be rewritten as follows:




  
    
      
        
          
            d
            
              d
              t
            
          
          ⁢
          
            Z
            i
          
        
        =
        
          
            
              -
              
                2
                
                  
                    R
                    pi
                  
                  ⁢
                  
                    C
                    i
                  
                
              
            
            ⁢
            
              Z
              i
            
          
          +
          
            
              2
              
                C
                i
              
            
            ⁢
            
              Ψ
              i
            
          
        
      
    
    
      
        (
        25
        )
      
    
  




Let us consider the following Lyapunov function:




  
    
      
        
          V
          Lyap
        
        =
        
          
            1
            2
          
          ⁢
          
            E
            Zi
            2
          
        
      
    
    
      
        (
        26
        )
      
    
  




where Ezi=Zi−Zi—ref is the error in tracking the square of the DC voltage. The differential of the Lyapunov function is given by:





  
    
      
        
          
            d
            
              d
              t
            
          
          ⁢
          
            V
            Lyap
          
        
        =
        
          
            E
            Zi
          
          ⁡
          
            (
            
              
                
                  -
                  
                    2
                    
                      
                        R
                        pi
                      
                      ⁢
                      
                        C
                        i
                      
                    
                  
                
                ⁢
                
                  Z
                  i
                
              
              +
              
                
                  2
                  
                    C
                    i
                  
                
                ⁢
                
                  Ψ
                  i
                
              
              -
              
                
                  d
                  
                    d
                    t
                  
                
                ⁢
                
                  Z
                  i_ref
                
              
            
            )
          
        
      
    
    
      
        (
        27
        )
      
    
  




The above equation is defined as negative by choosing the following voltage control law in accordance with the invention:




  
    
      
        
          Ψ
          i
        
        =
        
          
            
              C
              i
            
            2
          
          ⁢
          
            (
            
              
                
                  2
                  
                    
                      R
                      pi
                    
                    ⁢
                    
                      C
                      i
                    
                  
                
                ⁢
                
                  Z
                  i
                
              
              +
              
                
                  d
                  
                    d
                    t
                  
                
                ⁢
                
                  Z
                  i_ref
                
              
              -
              
                
                  K
                  
                    1
                    ⁢
                    
                      
                    
                    ⁢
                    Zi
                  
                
                ⁢
                
                  E
                  Zi
                
              
              -
              
                
                  K
                  
                    2
                    ⁢
                    
                      
                    
                    ⁢
                    Zi
                  
                
                ⁢
                
                  sign
                  ⁡
                  
                    (
                    
                      E
                      Zi
                    
                    )
                  
                
              
            
            )
          
        
      
    
    
      
        (
        28
        )
      
    
  




where K1zi and K2zi are positive adjustment gains.

In the same way we may consider equation (23a) with the associated Lyapunov equation given by




  
    
      
        
          V
          ie
        
        =
        
          
            1
            2
          
          ⁢
          
            E
            ie
            2
          
        
      
    
    
      
        (
        29
        )
      
    
  




where Eie=ie−ie—ref. In the same way, the control law that causes the tracking error Eie to converge towards zero is given by





  
    
      
        
          V
          e
        
        =
        
          
            L
            
              2
              ⁢
              
                
              
              ⁢
              T
            
          
          ⁡
          
            (
            
              
                
                  k
                  
                    1
                    ⁢
                    
                      _
                      ⁢
                      ie
                    
                  
                
                ⁢
                
                  E
                  ie
                
              
              -
              
                
                  k
                  
                    2
                    ⁢
                    
                      _
                      ⁢
                      ie
                    
                  
                
                ⁢
                
                  sign
                  ⁡
                  
                    (
                    
                      E
                      
                        i
                        e
                      
                    
                    )
                  
                
              
              -
              
                
                  
                    R
                    
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      T
                    
                  
                  
                    L
                    
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      T
                    
                  
                
                ⁢
                
                  i
                  e
                
              
              +
              
                
                  1
                  
                    L
                    
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      T
                    
                  
                
                ⁢
                
                  V
                  res
                
              
              -
              
                
                  d
                  
                    d
                    t
                  
                
                ⁢
                
                  i
                  e_ref
                
              
            
            )
          
        
      
    
    
      
        (
        30
        )
      
    
  




where k1—ie and k2—ie are positive adjustment gains.

From equations (24), (23b) and (23c) we get:




  
    
      
        
          
            ∑
            
              i
              =
              1
            
            4
          
          ⁢
          
            Ψ
            i
          
        
        =
        
          
            
              V
              e
              d
            
            ⁢
            
              i
              e
              d
            
          
          +
          
            
              V
              e
              q
            
            ⁢
            
              i
              e
              q
            
          
        
      
    
    
      
        (
        31
        )
      
    
  




Therefore the current component d is calculated by:




  
    
      
        
          i
          e
          d
        
        =
        
          
            
              
                ∑
                
                  i
                  =
                  1
                
                4
              
              ⁢
              
                Ψ
                i
              
            
            -
            
              
                V
                e
                q
              
              ⁢
              
                i
                e
                q
              
            
          
          
            V
            e
            d
          
        
      
    
    
      
        (
        32
        )
      
    
  




Thus we have, below, the reference component d of the current ie—ref for the high speed control loop:




  
    
      
        
          i
          e_ref
          d
        
        =
        
          
            
              
                ∑
                
                  i
                  =
                  1
                
                4
              
              ⁢
              
                Ψ
                i
              
            
            -
            
              
                V
                e
                q
              
              ⁢
              
                i
                e
                q
              
            
          
          
            V
            e
            d
          
        
      
    
    
      
        (
        33
        )
      
    
  




Given that the current reference ie—refd is given by an external loop, the reference of AC current ie—refd can be calculated as follows:

ie—ref=ie—refd cos(ωt)−ie—refq sin(ωt)  (34)

where ω=2π f, f being the frequency of the power distribution network.

If we suppose that there are only two cells in the system, then there is an infinite number of solutions to generate {right arrow over (V)}e1 and {right arrow over (V)}e2 if the same module is considered for all of the vectors. FIG. 5 accordingly presents three solutions. It would be the same for N cells.
In order to improve the system of the invention, the module and phase components are delivered to the cells relative to a frame of reference (d, q) of each voltage {right arrow over (V)}ei. Considering the total voltage {right arrow over (V)}e=[VedVeq]T, the distributed voltages {right arrow over (V)}ei=[VeidVeiq]T cannot be calculated at once by fixing the same module for each vector and by deducing their corresponding phases. It is necessary to apply an additional constraint to the distributed voltages. There are two possible ways to do this, as follows.
Method 1
The same component d is used for all of the distributed voltages.




  
    
      
        
          V
          ei
          d
        
        =
        
          
            V
            
              e
              ⁢
              
                
              
              ⁢
              1
            
            d
          
          =
          
            
              V
              
                e
                ⁢
                
                  
                
                ⁢
                2
              
              d
            
            =
            
              ...
              =
              
                
                  V
                  eN
                  d
                
                =
                
                  
                    V
                    e
                    d
                  
                  N
                
              
            
          
        
      
    
    
      
        (
        35
        )
      
    
  




From equation (25) we get:




  
    
      
        
          V
          ei
          q
        
        =
        
          
            
              Ψ
              i
            
            -
            
              
                V
                ei
                d
              
              ⁢
              
                i
                e
                d
              
            
          
          
            i
            e
            q
          
        
      
    
    
      
        (
        36
        )
      
    
  




Substituting equations (33) and (35) in equation (36), we get the following distributed voltage component q:




  
    
      
        
          V
          ei
          q
        
        =
        
          
            
              1
              
                i
                e
                q
              
            
            ⁢
            
              Ψ
              i
            
          
          -
          
            
              
                
                  ∑
                  
                    i
                    =
                    1
                  
                  N
                
                ⁢
                
                  Ψ
                  i
                
              
              -
              
                
                  V
                  e
                  d
                
                ⁢
                
                  i
                  e
                  d
                
              
            
            
              Ni
              e
              q
            
          
        
      
    
    
      
        (
        37
        )
      
    
  




The distributed voltages to be supplied to the cells are expressed as:




  
    
      
        
          V
          ei
        
        =
        
          
            [
            
              
                
                  
                    V
                    ei
                    d
                  
                
              
              
                
                  
                    V
                    ei
                    q
                  
                
              
            
            ]
          
          =
          
            [
            
              
                
                  
                    
                      V
                      e
                      d
                    
                    N
                  
                
              
              
                
                  
                    
                      
                        1
                        
                          i
                          e
                          q
                        
                      
                      ⁢
                      
                        Ψ
                        i
                      
                    
                    -
                    
                      
                        
                          
                            ∑
                            
                              i
                              =
                              1
                            
                            N
                          
                          ⁢
                          
                            Ψ
                            i
                          
                        
                        -
                        
                          
                            V
                            e
                            d
                          
                          ⁢
                          
                            i
                            e
                            d
                          
                        
                      
                      
                        Ni
                        e
                        q
                      
                    
                  
                
              
            
            ]
          
        
      
    
    
      
        (
        38
        )
      
    
  




It should be noted that the above relationship is not applicable where ieq=0.
Method 2
The same component q is used for all the voltages.




  
    
      
        
          V
          ei
          q
        
        =
        
          
            V
            
              e
              ⁢
              
                
              
              ⁢
              1
            
            q
          
          =
          
            
              V
              
                e
                ⁢
                
                  
                
                ⁢
                2
              
              q
            
            =
            
              ...
              =
              
                
                  V
                  eN
                  q
                
                =
                
                  
                    V
                    e
                    q
                  
                  N
                
              
            
          
        
      
    
    
      
        (
        39
        )
      
    
  




From equation (25) we get:




  
    
      
        
          V
          ei
          d
        
        =
        
          
            
              Ψ
              i
            
            -
            
              
                V
                ei
                q
              
              ⁢
              
                i
                
                  e
                  ⁢
                  
                    
                  
                
                q
              
            
          
          
            i
            e
            d
          
        
      
    
    
      
        (
        40
        )
      
    
  




Substituting equations (33) and (39) in equation (40), we get:




  
    
      
        
          V
          ei
          d
        
        =
        
          
            
              V
              e
              d
            
            
              
                
                  ∑
                  
                    i
                    =
                    1
                  
                  N
                
                ⁢
                
                  Ψ
                  i
                
              
              -
              
                
                  V
                  e
                  q
                
                ⁢
                
                  i
                  e
                  q
                
              
            
          
          ⁢
          
            (
            
              
                Ψ
                i
              
              -
              
                
                  
                    V
                    e
                    q
                  
                  ⁢
                  
                    i
                    e
                    q
                  
                
                N
              
            
            )
          
        
      
    
    
      
        (
        41
        )
      
    
  




The distributed voltages to be supplied to the cells are expressed as follows:




  
    
      
        
          V
          ei
        
        =
        
          
            [
            
              
                
                  
                    V
                    ei
                    d
                  
                
              
              
                
                  
                    V
                    ei
                    q
                  
                
              
            
            ]
          
          =
          
            [
            
              
                
                  
                    
                      
                        V
                        e
                        d
                      
                      
                        
                          
                            ∑
                            
                              i
                              =
                              1
                            
                            N
                          
                          ⁢
                          
                            Ψ
                            i
                          
                        
                        -
                        
                          
                            V
                            e
                            q
                          
                          ⁢
                          
                            i
                            e
                            q
                          
                        
                      
                    
                    ⁢
                    
                      (
                      
                        
                          Ψ
                          i
                        
                        -
                        
                          
                            
                              V
                              e
                              q
                            
                            ⁢
                            
                              i
                              
                                e
                                ⁢
                                
                                  
                                
                              
                              q
                            
                          
                          N
                        
                      
                      )
                    
                  
                
              
              
                
                  
                    
                      V
                      e
                      q
                    
                    N
                  
                
              
            
            ]
          
        
      
    
    
      
        (
        42
        )
      
    
  




A complementary improvement may also be made.
If we consider equation (23a), and if we suppose that the components (d, q) of the distributed voltages Vei are not known precisely, but are estimated as {circumflex over (V)}eid and {circumflex over (V)}eiq, where:




  
    
      
        {
        
          
            
              
                
                  
                    V
                    ^
                  
                  ei
                  d
                
                =
                
                  
                    V
                    ei
                    d
                  
                  +
                  
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      V
                      ei
                      d
                    
                  
                
              
            
          
          
            
              
                
                  
                    V
                    ^
                  
                  ei
                  q
                
                =
                
                  
                    V
                    ei
                    q
                  
                  +
                  
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      V
                      ei
                      q
                    
                  
                
              
            
          
        
      
    
    
      
        (
        43
        )
      
    
  




where ΔVeid and ΔVeiq are estimating errors on each component. Substituting equations (24) and (43) in equation (25), we get:





  
    
      
        
          
            d
            
              d
              t
            
          
          ⁢
          
            Z
            i
          
        
        =
        
          
            
              -
              
                2
                
                  
                    R
                    
                      p
                      i
                    
                  
                  ⁢
                  
                    C
                    i
                  
                
              
            
            ⁢
            
              Z
              i
            
          
          +
          
            
              2
              
                C
                i
              
            
            ⁢
            
              Ψ
              i
            
          
          +
          
            
              2
              
                C
                i
              
            
            ⁢
            
              (
              
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    V
                    ei
                    d
                  
                  ⁢
                  
                    i
                    e
                    d
                  
                
                +
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    V
                    ei
                    e
                  
                  ⁢
                  
                    i
                    e
                    q
                  
                
              
              )
            
          
        
      
    
    
      
        (
        44
        )
      
    
  




Since the voltage components q are close to zero by comparison with the voltage components d, the estimating errors can be considered to be around zero ΔVeiq≅0. Equation (44) can then be rewritten as:




  
    
      
        
          
            d
            
              d
              t
            
          
          ⁢
          
            Z
            i
          
        
        =
        
          
            
              -
              
                2
                
                  
                    R
                    
                      p
                      i
                    
                  
                  ⁢
                  
                    C
                    i
                  
                
              
            
            ⁢
            
              Z
              i
            
          
          +
          
            
              2
              
                C
                i
              
            
            ⁢
            
              Ψ
              i
            
          
          +
          
            
              2
              
                C
                i
              
            
            ⁢
            Δ
            ⁢
            
              
            
            ⁢
            
              V
              ei
              d
            
            ⁢
            
              i
              e
              d
            
          
        
      
    
    
      
        (
        45
        )
      
    
  




Using the same procedure as for the low speed synthesis g of the voltage and using the same Lyapunov function, we get the following differential equation:




  
    
      
        
          
            d
            
              d
              t
            
          
          ⁢
          
            V
            Lyap
          
        
        =
        
          
            E
            Zi
          
          ⁡
          
            (
            
              
                
                  -
                  
                    2
                    
                      
                        R
                        pi
                      
                      ⁢
                      
                        C
                        i
                      
                    
                  
                
                ⁢
                
                  Z
                  i
                
              
              +
              
                
                  2
                  
                    C
                    i
                  
                
                ⁢
                
                  Ψ
                  i
                
              
              -
              
                
                  d
                  
                    d
                    t
                  
                
                ⁢
                
                  Z
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    _
                    ⁢
                    
                      
                    
                    ⁢
                    ref
                  
                
              
              +
              
                
                  2
                  
                    C
                    i
                  
                
                ⁢
                Δ
                ⁢
                
                  
                
                ⁢
                
                  V
                  ei
                  d
                
                ⁢
                
                  i
                  e
                  d
                
              
            
            )
          
        
      
    
    
      
        (
        46
        )
      
    
  




The above equation is defined as negative by choosing the following voltage control law and the estimating error of adaptive power control:




  
    
      
        
          Ψ
          i
        
        =
        
          
            
              C
              i
            
            2
          
          ⁢
          
            (
            
              
                
                  2
                  
                    
                      R
                      pi
                    
                    ⁢
                    
                      C
                      i
                    
                  
                
                ⁢
                
                  Z
                  i
                
              
              +
              
                
                  d
                  
                    d
                    t
                  
                
                ⁢
                
                  Z
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    _
                    ⁢
                    
                      
                    
                    ⁢
                    ref
                  
                
              
              -
              
                
                  K
                  
                    1
                    ⁢
                    Zi
                  
                
                ⁢
                
                  E
                  Zi
                
              
              -
              
                
                  K
                  
                    2
                    ⁢
                    Zi
                  
                
                ⁢
                
                  sign
                  ⁡
                  
                    (
                    
                      E
                      Zi
                    
                    )
                  
                
              
            
            )
          
        
      
    
    
      
        (
        47
        )
      
    
  
  
    
      
        
          Δ
          ⁢
          
            
          
          ⁢
          
            V
            ei
            d
          
        
        =
        
          
            k
            PLLi
          
          ⁢
          
            2
            
              
                C
                i
              
              ⁢
              
                i
                e
                d
              
            
          
          ⁢
          
            
              ∫
              0
              t
            
            ⁢
            
              
                E
                Zi
              
              ⁢
              
                d
                t
              
            
          
        
      
    
    
      
        (
        48
        )
      
    
  




where K1zi, K2zi and KPLLi are adjustment gains. Equation (48) can be configured without detriment to the analysis of overall stability, as follows:




  
    
      
        
          Δ
          ⁢
          
            
          
          ⁢
          
            V
            ei
            d
          
        
        =
        
          
            K
            PLLi
          
          ⁢
          
            
              ∫
              0
              t
            
            ⁢
            
              
                E
                Zi
              
              ⁢
              
                d
                t
              
            
          
        
      
    
    
      
        (
        49
        )
      
    
  




where KPLLi are adjustment gains. The above quantity may be added directly to the components d of each distributed voltage corresponding to each cell.

Instead of correcting each distributed voltage by its components d, it is possible to modify the component d of the total voltage directly. The voltage components obtained are then used for calculating the current reference ie—refd and the distributed voltages.
The component d of total voltage is recalculated by using the output d of the PLL loop VePLLd and the correction term ΔVed, to give:




  
    
      
        
          
            
              
                V
                e
                d
              
              =
              
              ⁢
              
                
                  ∑
                  
                    i
                    =
                    1
                  
                  N
                
                ⁢
                
                  (
                  
                    
                      V
                      ei
                      d
                    
                    +
                    
                      Δ
                      ⁢
                      
                        
                      
                      ⁢
                      
                        V
                        ei
                        d
                      
                    
                  
                  )
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    ∑
                    
                      i
                      =
                      1
                    
                    N
                  
                  ⁢
                  
                    V
                    ei
                    d
                  
                
                +
                
                  
                    ∑
                    
                      i
                      =
                      1
                    
                    N
                  
                  ⁢
                  
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      V
                      ei
                      d
                    
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  V
                  ePLL
                  d
                
                +
                
                  
                    ∑
                    
                      i
                      =
                      1
                    
                    N
                  
                  ⁢
                  
                    (
                    
                      
                        K
                        PLLi
                      
                      ⁢
                      
                        
                          ∫
                          0
                          t
                        
                        ⁢
                        
                          
                            E
                            Zi
                          
                          ⁢
                          
                            d
                            t
                          
                        
                      
                    
                    )
                  
                
              
            
          
        
      
    
    
      
        (
        50
        )
      
    
  




If we suppose that the same adaptive gain is used for each cell, we get:




  
    
      
        
          V
          e
          d
        
        =
        
          
            V
            ePLL
            d
          
          +
          
            
              K
              PLL
            
            ⁢
            
              
                ∑
                
                  i
                  =
                  1
                
                N
              
              ⁢
              
                (
                
                  
                    ∫
                    0
                    t
                  
                  ⁢
                  
                    
                      E
                      Zi
                    
                    ⁢
                    
                      d
                      t
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        51
        )
      
    
  




EXAMPLE
In one example of the embodiment with a system having four cells 12, the control system of the invention is shown in FIG. 6. The voltages distributed from the cells 12, with H (four-arm) bridges, are generated using method 1 described above, because method 2 is very sensitive to harmonics and must be used where ieq≅0. FIG. 6 shows the devices already shown in FIG. 4.
The system of the invention, which here has the reference numeral 30, is connected to the input of a pulse width modulation (PWM) generator 31, and is connected at its output to measuring sensors 32 and to a data collecting device 33. Here, the control law 34 includes correction of phase locked loop (PLL) error at 35, and an algorithm 36 for generating distributed voltages [veidveiq]T.
This example enables the performance of the control system of the invention to be measured. In order to show how strongly the control law of the invention performs, various values of capacitors are used for the direct current (DC) bus in the STATCOM system (C1=1.2 CN, C2=CN, C3=CN and C4=0.8 CN), and the control law is deduced with a common nominal value CN of 18 mF for all the capacitors. Three types of result are obtained in this way, as follows.
1. Control Law with Phase Locked Loop (PLL) Correction
FIGS. 8A, 8B and 8C show AC current control with its corresponding components d (in FIG. 8A) and q (in FIG. 8B):

    
    
        At 0.5 second, 1 p.u. of reactive current is demanded;
        At 1.7 second, the control reference is changed from 1 p.u. to −1 p.u. of reactive current.
    
    


FIG. 8C shows the values of current that are calculated and measured. FIGS. 9A, 9B and 9C are a “zoom” version (i.e. a version with the time scale stretched) of a portion of FIGS. 8A, 8B and 8C respectively, around the 0.5 second mark.
The corresponding AC voltage and the control inputs of the DC cells are shown in FIGS. 10A and 10B respectively. FIGS. 11A and 11B are a “zoom” version of FIGS. 10A and 10B, again around the 0.5 second mark.
FIGS. 12A and 12B show, respectively, the DC voltages and the filtered DC voltages with the reference voltage. FIGS. 13A and 13B are a “zoom” version of FIGS. 12A and 12B, once again around the 0.5 second mark.
It may be noticed that the results thus obtained are very good in terms of dynamic (speed of response) and in terms of the strength of the control law concerned.
2. Control Law with no PLL ConnectionU
Making use of this control law gives rise to errors during normal operation for the DC bus voltages. It can be shown that, even if the phase locked loop (PLL) has low performance, the proposed connection of the outputs of the phase locked loop is able to compensate for disturbances and uncertainties in control. FIGS. 14A, 14B and 14C show, respectively, DC voltages, filtered DC voltages with a reference voltage (with no PLL correction), and corresponding AC H-bridge voltages.
3. Comparison with a Prior Art Control Law
In order to measure the performance of the control law of the invention, this law is compared with a control law that is set out in the document referenced [1].
In that document, the control of the cells is worked out around a mean value, and a control correction is added to each cell if the value is different from the reference value. This justifies the dispersions obtained in the DC voltages around the mean value illustrated in FIGS. 16A, 16B and 16C, which can create instability in the control law. Current control, using the technique described in said document, is illustrated in FIG. 15A (current iq) and FIG. 15B (current id).
On the contrary, as illustrated in FIGS. 18A-18C, the DC voltages are well controlled in using the invention control law. The FIG. 18A illustrates DC voltages, the FIG. 18B illustrates filtered DC voltages and the reference voltage, the FIG. 18C illustrates the mean DC voltage and the reference voltage. The corresponding current control is illustrated in FIGS. 17A (current iq) and 17B (current id).
Example of a Three-phase Version
FIG. 7 shows a three-phase chain link STATCOM configuration making use of the method of the invention, and shows single-phase STATCOMS, namely STATCOM 1, STATCOM 2 and STATCOM 3, together with filters 40, 41 and 42 and a computer 45 for receiving the measurements.
REFERENCES


    [1] “Industrial D-STATCOM Chain Link Modeling and Control” by A. Benchaib, E. Lee-Kwet-Sun, J.-L. Thierry and G. de Preville (12th European Conference on Power Electronics and Applications, EPE 2007, Aalborg, Denmark, 2-5 Sep. 2007).
    [2] “Analysis and Performance Investigation of a Cascaded Multilevel STATCOM for Power System Voltage Regulation” by Nitrus Voraphonpiput, Teratam Bunyagul and Somchai Chatratana (Energy for Sustainable Development, March 2006).

## Claims

1. A control system for at least one voltage converter having a plurality of cells in series wherein a module and phase components of an input voltage Vei relative to a frame of reference (d, q) are delivered to each cell, comprising an AC part and a DC part, characterized in that it includes a high speed current control loop relating to the AC part that delivers a control voltage Ve, and a lower speed voltage control loop relating to the cells that delivers an active power ψi, such that:




  
    Ψ
    i
  
  =
  
    
      
        C
        i
      
      2
    
    ⁢
    
      (
      
        
          
            2
            
              
                R
                pi
              
              ⁢
              
                C
                i
              
            
          
          ⁢
          
            Z
            i
          
        
        +
        
          
            d
            
              d
              t
            
          
          ⁢
          
            Z
            
              i
              ⁢
              
                
              
              ⁢
              _
              ⁢
              
                
              
              ⁢
              ref
            
          
        
        -
        
          
            K
            
              1
              ⁢
              Zi
            
          
          ⁢
          
            E
            Zi
          
        
        -
        
          
            K
            
              2
              ⁢
              Zi
            
          
          ⁢
          
            sign
            ⁡
            
              (
              
                E
                Zi
              
              )
            
          
        
      
      )
    
  




where:
K1zi and K2zi are positive adjustment gains;
Ci is a continuous capacitance of a capacitor Ci of each cell;
Rpi is losses associated with each cell;
Zi—ref is a referenced value of Zi=(UDCi)2, UDCi being a direct voltage across the capacitor Ci; and
EZi is such that EZi=Zi−Zi—ref,
in that a link between these two control loops is obtained via a consumption of active power in the cells representing an output of the lower speed control loop,
in that it comprises means for obtaining a current reference ie—refd for the high speed control loop, by effecting a summation Σψi of the active power consumptions in N cells with





  
    
      i
      
        e
        ⁢
        
          
        
        ⁢
        _
        ⁢
        
          
        
        ⁢
        ref
      
      d
    
    =
    
      
        
          
            ∑
            
              i
              =
              1
            
            N
          
          ⁢
          
            Ψ
            i
          
        
        -
        
          
            V
            e
            q
          
          ⁢
          
            i
            
              e
              ⁢
              
                
              
              ⁢
              _
              ⁢
              
                
              
              ⁢
              ref
            
            q
          
        
      
      
        V
        e
        d
      
    
  
  ,



 and by making use of a phase locked loop on the control voltage Ve, and including a module for error correction in the phase locked loop, such that the control voltage Ve, which is a total output voltage component (Ve) of the AC part, is given by the following equation:





  
    V
    e
  
  =
  
    
      V
      ePLL
    
    +
    
      
        K
        PLL
      
      ⁢
      
        
          ∑
          
            i
            =
            1
          
          N
        
        ⁢
        
          (
          
            
              ∫
              0
              1
            
            ⁢
            
              
                E
                Zi
              
              ⁢
              
                d
                t
              
            
          
          )
        
      
    
  




where:
VePLL is an output of the phase locked loop;





  
    K
    PLL
  
  ⁢
  
    
      ∑
      
        i
        =
        1
      
      N
    
    ⁢
    
      (
      
        
          ∫
          0
          1
        
        ⁢
        
          
            E
            Zi
          
          ⁢
          
            d
            t
          
        
      
      )
    
  



 is a correction term;
KPLL is an adjustment gain; and
EZi is an error in tracking a square of the DC voltage.

2. A system according to claim 1, including a phase locked loop on the total output voltage component (Ve) of the AC part.

3. A system according to claim 1 including a module for delivering to the cells module and phase components relative to a frame of reference (d, q).

4. A system according to claim 1, comprising a chain link STATCOM converter on each phase of a three-phase power distribution network.

5. A method of controlling at least one voltage converter having a plurality of cells in series, wherein a module and phase components of an input voltage Vei relative to a frame of reference (d, q) are delivered to each cell, comprising an AC part and a DC part, characterized in that the AC input voltage (Vei) of each cell is determined directly by a use of a high speed current control loop relating to the AC part that delivers a control voltage Ve and a lower speed voltage control loop relating to the cells that delivers an active power ψi, such that:




  
    Ψ
    i
  
  =
  
    
      
        C
        i
      
      2
    
    ⁢
    
      (
      
        
          
            2
            
              
                R
                pi
              
              ⁢
              
                C
                i
              
            
          
          ⁢
          
            Z
            i
          
        
        +
        
          
            d
            
              d
              t
            
          
          ⁢
          
            Z
            
              i
              ⁢
              
                
              
              ⁢
              _
              ⁢
              
                
              
              ⁢
              ref
            
          
        
        -
        
          
            K
            
              1
              ⁢
              Zi
            
          
          ⁢
          
            E
            Zi
          
        
        -
        
          
            K
            
              2
              ⁢
              Zi
            
          
          ⁢
          
            sign
            ⁡
            
              (
              
                E
                Zi
              
              )
            
          
        
      
      )
    
  




where:
K1zi and K2zi are positive adjustment gains;
Ci is a continuous capacitance of a capacitor Ci of each cell;
Rpi is losses associated with each cell;
Zi—ref is a referenced value of Zi=(UDCi)2, UDCi being a direct voltage across the capacitor Ci; and
EZi is such that EZi=Zi−Zi—ref,
in that a link between the two control loops is obtained via a consumption of active power in the cells representing an output of the lower speed control loop, and in that a current reference ie—refd is derived for the high speed control loop, by effecting a summation Σψi of the active power consumptions in N cells with





  
    
      i
      
        e
        ⁢
        
          
        
        ⁢
        _
        ⁢
        
          
        
        ⁢
        ref
      
      d
    
    =
    
      
        
          
            ∑
            
              i
              =
              1
            
            N
          
          ⁢
          
            Ψ
            i
          
        
        -
        
          
            V
            e
            q
          
          ⁢
          
            i
            
              e
              ⁢
              
                
              
              ⁢
              _
              ⁢
              
                
              
              ⁢
              ref
            
            q
          
        
      
      
        V
        e
        d
      
    
  
  ,



 and by making use of a phase locked loop on the control voltage Ve; and wherein use is made of error correction in the phase locked loop, such that the control voltage, which is a total output voltage component (Ve) of the AC part, is given by the following equation:





  
    V
    e
  
  =
  
    
      V
      ePLL
    
    +
    
      
        K
        PLL
      
      ⁢
      
        
          ∑
          
            i
            =
            1
          
          N
        
        ⁢
        
          (
          
            
              ∫
              0
              1
            
            ⁢
            
              
                E
                Zi
              
              ⁢
              
                d
                t
              
            
          
          )
        
      
    
  




where:
VePLL is an output of the phase locked loop;





  
    K
    PLL
  
  ⁢
  
    
      ∑
      
        i
        =
        1
      
      N
    
    ⁢
    
      (
      
        
          ∫
          0
          1
        
        ⁢
        
          
            E
            Zi
          
          ⁢
          
            d
            t
          
        
      
      )
    
  



 is a correction term;
KPLL is an adjustment gain; and
EZi is an error in tracking a square of the DC voltage.

6. A method according to claim 5, including the use of a phase locked loop on the total output voltage component (Ve) of the AC part.

7. A method according to claim 5, wherein a chain link STATCOM converter is used on each phase of a three-phase power distribution network.