Next Article in Journal
Exploring Precursors of Renewable Energy Portfolio Diversification Using TPB
Previous Article in Journal
Energy-Efficient Train Driving Based on Optimal Control Theory
Previous Article in Special Issue
Assessment of Converter Performance in Hybrid AC-DC Power System under Optimal Power Flow with Minimum Number of DC Link Control Variables
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bilinear Quadratic Feedback Control of Modular Multilevel Converters

by
Guacira Costa de Oliveira
1,2,
Renato Machado Monaro
1,
Gilney Damm
3,*,
Filipe Perez
4 and
Miguel Jimenez Carrizosa
5
1
Polytechnic School, University of São Paulo (USP), Avenida Professor Luciano Gualberto, 380-Butantã, Sao Paulo 05508-010, SP, Brazil
2
Laboratoire des Signaux et Systèmes, CNRS-CentraleSupélec-Université Paris Saclay, 91190 Gif-sur-Yvette, France
3
COSYS-IMSE, Univ Gustave Eiffel, IFSTTAR, F-77447 Marne-la-Vallée, France
4
SuperGrid Institute, 23 Rue de Cyprian, 69100 Villeurbane, France
5
Centro de Electrónica Industrial, UPM, José Gutiérrez Abascal 2, 28006 Madrid, Spain
*
Author to whom correspondence should be addressed.
Energies 2023, 16(18), 6713; https://doi.org/10.3390/en16186713
Submission received: 4 August 2023 / Revised: 13 September 2023 / Accepted: 14 September 2023 / Published: 19 September 2023

Abstract

:
The present paper introduces the formulation and development of a bilinear quadratic control algorithm for Modular Multilevel Converter (MMCs), with a specific emphasis on achieving internal energy stabilization and balance within the converter. A bilinear average model of the MMC is employed, enabling the separation between the DC voltage and the voltage generated by submodules. The algorithm proposed in this study is formulated using bilinear theory and is founded on quadratic feedback control principles. The stability of the suggested controller is scrutinized using a meticulous mathematical approach based on Lyapunov theory. Subsequently, the theoretical findings are assessed using a comprehensive MMC switching model implemented in Matlab Simscape Electrical. The utilization of a phase-shift PWM technique, accompanied by a sorting algorithm, is considered in the study. Additionally, a comparison between the proposed bilinear controller and a standard PI controller is conducted. The outcomes demonstrate that the proposed controller effectively facilitates the regulation of circulating and AC currents, along with managing the internal energy of MMCs. Consequently, this achievement makes a noteworthy contribution to the field, as it introduces an innovative bilinear control approach capable of stabilizing all the state variables of MMCs converters using a single control law.

1. Introduction

The current electrical grid comprises a wide variety of energy generation sources (such as hydro, nuclear, wind, or solar plants) connected to consumer centers, regardless of distance [1]. Nowadays, the scenario has become even more complex because renewable generation is intrinsically intermittent and unpredictable. Consequently, managing power flow throughout the entire grid has become a challenge [2]. In this context, High Voltage Direct Current (HVDC) technology offers technical advantages such as increased transmission capacities, reduced transmission losses, independence from the system’s frequency, and greater cost-effectiveness compared to AC after a break-even distance [3].
The most recent advancement in the context of converters for HVDC is the Modular Multilevel Converter (MMC) topology. This topology combines the advantages of Voltage Source Converters (VSCs) with increasing voltage levels and lower harmonic production [4,5,6]. The three-phase MMC topology consists of two arms per phase: upper (u) and lower (l) arms, referring to the leg of the converter that generates the voltage over and under the neutral point, as depicted in Figure 1. In each arm, there are N submodules (SMs) connected in series, which are switched sequentially to produce N + 1 voltage levels on the Alternating Current (AC) side of the converter. Each SM is a half-bridge with one capacitor that stores energy, as shown in detail in the figure. Moreover, the SM switching produces a harmonic current that circulates between the arms ( i c i r , j ), composed of second-order harmonics as well [7,8]. The MMC has several advantages due to its modular design, including lower switching losses, lower harmonic contents in the output AC voltage, and a smaller station footprint. These advantages give it potential for applications in a wide voltage range [3,9].
Significant progress has recently been made in MMC control [10,11,12,13,14]. Approaches based on PI controllers for energy control, voltage oscillation reduction, and DC bus voltage stabilization were proposed by [15,16]. Moreover, in [17], current and energy control, including energy balancing, of a periodic bilinear time-varying state-space MMC model, is achieved by a periodic linear quadratic regulator. In [18], a discrete-time bilinear model of an MMC is developed and controlled using a sum-of-squares decomposition method. Additionally, [19] proposes an analytical filter to eliminate the intrinsic oscillations in the SM capacitor’s energy and to improve the dynamic response of the system. Similarly, in [20], a proportional resonant controller aims to regulate the positive and negative sequences of the circulating converter current. Controlling the internal dynamics of an MMC is important to ensure the stability and robustness of the system. Some control schemes that do not consider these internal dynamics may still stabilize the system asymptotically thanks to the linearization in the modulation step. However, these control schemes are less robust because they are prone to poorly damped oscillations on the DC side of the converter. Therefore, it is required to control the MMC’s internal dynamics to prevent these issues [9].
In [21], an improved average model is proposed, where a switching function model represents each converter’s arm that accurately includes each SM’s capacitor dynamics. Furthermore, in [22], a continuous model is used, where a variable-voltage source function describes the arms. A commonality among these models is the increase in system order, leading to heightened complexity in the proposed solutions.
The aforementioned control strategies are mainly based on cascaded control of the arm’s energies and the arm’s current control loops. As an alternative to these approaches, the objective of this research is to propose a control strategy capable of individually controlling each of the state variables using a single control law (AC current, circulating current, total energy, and energy distribution among all upper/lower arms). It also aims to provide rigorous proof demonstrating that the proposed control strategy ensures asymptotic stability within the entire region where a Lyapunov function holds (similar to the development in [23] for low-voltage grids). This approach offers the advantage of enabling an explicit computation of the operational region for the proposed controller. The main contributions of this research are summarized as follows:
  • a direct control strategy for all state variables of the converter (SM energies and currents);
  • a Lyapunov-based stability analysis to demonstrate the suitability of the proposed scheme for asymptotically stabilizing the complex dynamics of the MMC using a single control law.
An average model [18] was used to develop the control, while performance validation was carried out using a detailed MMC switching model implemented in Matlab Simscape Electrical. The switching model includes a low-level controller, which effectively achieves SM voltage balancing. This low-level control has implemented a sorting algorithm based on the standard technique described in [7].

2. Model of a Modular Multilevel Converter

An illustration of an MMC is shown in Figure 1. The N SMs connected in an arm produce a voltage equal to v m , j , as shown in detail for upper arm phase a ( v u , a ). It is a three-phase converter and each phase (j = a, b and c) is connected with a leg that consists of upper and lower arms (m = u or l). There are two operational modes for the SMs; mode OFF implies 0 V as output and mode ON entails V S M = V D C / N as output when the system is well-balanced. It is considered that an inductor L is connected in series to each arm of the MMC to limit the short-circuit current and high-frequency harmonics. Additionally, for the average model, the SM switching losses can be represented by an equivalent resistor R. On the AC side, an R L filter is considered per phase ( R c and L c ). Lastly, the voltage at the Direct Current (DC) side is set as V D C / 2 by the converter on the other side of the transmission line, while on the AC side, v f , j represent the voltage at each phase on the Point of Common Coupling (PCC) with frequency ω .
Based on the MMC circuit shown in Figure 1, an average model can be defined for each phase j by applying Kirchhoff’s voltage law as follows:
V D C 2 = v u , j R i u , j L i ˙ u , j + R c i v , j + L c i ˙ v , j v f , j j = a , b and c
V D C 2 = v l , j + R i l , j + L i ˙ l , j + R c i v , j + L c i ˙ v , j v f , j
The converter currents can be rewritten in terms of i c i r , j and i v , j in Equation (3), which are the circulating and AC side currents, respectively [20].
i c i r , j = i u , j + i l , j 2 , i v , j = i u , j + i l , j
Adding (1) and (2) and replacing Equation (3), one obtains the current dynamic (4) for the phase j from the AC side.
i ˙ v , j = R e q L e q i v , j + v u , j v l , j L e q + 2 v f , j L e q
where
R e q = R + 2 R c , L e q = L + 2 L c
Considering Park’s transformation matrix, the system in dq0 reference is obtained. Taking into account that, in the equilibrium of balanced systems, the zero component is equal to zero, Equation (4) can be rewritten as system (5), which describes the AC current state variables that should be controlled as follows:
i ˙ v , d = R e q L e q i v , d + ω · i v , q + v u , d v l , d L e q + 2 v f , d L e q i ˙ v , q = ω · i v , d R e q L e q i v , q + v u , q v l , q L e q + 2 v f , q L e q
By subtracting (2) from (1) and replacing Equation (3), the circulating current dynamic is obtained:
i ˙ c i r , j = R L i c i r , j v u , j + v l , j 2 L + V D C 2 L
In the dq0 frame, (6) can be rewritten as system (7), where v u , 0 + v l , 0 = v d 0 means the DC component in the converter legs. Then, system (7) describes the circulating current dynamics in the dq0 reference that should be controlled.
i ˙ c i r , d = R L i c i r , d + ω · i c i r , q v u , d + v l , d 2 L i ˙ c i r , q = ω · i c i r , d R L i c i r , q v u , q + v l , q 2 L i ˙ c i r , 0 = R L i c i r , 0 v d 0 2 L + V D C 2 L
The operation of an MMC implies voltage changes in the SM’s capacitor. In this way, SMs will be used to synthesize the output AC voltage and, as a consequence, to drive the output currents. Based on this, it is possible to directly state that the voltage on SMs’ capacitor means a certain amount of stored energy per SM. In this paper, the total amount of energy stored in the entire MMC is considered as a state, as well as the difference between the energy stored on all upper and lower arms. In this approach, proposed by [18], using the MMC’s average model, the total stored energy ( W h ) is achieved by the converter power balance, which represents the derivative of the converter’s energy.
W h ˙ = 3 4 v u , d i v , d + 3 2 v u , d i c i r , d 3 4 v u , q i v , q + 3 2 v u , q i c i r , q + 3 4 v l , d i v , d + 3 2 v l , d i c i r , d + 3 4 v l , q i v , q + 3 2 v l , q i c i r , q + 3 i c i r , 0 v d , 0
Similarly, the derivative of the energy balance ( W v ) is defined as the difference between the stored energy on upper and lower SMs, as shown below:
W v ˙ = 3 4 v u , d i v , d + 3 2 v u , d i c i r , d 3 4 v u , q i v , q + 3 2 v u , q i c i r , q 3 4 v l , d i v , d 3 2 v l , d i c i r , d 3 4 v l , q i v , q 3 2 v l , q i c i r , q
Finally, considering the state vector as x = [ i v , d i v , q i c i r , d i c i r , q i c i r , 0 W h W v ] T and the controller input as u = [ v u , d v u , q v l , d v l , q v d 0 ] T , the three-phase MMC can be represented as the following bilinear form:
x ˙ = A x + k = 1 ρ B k u k x + b k u k + z
where, considering ρ = 5 to represent a three-phase MMC,
A = R e q L e q ω 0 0 0 0 0 ω R e q L e q 0 0 0 0 0 0 0 R L ω 0 0 0 0 0 ω R L 0 0 0 0 0 0 0 R L 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
B 1 = 0 0 0 0 0 0 0 0 0 0 3 4 0 3 2 0 0 3 4 0 3 2 0 0
B 2 = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 3 4 0 3 2 0 0 0 0 3 4 0 3 2 0 0 0
B 3 = 0 0 0 0 0 0 0 0 0 0 3 4 0 3 2 0 0 3 4 0 3 2 0 0
B 4 = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 3 4 0 3 2 0 0 0 0 3 4 0 3 2 0 0 0
B 5 = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0
b 1 = 1 L e q 0 1 2 L 0 0 0 0 T
b 2 = 0 1 L e q 0 1 2 L 0 0 0 T
b 3 = 1 L e q 0 1 2 L 0 0 0 0 T
b 4 = 0 1 L e q 0 1 2 L 0 0 0 T
b 5 = 0 0 0 0 1 2 L 0 0 T
z = 2 v f , d L e q 2 v f , q L e q 0 0 V D C 2 L 0 0 T

3. Bilinear Control Theory

MMCs, as the majority of power electronic converters, can be represented by a bilinear model [24], because their states multiply control inputs as in (10). In the following, we consider a bilinear system with state variables x R n and input u R ρ .

3.1. Bilinear System Stability—A General Case

There are several ways to study the stability of bilinear systems [25]; however, a direct way to carry out this task is to apply standard Lyapunov theory (see [26,27]).
In this context, the goal is to find a positive definite Lyapunov function candidate V, whose time derivative will be negative definite. In this way, the candidate V g : R n R > 0 could be defined as a quadratic function as
V g ( x ) = x T P x
with a symmetric positive definite matrix P R n × n , and we can compute the derivative of the Lyapunov function:
V ˙ g ( x ) = ( A x + k = 1 ρ B k u k x + b k u k ) T P x + x T P ( A x + k = 1 ρ B k u k x + b k u k ) = x T [ A T P + P A ] x + 2 x T P k = 1 ρ u k ( B k x + b k )
where
( B k u k x + b k u k ) T P x = x T P ( B k u k x + b k u k )
Under the control law u k = α s [ B k x + b k ] T P x , with a positive constant α s and the further assumptions that for k = 1 , . . . , ρ , it is true that ( B k x + b k ) T P x 0 x 0 and that x T [ A T P + P A ] x 0 (see [28]), Equation (12) will enter a closed loop:
V ˙ g ( x ) ϕ | | x | | 2 2 α s k = 1 ρ ( x T P ( B k x + b k ) ) 2 V ˙ g ( x ) < 0
where the constant ϕ is greater than or equal to zero.
Considering α s > 0 , Lyapunov’s theory assures that the full system is asymptotically stable.

3.2. Bilinear System Stability—Studied Case

Following [29], with the purpose of maintaining the state variable x on a desired operation point x ¯ , we perform a change of variables x = x ˜ + x ¯ , and u = u ˜ + u ¯ .
Applying the change of variables in Equation (10), one obtains (14):
x ˜ ˙ + x ¯ ˙ = A ( x ˜ + x ¯ ) + k = 1 ρ { ( u ˜ k + u ¯ k ) ( B k ( x ˜ + x ¯ ) + b k ) } + z
which can be rewritten into two parts, where the input reference values ( u ¯ k ) on the equilibrium can be defined as
x ¯ ˙ = A x ¯ + k = 1 ρ u ¯ k ( B k x ¯ + b k ) + z = 0
and then
x ˜ ˙ = A x ˜ + k = 1 ρ u ¯ k B k x ˜ + k = 1 ρ ( B k x ˜ + B k x ¯ + b k ) u ˜ k
The evolution of the error states is analyzed by the Lyapunov function and its derivative:
V ( x ˜ ) = x ˜ T P x ˜
V ˙ ( x ˜ ) = x ˜ T ( A ˜ T P + P A ˜ ) x ˜ 1 s t T e r m + + 2 [ k = 1 m u ˜ k ( B k x ˜ + B k x ¯ + b k ) ] T P x ˜ 2 n d T e r m
where A ˜ = A + k = 1 ρ u ¯ k B k
u ˜ k = α k [ B k x ˜ + B k x ¯ + b k ] T P x ˜
u k = u ˜ k + u ¯ k
Remark 1.
Because of the MMC model characteristics (shape of A matrix), not all state variables are present in the first part of Equation (18). In this way, this part of the Lyapunov function derivative would only be negative semi-definite. As a consequence, its 2nd part will provide the remaining states.
Acknowledging these facts, the following procedure based on Lemma 1 of [30] is necessary to design the controller and access the stability properties of the closed loop system.
Matrix A ˜ R n × n has two zero eigenvalues, λ n 1 ( A ˜ ) = 0 and λ n ( A ˜ ) = 0 , such that A ˜ admits the Jordan decomposition A ˜ = U Λ U 1 , where U C n × n is non-singular and Λ is a diagonal Jordan matrix of A ˜ . Considering Λ R : = d i a g ( e ( Λ i ) ) for all i [ 1 , n 2 ] , then since e ( λ i ) < 0 i [ 1 , n 2 ] , there exists a matrix Γ R R n 2 × n 2 , symmetric and positive definite, such that Γ R Λ R + Λ R T Γ R = Φ , where the diagonal matrices Φ R n 2 × n 2 and Γ R n × n are positive definite, and constants Γ 1 , Γ 2 > 0 as follows:
Φ : = Φ 1 0 0 0 0 0 0 Φ n 2
Γ : = 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 Γ 1 0 0 0 0 0 Γ 2
We solve the equation Γ Λ + Λ * Γ = Q , where
Q : = Φ 1 0 0 0 0 0 0 0 0 0 0 Φ n 2 0 0 0 0 0 0 0 0 0 0 0 0
Then, P = U T Γ U 1 is Hermitian positive definite, since Γ is symmetric positive definite and U 1 has full rank (see [30]). That implies that V ( x ) is a real positive definite function. Next, the derivative of the Lyapunov function V ˙ ( x ˜ ) is evaluated as
V ˙ ( x ˜ ) = c = 1 n 2 Φ c | | x ˜ c | | 2 + 2 ( k = 1 ρ α k ( ( B k x ˜ + B k x ¯ + b k ) T P x ˜ ) 2 )
with x ˜ c , Φ c , c = 1 n 2 , the n 2 error states not associated with a zero eigenvalue and their exponential decay rates, and where we can show that for typical values k = 1 ρ α k ( ( B k x ˜ + B k x ¯ + b k ) T P x ˜ ) 2 0 x ˜ 0 , and as a consequence, V ˙ ( x ˜ ) < 0 , x ˜ 0 .
Theorem 1.
The MMC converter described by the system of equations (10), under the bilinear control law described in Equation (19), and with desired equilibrium point satisfying Equation (15), and the assumption that for k = 1 , . . . , ρ , it is true that ( B k x ˜ + b k ) T P x ˜ 0 x ˜ 0 , will be asymptotically stabilized towards this equilibrium point, inside the whole domain where the assumption is valid.
Proof. 
The proof is standard and is obtained by virtue of the positive definite Lyapunov function (17) of the error system, in which closed-loop derivative (24) is negative definite. Then, by Lyapunov theory, the considered equilibrium point is asymptotically stable. □
Remark 2.
It is important to remark that (24) is negative definite for any tuning parameter value α k > 0 (including all α k equal to one single value). These parameters indeed set the speed of convergence of the terms α k ( ( B k x ˜ + B k x ¯ + b k ) T P x ˜ ) 2 0 . The larger those values, the faster the overall system will converge to its equilibrium point, in a trade-off with the control efforts. It is then necessary to proceed with an optimization strategy for this tuning, even though any chosen value will keep the system stable. Such an optimization study will be carried out in future works.
In order to use the bilinear control defined above, the references for the state variables and the steady-state inputs must be derived from the desired outputs. Consequently, Equation (15) is solved to obtain u ¯ , as shown in (25) for c [ 1 , n 2 ] .
u ¯ = γ x ¯ c + ζ
γ = R e q 2 ω L e q 2 R w 0 ω L + L c 2 R e q 2 w L R 0 R e q 2 ω L + L c 2 R ω 0 ω L + L c 2 R e q 2 ω L R 0 0 0 0 0 2 R
ζ = v f d v f q v f d v f q v D C T
Following the AC active and reactive power definitions in Equations (26) and (27) below, the desired values for the AC currents are defined in Equations (28) and (29).
P e = 3 2 v f , d i v , d
Q e = 3 2 v f , d i v , q
x ¯ 1 = i ¯ v , d = 2 P ¯ e 3 v f , d
x ¯ 2 = i ¯ v , q = 2 Q ¯ e 3 v f , d
Since the components d q of circulating current in an MMC do not produce useful power, and in addition, they increase the converter losses, they should be minimized. For this reason, their references are taken as zero:
x ¯ 3 = i ¯ c i r , d = 0 , x ¯ 4 = i ¯ c i r , q = 0
Reference for zero component of the circulating current comes from the equilibrium point ( x ¯ ˙ = 0 ) of Equation (8), replacing (25) and (30) to obtain Equation (31) as follows:
x ¯ 5 = i ¯ c i r , 0 = V D C V D C 2 4 R ( x ¯ 1 u ¯ 1 + x ¯ 2 u ¯ 2 ) 4 R
One of the roots is unfeasible since it brings the MMC to a very high dissipating mode, so the other root is chosen.
The reference to SM energy is linked with SMs’ capacitor voltage, and this can be expressed as in the following equation:
V D C 2 R i c i r , j 2 L i ˙ c i r , j = 2 N V S M
applying Park’s transformation in Equation (32), the equilibrium point for capacitors’ voltage is
V S M = V D C 2 R x ¯ 5 2 N
Then, the stored energy per SM, with capacitor C S M , can be shown as
W S M = 1 2 C S M V S M 2
By Equation (34), considering all active SMs, the total energy of the converter is obtained:
x ¯ 6 = W h ¯ = 6 N W c = 3 C S M 4 N ( V D C 2 R x ¯ 5 ) 2
Because of the symmetric shape of converters, upper and lower arms should have the same stored energy, W c u p and W c l o w , respectively, during steady state, so zero is chosen as the reference for the energy balance:
3 N W c u p 3 N W c l o w = x ¯ 7 = W ¯ v = 0
The state reference and the stability condition defined above fuel the simulation test section.

4. Validation of the Proposed Control Strategy

The dynamic performance of the proposed bilinear control strategy applied to an MMC switching model (shown in Figure 1) was verified by computer simulations using the Matlab Simscape Electrical environment. The system’s parameters and the controller gains are presented in Table 1. The feedback gains were chosen by a trade-off between speed of response and control effort. In this way, the controlled variables could have the speed dynamics of classical power systems. As explained in Remark 2, any positive value chosen for the parameters α k assures the stability of the equilibrium point. An optimization study would be necessary to set optimal values, which will be carried out in future works. In the present case, all parameters were set to the same value, with already good performance (even if not an optimal performance).
The performance of the proposed bilinear control strategy for the MMC was evaluated through step changes in the converter operating point of variables such as active and reactive power ( P e and Q e ), total stored energy ( W h ), and on the energy balance from the upper and lower arms ( W v ). The tests were carried out in three-step changes: for active and reactive power, for total stored energy, and for the energy balance. For comparison purposes, a PI control based on [16] was implemented. The PI control was tuned according to [31]. The diagram in Figure 2 summarizes the implementation of PI (white blocks) and bilinear (black blocks) controllers, while gray blocks are common for both controllers.

4.1. Step Changes on Active and Reactive Power References

The performance of the proposed bilinear control strategy for the MMC was first evaluated through step changes on active ( P e ) and reactive power ( Q e ). Active and reactive power are followed by state variables in Figure 3. The system is initialized, transmitting zero power.
The considered controllers, bilinear, and PI, provide to the states i c i r d and i c i r 0 different references. So, for those two states, it is written Ref. Bilinear to indicate the bilinear controller’s reference and Ref. PI to identify the PI’s references. The remaining states have the same references, indicated as Ref.
The second and third graphics represent direct ( i v , d ) and quadrature ( i v , q ) AC currents. These two state variables have a direct relation with the required power, shown in the first graphic.
At t = [ 0.05 0.15 0.25 0.45 0.55 ] s, the changes in the active power reference occur. In this test, the bilinear and PI controllers were capable of tracking the i v , d reference in less than 4 ms ; meanwhile, state i v , q is lightly affected by the changes in P e . The i v , q follows the change that occurs in reactive power at t = [ 0.1 0.30 0.40 ] s. It can be observed that the controller was capable of bringing i v , q to its reference i ¯ v , q , for all operating points. Furthermore, for both controllers PI and Bilinear, it is possible to see the coupling effect between i v , d and i v , q , while the change in one state implies a transient deviation (from the reference) in the other.
Concerning the state variable i c i r , d , there are distinct references for each controller as a result of each implementation. Under PI, the reference changes from zero for each active power step. However, PI control drives the state variables to their references. In addition to that, Bilinear keeps i c i r , d in the i ¯ c i r , d = 0 . The controllers hold the quadratic component ( i c i r , q ) in the desired reference for most of the operating points. When Q e changes from positive to zero and negative values, the state variables present a slower response under bilinear control (from t = 0.30 s).
Directly related to the amount of DC power that is transferred, changes in i c i r , 0 (6th graph) follow the changes in i v , d (directly linked with AC active power). Here again, the control also performs well and the PI time response is around 100 ms and around 10 ms for bilinear. In the seventh graph, the total converter energy is shown, and it swings for each new power operation change. It is clear that this state is slower than the others, as it takes around 200 ms to track the reference again under bilinear control.
The eighth graph shows the energy balance between the upper and lower arms. Controls keep the state variable around the reference, however with a slow speed of convergence. Under bilinear control, energy balancing ( W v ) requires more than 150 ms to track its reference, while at the same time, under PI control, state variables require different speeds of convergence according to the operation points.
It is important to highlight that Figure 3 summarizes a number of operating points in a matter of 0.7 s to explore possible control constraints.
Table 2 summarizes the comparison control performance between the proposed bilinear control and the standard PI control. Other important control aspects such as stability, region of attraction, parameter tuning, and robustness are described next, where the bilinear control presents great advantages compared with PI control.

4.2. Varying Total Converter Energy

Figure 4 aims to demonstrate the controllers’ management over the converter energy. In the second graphic, i v , d is shown, followed by i v , q , i c i r , 0 , W h , and W v .
For this test, the converter operating point was set to 70 % of the nominal power, at 0.05 s. The converter energy reference was changed by ± 10 % from the nominal value ( 3.6 MJ) in this test. The seventh graphic shows the total converter energy ( W h ), where either bilinear or PI controllers drive the energy to its reference value. The tracking process takes less than 50 ms for both controllers. The speed of convergence can be modified by modifying the controller gains according to the grid requirements, respecting the converter capability curve. The tracking process of total energy ( W h ) implies an overshoot in state variable i c i r , 0 (sixth graphic), because this state variable handles the DC output power, which is directly linked with stored energy. The PI behavior can also be regulated by gain tuning.

4.3. Disturbance in Energy Balance

One of the goals of the proposed control strategy is to keep the energy balance between the upper and lower arms of the MMC ( W v ). On the previous tests (four-quadrant operation), W v was kept around its reference value ( W v ¯ = 0 ). Although the MMC control strategies in general aim at keeping W v = 0 , it would be interesting to know how the operation with W v 0 would influence the proposed control strategy and the behavior of the considered states. With this purpose, a disturbance is considered as an imbalance between upper and lower converter arms. A positive value of W v means that there is more energy in the upper arms than in the lower, and a negative value means more energy in the lower arms. The eighth graphic in Figure 4 shows W v . Each state reacts in a distinct way, as follows:
  • The fourth and fifth graphics show circulating current i c i r , d and i c i r , q . Their behaviors are proportional to changes in W v , more pronounced in i c i r , d than in i c i r , q . Furthermore, to increase the upper arms’ energy (at t = 1 s), circulating currents i c i r , d q become negative (which means a change in the current direction), and to balance upper and lower energies, currents become positive (at t = 1.5 s). They also become positive to balance lower arms at t = 0.9 s. PI and bilinear have similar performance.
  • Under the bilinear controller, state variables i c i r , 0 and W h show overshoot for each change in W v , while under PI, these changes do not produce overshoot.
  • Lastly, the control response for an unbalanced condition of W v is shown in the eighth graphic. It shows that bilinear control reacts at around 100 ms to drive the state to a desired reference, while PI presents a faster response.
The test presented above pointed out for bilinear control that circulating current components, i c i r , d q 0 , and total energy, W h , actuate to manage the energy flow between the upper and lower converters’ arms.

4.4. Robustness Analysis

A set of scenarios are created where parameter uncertainties are considered. For each scenario, one parameter (R, L, or C S M ) assumes ± 10 % or ± 20 % from its nominal value. Meanwhile, the bilinear controller is computed using the parameters’ nominal values. At t = 0.01 s, the operating point is changed.
Figure 5 and Figure 6 consider uncertainties in arm resistance R and inductance L. There, one may see that the bilinear control still keeps stability with a very small difference from the nominal value, even with important parameter errors. Considering the variation in resistance, the controller was able to track the references, except for state W h . One explanation for this result is that reference calculation for i c i r c 0 and W h , Equations (31) and (35), is based on resistance value. For this reason, the controller is tracking a slightly different reference.
In the case of errors in L, one may observe a very small (≃1%) stationary error in i v d and W h ; meanwhile, the other states are perfectly driven to their references in all proposed scenarios.
For C uncertainties (Figure 7), the bilinear controller tracked the references for all states well. It should be highlighted that the capacitor voltage’s initial condition is the same for all the scenarios; therefore, the initial energy in each sub-module varies according to its capacitance value. That is the reason that the energy initial condition in Figure 7 is different for each value of C. However, it is interesting to note that the proposed controller steered W h to its reference for these different initial conditions. Similar behavior is observed for the scenario when the DC voltage is different than the nominal value; see Figure 8.
Figure 9 presents state variables facing errors in AC voltage. There is a time-varying response during transitory, with a fast convergence to the desired value. The states i v d and i v q presented minor stationary error. This can be explained by the calculation of their references ( i ¯ v d and i ¯ v q ), Equations (28) and (29), which takes into consideration the nominal values of the AC voltage.

4.5. SM Voltage Balancing

Figure 10 shows the voltage balancing of phase C lower arm sub-modules. Figure 10a shows that the voltage average value is almost constant during a 2-second simulation. Figure 10b shows in detail the voltage on each sub-module. It is noteworthy that the proposed controller regulates the total energy and energy balance between the lower and upper parts of the converter. The voltage balance among sub-modules is performed by the sorting algorithm after the PWM modulation stage. This sorting algorithm [7], however, is not a result of the present paper.
These simulations used a simulated model which is much more realistic than the one used for control design. Nevertheless, the controller was able to stabilize the system to its reference, as expected from the theoretical analysis. This illustrates the robustness of the proposed scheme.

5. Conclusions

This paper presents a bilinear control for an MMC converter, which directly controls all its state variables, AC, DC, and circulating currents and MMC energy, using a single control law. A rigorous mathematical proof is given for its stability, which is based on Lyapunov’s theory. The controller provides asymptotic stabilization for the three-phase MMC, and the use of a Lyapunov function implies formal verification of stability and an explicit region of attraction for the considered model. Moreover, a sensitivity analysis is shown, presenting very good robustness properties for realistic parameter errors. In addition, the small number of tuning parameters, and the existence of the Lyapunov function as a guideline, make its tuning much easier than for a set of cascaded PIs, which is the standard in industry.
A bilinear control technique is simulated using a detailed switching model in the Matlab Simscape Electrical environment. The theoretical and simulated results show that the proposed controller suits MMC time response and overshoot requirements, and provides an effective trade-off between power and energy control in the converter. Moreover, the energy management provided by the bilinear controller opens the way for new ancillary services for grid support, both for AC and DC grids, such as fast frequency response or synthetic inertia.
Comparing the bilinear control with respect to standard PI control, there is a rigorous stability proof for the first, with a well-established operating range, which is not true for PI. In the simulations, the controllers showed similar performance and robustness, with an improvement in the case of the bilinear control. Future work will use experimental setups to better investigate the properties of this control system. Likewise, optimization of the control gains will be investigated to further improve the performance of the controller.

Author Contributions

Conceptualization, all authors; methodology, G.C.d.O., G.D. and R.M.M.; validation, G.C.d.O. and F.P.; formal analysis, G.C.d.O. and G.D.; investigation, all authors; writing—original draft preparation, G.C.d.O.; writing—review and editing, G.D., M.J.C., F.P. and R.M.M.; supervision, G.D. and R.M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Erasmus Mundus SMART2 grant number 552042-EM-1-2014-1-FR-ERA MUNDUS-EMA2 coordinated by CentraleSupelec and by CAPES Foundation-Coordenação de Aperfeiçoamento de Pessoal de Nível Superior-Brasil.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Farsi, M.; Liu, J. Nonlinear Optimal Feedback Control and Stability Analysis of Solar Photovoltaic Systems. IEEE Trans. Control Syst. Technol. 2020, 28, 2104–2119. [Google Scholar] [CrossRef]
  2. Schmuck, C.; Woittennek, F.; Gensior, A.; Rudolph, J. Feed-Forward Control of an HVDC Power Transmission Network. IEEE Trans. Control Syst. Technol. 2014, 22, 597–606. [Google Scholar] [CrossRef]
  3. Arrillaga, J.; Liu, Y.; Watson, N. Flexible Power Transmission: The HVDC Options; Wiley: Hoboken, NJ, USA, 2007. [Google Scholar]
  4. Lesnicar, A.; Marquardt, R. An innovative modular multilevel converter topology suitable for a wide power range. In Proceedings of the 2003 IEEE Bologna Power Tech Conference Proceedings, Bologna, Italy, 23–26 June 2003; Volume 3, p. 6. [Google Scholar] [CrossRef]
  5. Jankovic, M.; Costabeber, A.; Watson, A.; Clare, J.C. Arm-Balancing Control and Experimental Validation of a Grid-Connected MMC With Pulsed DC Load. IEEE Trans. Ind. Electron. 2017, 64, 9180–9190. [Google Scholar] [CrossRef]
  6. Quester, M.; Loku, F.; El Azzati, O.; Noris, L.; Yang, Y.; Moser, A. Investigating the Converter-Driven Stability of an Offshore HVDC System. Energies 2021, 14, 2341. [Google Scholar] [CrossRef]
  7. Jovcic, D.; Ahmed, K. High Voltage Direct Current Transmission: Converters, Systems and DC Grids; Wiley: Hoboken, NJ, USA, 2015. [Google Scholar]
  8. Harnefors, L.; Antonopoulos, A.; Norrga, S.; Angquist, L.; Nee, H.P. Dynamic Analysis of Modular Multilevel Converters. IEEE Trans. Ind. Electron. 2013, 60, 2526–2537. [Google Scholar] [CrossRef]
  9. Zhang, H.; Belhaouane, M.M.; Colas, F.; Kadri, R.; Gruson, F.; Guillaud, X. On comprehensive description and analysis of MMC control design: Simulation and experimental study. IEEE Trans. Power Deliv. 2020, 36, 244–253. [Google Scholar] [CrossRef]
  10. de Oliveira, G.C.; Damm, G.; Monaro, R.M.; Lourenço, L.F.; Carrizosa, M.J.; Lamnabhi-Lagarrigue, F. Nonlinear control for modular multilevel converters with enhanced stability region and arbitrary closed loop dynamics. Int. J. Electr. Power Energy Syst. 2021, 126, 106590. [Google Scholar] [CrossRef]
  11. Carrizosa, M.J.J.; Bergna, G.; Arzande, A.; Damm, G.; Alou, P.; Benchaib, A.; Lamnabhi-Lagarrigue, F. Stability of DC/DC three terminals converter using Modular Multilevel Converters for HVDC systems. In Proceedings of the 17th European Conference on Power Electronics and Applications (EPE-ECCE Europe 2015), Genève, Switzerland, 8–10 September 2015. [Google Scholar] [CrossRef]
  12. Bergna-Diaz, G.; Freytes, J.; Guillaud, X.; D’Arco, S.; Suul, J.A. Generalized Voltage-Based State-Space Modeling of Modular Multilevel Converters With Constant Equilibrium in Steady State. IEEE J. Emerg. Sel. Top. Power Electron. 2018, 6, 707–725. [Google Scholar] [CrossRef]
  13. Bergna-Diaz, G.; Suul, J.A.; Berne, E.; Vannier, J.C.; Molinas, M. Optimal Shaping of the MMC Circulating Currents for Preventing AC-Side Power Oscillations From Propagating Into HVdc Grids. IEEE J. Emerg. Sel. Top. Power Electron. 2019, 7, 1015–1030. [Google Scholar] [CrossRef]
  14. Zhang, Z.; Jin, Y.; Xu, Z. Modeling and Control of Modular Multilevel Matrix Converter for Low-Frequency AC Transmission. Energies 2023, 16, 3474. [Google Scholar] [CrossRef]
  15. Bergna, G.; Vannier, J.; Lefranc, P.; Arzande, A.; Berne, E.; Egrot, P.; Molinas, M. Modular multilevel converter leg-energy controller in rotating reference frame for voltage oscillations reduction. In Proceedings of the 2012 3rd IEEE International Symposium on Power Electronics for Distributed Generation Systems (PEDG), IEEE, Aalborg, Denmark, 25–28 June 2012; pp. 698–703. [Google Scholar]
  16. Saad, H.; Guillaud, X.; Mahseredjian, J.; Dennetiere, S.; Nguefeu, S. MMC capacitor voltage decoupling and balancing controls. IEEE Trans. Power Deliv. 2014, 30, 704–712. [Google Scholar] [CrossRef]
  17. Münch, P.; Görges, D.; Izák, M.; Liu, S. Integrated current control, energy control and energy balancing of Modular Multilevel Converters. In Proceedings of the IECON 2010—36th Annual Conference on IEEE Industrial Electronics Society, Glendale, AZ, USA, 7–10 November 2010; pp. 150–155. [Google Scholar] [CrossRef]
  18. Vatani, M.; Hovd, M.; Saeedifard, M. Control of the Modular Multilevel Converter Based on a Discrete-Time Bilinear Model Using the Sum of Squares Decomposition Method. IEEE Trans. Power Deliv. 2015, 30, 2179–2188. [Google Scholar] [CrossRef]
  19. Shinoda, K.; Benchaib, A.; Dai, J.; Guillaud, X. Energy control of modular multilevel converter with a novel analytic filter. In Proceedings of the 2016 18th European Conference on Power Electronics and Applications (EPE’16 ECCE Europe), Karlsruhe, Germany, 5–9 September 2016; pp. 1–10. [Google Scholar] [CrossRef]
  20. Ma, Y.; Fan, L. Circulating current and DC current ripple control in MMC under unbalanced grid voltage. In Proceedings of the 2015 North American Power Symposium (NAPS), Charlotte, NC, USA, 4–6 October 2015; pp. 1–6. [Google Scholar] [CrossRef]
  21. Saad, H.; Dennetière, S.; Mahseredjian, J.; Delarue, P.; Guillaud, X.; Peralta, J.; Nguefeu, S. Modular Multilevel Converter Models for Electromagnetic Transients. IEEE Trans. Power Deliv. 2014, 29, 1481–1489. [Google Scholar] [CrossRef]
  22. Antonopoulos, A.; Angquist, L.; Nee, H. On dynamics and voltage control of the Modular Multilevel Converter. In Proceedings of the 2009 13th European Conference on Power Electronics and Applications, Barcelona, Spain, 8–10 September 2009; pp. 1–10. [Google Scholar]
  23. Iovine, A.; Siad, S.B.; Damm, G.; De Santis, E.; Di Benedetto, M.D. Nonlinear Control of a DC MicroGrid for the Integration of Photovoltaic Panels. IEEE Trans. Autom. Sci. Eng. 2017, 14, 524–535. [Google Scholar] [CrossRef]
  24. Bacha, S.; Munteanu, I.; Bratcu, A. Power Electronic Converters Modeling and Control: With Case Studies; Advanced Textbooks in Control and Signal Processing; Springer: London, UK, 2013. [Google Scholar]
  25. Mohler, R.R. Bilinear Control Processes: With Applications to Engineering, Ecology, and Medicine; Mathematics in Science and Engineering; Elsevier Science: Amsterdam, The Netherlands, 1974. [Google Scholar]
  26. Pardalos, P.M.; Yatsenko, V. Optimization and Control of Bilinear Systems; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  27. Carrizosa, M.J.; Damm, G.; Benchaib, A.; Alou, P.; Netto, M.; Lamnabhi-Lagarrigue, F. Bilinear and nonlinear control algorithms for a DC/DC converter for Multi-Terminal HVDC networks. IFAC Proc. Vol. 2014, 47, 523–528. [Google Scholar] [CrossRef]
  28. Pardalos, P.; Yatsenko, V. Optimization and Control of Bilinear Systems: Theory, Algorithms, and Applications; Springer Optimization and Its Applications; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  29. Landau, I. On the Optimal Regulator Problem and Stabilization of Bilinear Systems; Laboratoire d’Automatique (CNRS): Grenoble, France, 1979. [Google Scholar]
  30. Panteley, E.; Loría, A. Synchronization and Dynamic Consensus of Heterogeneous Networked Systems. IEEE Trans. Autom. Control 2017, 62, 3758–3773. [Google Scholar] [CrossRef]
  31. Yazdani, A.; Iravani, R. Voltage-Sourced Converters in Power Systems; Wiley Online Library: Hoboken, NJ, USA, 2010; Volume 34. [Google Scholar]
Figure 1. Three-phase MMC.
Figure 1. Three-phase MMC.
Energies 16 06713 g001
Figure 2. PI and bilinear control diagram.
Figure 2. PI and bilinear control diagram.
Energies 16 06713 g002
Figure 3. Bilinear and PI controllers in four-quadrant operation.
Figure 3. Bilinear and PI controllers in four-quadrant operation.
Energies 16 06713 g003
Figure 4. Bilinear and PI controllers with energy step changes.
Figure 4. Bilinear and PI controllers with energy step changes.
Energies 16 06713 g004
Figure 5. The system state over variation of parameter R with bilinear controller.
Figure 5. The system state over variation of parameter R with bilinear controller.
Energies 16 06713 g005
Figure 6. States over variation of parameter L with bilinear controller.
Figure 6. States over variation of parameter L with bilinear controller.
Energies 16 06713 g006
Figure 7. States over variation of parameter C with bilinear controller.
Figure 7. States over variation of parameter C with bilinear controller.
Energies 16 06713 g007
Figure 8. Fluctuation of DC voltage with bilinear controller.
Figure 8. Fluctuation of DC voltage with bilinear controller.
Energies 16 06713 g008
Figure 9. Error on AC voltage with bilinear controller.
Figure 9. Error on AC voltage with bilinear controller.
Energies 16 06713 g009
Figure 10. SMs’ voltage balancing for lower arm of phase C.
Figure 10. SMs’ voltage balancing for lower arm of phase C.
Energies 16 06713 g010
Table 1. Parameters of simulated system and bilinear control’s gains.
Table 1. Parameters of simulated system and bilinear control’s gains.
ParameterValueParameterValue
S M M C 50 MVA α 1 0.5
V A C 30 kV α 2 0.5
V D C 180 kV α 3 0.5
L c 5 mH α 4 0.5
L 14 mH α 5 0.5
R c 0.03 Ω Γ 1 1
C S M 3 mF Γ 2 1
Freq 60 Hz Φ c 1
R 0.5 Ω N20
Table 2. Performance comparison between bilinear and PI control.
Table 2. Performance comparison between bilinear and PI control.
PIBilinearComments
i v d q Larger overshoots
and oscillations
Smaller overshootsSimilar settling
time for both
i c i r c d q Good reference tracking
(nonzero reference)
Steady state errors,
but closer to zero
Similar oscillation
i c i r c 0 Very good track
of reference
Larger oscillationsDifferent references
for each controller
W h Very good track
of reference
Steady state errorsBilinear with
oscillating dynamics
W v Larger overshootsSmaller overshootsBoth present steady
state errors
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Costa de Oliveira, G.; Machado Monaro, R.; Damm, G.; Perez, F.; Jimenez Carrizosa, M. Bilinear Quadratic Feedback Control of Modular Multilevel Converters. Energies 2023, 16, 6713. https://doi.org/10.3390/en16186713

AMA Style

Costa de Oliveira G, Machado Monaro R, Damm G, Perez F, Jimenez Carrizosa M. Bilinear Quadratic Feedback Control of Modular Multilevel Converters. Energies. 2023; 16(18):6713. https://doi.org/10.3390/en16186713

Chicago/Turabian Style

Costa de Oliveira, Guacira, Renato Machado Monaro, Gilney Damm, Filipe Perez, and Miguel Jimenez Carrizosa. 2023. "Bilinear Quadratic Feedback Control of Modular Multilevel Converters" Energies 16, no. 18: 6713. https://doi.org/10.3390/en16186713

APA Style

Costa de Oliveira, G., Machado Monaro, R., Damm, G., Perez, F., & Jimenez Carrizosa, M. (2023). Bilinear Quadratic Feedback Control of Modular Multilevel Converters. Energies, 16(18), 6713. https://doi.org/10.3390/en16186713

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop