Advanced
Finite Control Set Model Predictive Control of AC/DC Matrix Converter for Grid-Connected Battery Energy Storage Application
Finite Control Set Model Predictive Control of AC/DC Matrix Converter for Grid-Connected Battery Energy Storage Application
Journal of Power Electronics. 2015. Jul, 15(4): 1006-1017
Copyright © 2015, The Korean Institute Of Power Electronics
  • Received : October 12, 2014
  • Accepted : March 06, 2015
  • Published : July 20, 2015
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Bo Feng
State Key Laboratory of Advanced Electromagnetic Engineering and Technology, School of Electrical and Electronic Engineering, Huazhong University of Science and Technology, Wuhan, China
fengbohust@163.com
Hua Lin
State Key Laboratory of Advanced Electromagnetic Engineering and Technology, School of Electrical and Electronic Engineering, Huazhong University of Science and Technology, Wuhan, China

Abstract
This paper presents a finite control set model predictive control (FCS-MPC) strategy for the AC/DC matrix converter used in grid-connected battery energy storage system (BESS). First, to control the grid current properly, the DC current is also included in the cost function because of input and output direct coupling. The DC current reference is generated based on the dynamic relationship of the two currents, so the grid current gains improved transient state performance. Furthermore, the steady state error is reduced by adding a closed-loop. Second, a Luenberger observer is adopted to detect the AC input voltage instead of sensors, so the cost is reduced and the reliability can be enhanced. Third, a switching state pre-selection method that only needs to evaluate half of the active switching states is presented, with the advantages of shorter calculation time, no high dv / dt at the DC terminal, and less switching loss. The robustness under grid voltage distortion and parameter sensibility are discussed as well. Simulation and experimental results confirm the good performance of the proposed scheme for battery charging and discharging control.
Keywords
I. INTRODUCTION
With the spring up of renewable energy generation and smart grid, the energy storage technology, especially the battery energy storage system (BESS), has found increased applications. Voltage source converter (VSC) is widely used as an interface between grid and battery. However, the battery usually has a lower voltage than the grid. Therefore, a large number of battery cells must be connected in series or an additional bi-directional DC/DC converter is required to achieve voltage level matching [1] , [2] .
AC/DC matrix converter is a novel converter topology derived from the well-known AC/AC matrix converter [3] , [4] . It has been studied in several applications, such as DC motor drive [5] , power supply [6] , and battery charger [7] , [8] . It inherits many advantages from the AC/AC matrix converter, such as: 1) compact structure without bulky electrolytic capacitor, 2) sinusoidal AC current and controllable power factor, and 3) four-quadrant operation capability [9] - [11] . More importantly, it is intrinsically a buck converter and features a tight DC voltage regulation in a wide range. Thus, it is an attractive alternative to the VSC-based topologies.
In grid-connected BESS application, grid current is the primary control objective, while the reported studies on AC/DC matrix converter mainly focused on DC side control. Traditionally, the converter is controlled with the aid of a modulator. In [7] and [8] , space vector modulation (SVM) is utilized. A dual-loop control scheme with outer DC voltage and inner DC current loops is used. The AC input current is set to be in phase with the grid voltage, while its amplitude is indirectly controlled. In [12] , the DC voltage controller sets the d -axis AC input current reference and then the AC input current is directly controlled using a one-cycle control strategy without a modulator. In these strategies, the AC side performance greatly depends on the DC side linear controllers. Moreover, the AC input current is controlled rather than the grid current. Because of the presence of AC side filter capacitors, the amplitude and phase of the grid current cannot be controlled precisely.
Finite control set model predictive control (FCS-MPC, hereafter referred to as MPC) provides a way to regulate the grid current directly and quickly. It takes advantage of the discrete and nonlinear nature of switching power converters and proves to be a simple and effective alternative to traditional control schemes. It offers many benefits, such as fast dynamic response, easy inclusion of nonlinearities and constraints, and does not need a modulator [13] , [14] . The MPC scheme uses a discrete model to predict the system behavior in every possible switching state. Then, the predicted values are used to evaluate a cost function that can be customized flexibly according to specific requirements. Finally, the switching state that leads to the minimum cost function is selected and applied in the next sampling instant.
The MPC strategies for AC/AC matrix converter and current source rectifier can be used as references for the control of the AC/DC matrix converter because these converters have similar structures. Among them, the MPC with imposed sinusoidal source current shows excellent currents at both output and input, which is very attractive [15] - [17] . In this method, besides the load current control, the source current is also forced to follow a sinusoidal waveform. The source current reference is obtained from the load current reference according to the steady state input and output power balance. However, the difference between the dynamic responses of the two currents is not considered. In addition, the generated reference relates to the converter efficiency that varies with the operation condition. The accuracy of this parameter will affect the steady state tracking error, as analyzed in a study on predictive controlled back-to-back converter [18] .
MPC requires the information of both AC source and input voltages [15] - [17] , whereas the conventional control schemes only need either of the voltages. Voltage sensors increase the system size and cost and raises the risk of control failure caused by sensor malfunction. For the voltage source inverter, an observer can be used to estimate the load current without increasing the number of sensors [19] . However, using an observer to estimate the AC input voltage in an AC/DC matrix converter has not been mentioned in literature.
Evaluating all the switching states will lead to several problems. First, the calculation process is time-consuming, which limits the sampling frequency and eventually degrades the control performance. Many methods have been proposed to simplify the MPC algorithm for two- and multi-level VSCs [20] - [22] . Second, the output side will bear too large dv / dt [23] . Third, the switching loss may be higher. Therefore, discussing the switching state selection method for the MPC of the AC/DC matrix converter is worthwhile.
The contribution of this paper is the MPC implementation on an AC/DC matrix converter for grid-connected BESS application. To address the above issues, three aspects are investigated. First, the grid and DC currents are both controlled using a single cost function. The reference of the DC current is generated considering the transient state relationship of the two currents, and a proportional-integral (PI) controller is added to compensate for the steady state error. Second, the methods to estimate the AC input voltage are discussed. A Luenberger observer, which reduces cost and enhances reliability, is adopted. Third, to reduce the calculation amount and switching loss, as well as eliminate high dv / dt , a switching states pre-selection method based on the input current sector is proposed. Moreover, the robustness under grid voltage distortion and parameter sensitivity are discussed through simulation. Experiments are carried out on an AC/DC matrix converter prototype with a battery load to evaluate the proposed methods.
II. PREDICTIVE MODEL
The AC/DC matrix converter circuit is shown in Fig. 1 . The power stage consists of six bidirectional switches, and each switch is implemented using two insulated gate bipolar transistors (IGBTs) with freewheeling diodes in common emitter configuration. The AC side LC filter is generally required for matrix converters to prevent over voltage and attenuate high-frequency input current harmonics. Subscripts “ s ,” “ i ,” and “ dc ” denote grid side, AC input side, and DC side quantities, respectively. Rf , Lf , and Cf are the AC filter parameters. Ro , Lo , and Co are the DC filter parameters. Sij ( i = P , N and j = a , b , c ) are the switching signals, Sij = 1, 0. Complying with the restriction that no short-circuit at the AC side and no open-circuit at the DC side, nine switching states are feasible in the AC/DC matrix converter.
PPT Slide
Lager Image
AC/DC matrix converter circuit.
The AC input currents and the DC voltage are derived as
PPT Slide
Lager Image
PPT Slide
Lager Image
To reduce the amount of calculation, the MPC algorithm is implemented in the stationary αβ frame. At the AC side, the second order filter dynamics can be described as
PPT Slide
Lager Image
PPT Slide
Lager Image
where the bold letters represent the matrix of α - and β - components.
The state-space representation of the continuous-time system is
PPT Slide
Lager Image
with x = [ is ui ] T and u = [ ii us ] T .
Considering the sampling period Ts , the discrete-time model of the AC side filter can be obtained by
PPT Slide
Lager Image
where
PPT Slide
Lager Image
According to Equ. (6), the grid current is at time instant k + 1 is expressed as
PPT Slide
Lager Image
where Λ (m, n) is the (m, n) element of matrix Λ .
To avoid the performance deterioration caused by the delay effect in digital implementation, a two-step-ahead prediction method [13] is adopted. By shifting Equ. (7) one step forward, the grid current at time instant k + 2 is obtained by Equ. (8). For simplicity, us ( k + 1) ≈ us ( k ) because of its small change in one sampling period.
PPT Slide
Lager Image
The DC side behavior can be described by Equ. (9). Because only high-frequency currents flow through the DC filter capacitor, this capacitor is ignored in the modeling.
PPT Slide
Lager Image
Because Equ. (9) is a first-order system, the DC current can be predicted simply using the forward Euler approximation as Equs. (10) and (11), where uB ( k + 1) ≈ uB ( k ).
PPT Slide
Lager Image
PPT Slide
Lager Image
Fig. 2 provides an overview of the proposed MPC strategy. Details of the cost function design, AC input voltage estimation, and switching states pre-selection based on sector judgment are discussed in the following sections.
PPT Slide
Lager Image
Control scheme of the proposed MPC strategy.
III. COST FUNCTION DESIGN
- A. Cost Function
The goal of a grid-connected BESS is to control the grid current according to the grid command. However, in an AC/DC matrix converter, the AC and DC sides are directly connected through bidirectional switches without energy storage capability, so they should be controlled together. Fig. 3 presents the simulation waveforms of the predictive control with a cost function involving only the grid current. Though the grid current follows a 5A reference, the DC current is in reverse and far beyond the limit.
PPT Slide
Lager Image
Simulation waveforms with only grid current control.
Thanks to the flexibility of controlling multiple variables in MPC, the grid and DC currents can be controlled simultaneously by including them in a single cost function as follows:
PPT Slide
Lager Image
where Is is the grid current amplitude (negative Is corresponds to the discharging mode) and λ is the weighting factor, which assigns more priority to one variable or another.
Usually, λ is determined using repetitive simulations to achieve satisfactory results for the control objective. The total harmonic distortion (THD) and the amplitude error Is_err of the grid current are used to evaluate λ because this paper aims at grid side control. Fig. 4 shows the simulation results of the two indices in terms of λ , under a –5A grid current. Obviously, a smaller λ leads to a better grid current because the grid current is assigned more importance. However, when λ is less than 0.8, the grid current quality deteriorated quickly. This is because the DC current control is greatly weakened, which in turn affects the grid current. A too small λ will even cause instability. Considering the grid current performance and system stability, λ is set to 1.
PPT Slide
Lager Image
Simulation results of grid current THD and amplitude error Is_err in terms of weighting factor λ.
According to the AC and DC side power balance, the ratio of grid current amplitude to DC current is equal to the ratio of battery voltage to grid voltage amplitude Us . Then, the cost function can be simplified into Equ. (13). The coefficient of the DC term is approximately constant because the variations on the grid and battery voltages are very small.
PPT Slide
Lager Image
- B. DC Current Reference Generation
A problem that arises when applying the cost function is the unavoidable coupling of the grid and DC currents. To control the grid current to track the grid command properly, the DC current reference cannot be arbitrarily chosen. Therefore, the generation of this reference is a key issue.
1) Method 1: The MPC with imposed sinusoidal source current method offers an effective method to calculate the AC source current reference from the desired AC or DC load current based on active power balance. This idea is also applicable in this paper, except that the grid current becomes the primary target and the DC current the subordinate control objective.
The accurate AC input side active power model is
PPT Slide
Lager Image
At the DC output side, the active power is expressed as
PPT Slide
Lager Image
Considering the converter loss in a real system, the relationship between the input and output active power is
PPT Slide
Lager Image
where η is the efficiency in charging mode or the reciprocal of the efficiency in discharging mode.
From Equs. (14)–(16), the DC current reference idc * is determined by the amplitude of grid current reference is * as
PPT Slide
Lager Image
with A = 1−8 ωs 2 Lfcf .
Equ. (17) only considers the steady state relationship of the two currents and ignores their dynamic behaviors, thus a step change in is * will immediately generate a step change in idc * . Using MPC, fast regulation of DC current is possible. However, the difference between the AC and DC side filter parameters indicates that the DC current essentially has a much slower dynamic response than the grid current. A sudden change in the DC current may incur a large overshoot or undershoot in the grid current. Therefore, forcing both currents to track simultaneous step-changed references is inappropriate. In fact, the change rate of idc * should be limited.
Moreover, the obtained reference relies on accurate system parameters, especially on the converter efficiency. η varies with the operation condition (power level, charging or discharging, and so on) and its value should be adjusted empirically, which causes inconvenience in practical implementation. The mismatch of η will lead to inaccurate idc * and finally affect the steady state tracking error of the grid current.
2) Method 2: The PI control enables the perfect tracking of constant reference, so the DC current reference can also be obtained using a PI control loop. Then, the steady state error caused by inaccurate efficiency parameter is eliminated.
The DC current reference can be expressed as
PPT Slide
Lager Image
where kp and ki are the PI parameters.
Compared with Method 1, the dynamics of idc * is greatly reduced because idc * is completely determined by the error of grid current amplitude. Besides, tuning the PI parameters is a difficult task.
3) Method 3: To address the transient and steady state issues in Method 1 well, a novel strategy is proposed.
The AC filter has much faster dynamics than the DC filter, so the AC side dynamics can be neglected. Then, the following equation is established:
PPT Slide
Lager Image
Applying a small-signal perturbation to Equ. (19), the transfer function from the grid current to the DC current is obtained by:
PPT Slide
Lager Image
where Idc is the DC current at the operating point.
Equ. (20) indicates the transient state relationship of the two currents, that is, the change in idc lags behind the change in Is . This should be considered when generating idc * .
Moreover, idc * is still derived in an open-loop manner like Method 1, thus also faces the steady state error problem. In view of this, a PI closed-loop is added as a supplement.
The DC current reference is finally expressed as
PPT Slide
Lager Image
where τ = LoIdc / uB .
Fig. 5 shows the DC current reference generation Methods 1–3. In Method 3, idc * consists of two parts: the inertia term directly obtained from Is * and the PI term related to the actual error Is * Is . The inertia term, which is also a feed-forward of the grid current, plays the dominant role in both transient and steady states. By introducing 1/( τs + 1), the change rate of idc * is limited, thus the overshoot of the grid current in transient state is reduced compared with that in Method 1. The PI parameters can be much smaller than those in Method 2 because the PI term is only utlized to deal with the tracking error in the steady state.
PPT Slide
Lager Image
DC current reference generation. (a) Method 1. (b) Method 2. (c) Method 3.
IV. AC INPUT VOLTAGE ESTIMATION
According to Equ. (8), among all the variables involved in predicting is ( k + 2), us ( k ) is needed for synchronization and is ( k ) is the control objective, thus the two variables should be measured. ii ( k ) cannot be measured but can be calculated using the optimal switching state calculated in the last time instant, and ii ( k + 1) relates to the optimal switching state for the next time instant. The input voltage can be measured. However, ui ( k ) and ui ( k + 1) are estimated to prevent increasing the number of sensors.
A simple way to estimate the AC input voltage is using the derivative Equ. (3). The AC input voltage at time instant k can be calculated as
PPT Slide
Lager Image
where “^” indicates an estimated value.
Then, the future AC input voltage can be derived using a simple extrapolation of the past estimated values as follow.
PPT Slide
Lager Image
However, this approach operates in an open-loop manner and is very sensitive to measurement noises because of the use of grid current derivation. Therefore, the method using a closed-loop observer is preferred.
The Luenberger observer uses the system model to predict the system state and then feeds back the difference between the measured and estimated output and corrects the model with this error signal to minimize the divergence between the model and the actual system. It can correct the model inaccuracy, eliminate the disturbance effect and noise, and has high reliability. Moreover, it presents a predictive nature so that a measurement at time instant k results in an observation of the state variables valid in time instant k + 1, and this is exactly what the two-step-ahead prediction needs.
Suppose the system output is the grid current,
PPT Slide
Lager Image
The full-order observer is used to estimate the state vector x . It consists of an open-loop model of the system and a correcting term related to the measured output. The observer dynamics can be represented by
PPT Slide
Lager Image
where
PPT Slide
Lager Image
and G = [ h 1 h 2 ] T is the gain matrix.
The error vector of the observer is defined as
PPT Slide
Lager Image
thus its dynamics can be expressed as
PPT Slide
Lager Image
To make the observer error converge to zero with time for any initial value, the observer gain matrix should be selected such that A-GC is asymptotically stable. In other words, all the eigenvalues of A-GC should have negative real parts. The characteristic equation of the system is
PPT Slide
Lager Image
The gain matrix can be determined by applying the pole placement method. Suppose the observer has a pair of conjugate complex poles a ± jb , then the desired characteristic equation can be expressed as
PPT Slide
Lager Image
Equs. (27) and (28) should be identical, and the elements of G can be obtained by matching the coefficients of each power of s as follows:
PPT Slide
Lager Image
A-GC determines the observer dynamics. The gain matrix G should be designed such that the observer poles are farther from the imaginary axis than the original system. If the observer poles are closer to the imaginary axis, obtaining the actual value of the state variable will take a long time, whereas too far away from the imaginary axis, the system will be very sensitive to the noises and might be unstable. Therefore, a compromise should be made between bandwidth and noise rejection.
The discrete form of the observer is
PPT Slide
Lager Image
where
PPT Slide
Lager Image
and
PPT Slide
Lager Image
Then, the AC input voltage at time k + 1 can be calculated using the present measured grid voltages and currents.
V. ACTIVE SWITCHING STATE SELECTION METHOD
The current prediction and cost function evaluation for all feasible switching states need a large amount of calculation, especially when multiple constrains are considered. If the calculations take too much time, a long sampling period should be used, which will degrade the control performance.
On the other hand, negative DC terminal voltage is allowed in an AC/DC matrix converter. For instance, in the charging mode, if the actual DC current is larger than the reference, the switching state corresponding to the negative DC terminal voltage may be selected because it can decrease the current quickly. However, jumping between positive and negative voltages will lead to large DC side dv / dt and high switching loss. The increase in switching loss results from the higher switching frequency because two commutations take place in one sampling period and higher switched voltage that is approximately proportional to the switching loss [24] , [25] .
Evaluating a subset instead of all the switching states will help shorten the calculation time. Furthermore, if the subset is selected properly, the negative DC terminal voltage can be avoided.
The input current sector concept in SVM is introduced. Take sector 1 as an example, where the sector angle θ lies in (–π/6, π/6) [7] . Fig. 6 (a) shows the AC input currents and the corresponding line voltages in charging mode. Sector 1 has three positive line voltages that can be applied at the DC terminal, namely, uiab , uiac, and uicb in the first half of the sector (–π/6, 0) and uiab , uiac, and uibc in the second half of the sector (0, π/6). Thus, the optimization process needs to evaluate four active switching states ab , ac , cb , and bc , and the negative voltage still exists if bc is selected in the first half of the sector or cb is selected in the second half of the sector.
PPT Slide
Lager Image
AC input current sectors. (a) AC input currents and line voltages. (b) Proposed sector definition.
If the π/3 sector interval is shifted to the right by π/6, then uiab , uiac, and uibc are all positive in the newly defined sector 1 (0, π/3). As a result, the number of candidate active switching states is reduced to only three and the negative DC terminal voltage is eliminated. The new sector definition is shown in Fig. 6 (b).
Ideally, the evaluation time of the active switching states can be reduced by half. However, judging the input current sector introduces extra computation; therefore, this procedure should also be optimized.
The actual AC input currents are PWM waveforms, so the fundamental component should be derived first.
For sinusoidal and balanced three-phase voltages, the following equations can be deduced [26]
PPT Slide
Lager Image
Substituting Equ. (31) into Equ. (4), the fundamental AC input currents i and i are derived by Equ. (32). i and i can be readily calculated using the measured grid current and the estimated input voltage.
PPT Slide
Lager Image
The conventional sector judgment method requires complex trigonometric and modulus manipulations [27] . Because the value of sector angle is not necessary for MPC, a simple approach without sector angle calculation is adequate.
Fig. 6 (a) shows that in sector 1 (0, π/3), the input phase currents have the following relationship.
PPT Slide
Lager Image
By replacing iia , iib , and iic with i and i , Equ. (33) is rewritten as
PPT Slide
Lager Image
Similar relationships can be seen in the other sectors. P corresponds to the sector number and is defined as
PPT Slide
Lager Image
where
PPT Slide
Lager Image
The sector number calculation and the candidate active switching states are listed in Table I , where P 0 = sign( i ),
PPT Slide
Lager Image
SECTOR JUDGEMENT AND ACTIVE SWITCHING STATE SUBSET
PPT Slide
Lager Image
SECTOR JUDGEMENT AND ACTIVE SWITCHING STATE SUBSET
The sector number can be easily calculated using Equs. (32) and (35), and only a few addition and multiply operations, as well as a simple lookup table are used, which reduces the computational burden on the microprocessor.
VI. SIMULATION AND EXPERIMENTAL RESULTS
The proposed control methods are evaluated through Matlab simulation and experiment. Experimentally, the power stage is constructed with 12 discrete IGBTs of type IKW40N120. The MPC algorithm is executed in a DSP of type TMS320F28335, and the four-step current commutation is implemented using a FPGA of type EP2C8. The DC load is a 10 × 12 V/100 Ah lead-acid battery pack. Other system parameters are as follows: grid line voltage 200 Vrms/50 Hz; AC filter Lf = 1.2 mH, Cf = 10 μF, Rf = 0.1 Ω; DC filter Lo = 10 mH, Co = 20 μF, Ro = 0.1 Ω; sampling period Ts = 20 μs.
- A. Comparison of AC Input Voltage Estimation Methods
Fig. 7 shows the experimental waveforms using the simple estimation method based on derivative equation and with the Luenberger observer. The grid current is set to 3A charging. As observed, the derivative-based method results in a noisy input voltage that contains many harmonics, whereas the observer-based method obtains a clear input voltage with low THD. Because of the low pass filter effect of the observer, the estimated voltage presents less high frequency harmonics than the measured value, as shown in Fig. 7 (b). The grid current waveforms indicate that the derivative-based method leads to poor control performance, whereas the observer-based method is suitable for MPC control.
PPT Slide
Lager Image
Experimental waveforms using different AC input voltage estimation methods. (a) Derivative-based method. (b) Observer-based method. u [160 V/div], isa [6 A/div], time [10 ms/div].
- B. Comparison of Active Switching State Selection Methods
Fig. 8 shows the experimental waveforms using the conventional six active switching states method and with the proposed three active switching states method. For fair comparison, the DC current references are all generated using Equ. (17) and η is empirically adjusted to 0.94 in charging mode and 1.11 in discharging mode to obtain a ±5A grid current. The figure shows that with the conventional method, many negative voltage pulses are generated at the DC side and the converter suffers from large dv/dt , especially in the discharging mode. By contrast, selecting among three active switching states strongly reduces the negative voltage pulses at the DC side and only slightly affects the current quality.
PPT Slide
Lager Image
Experimental waveforms using different active switching state selection methods. (a) Conventional method, charging. (b) Proposed method, charging. (c) Conventional method, discharging. (d) Proposed method, discharging. uPN [200 V/div], usa [160 V/div], isa [10 A/div], idc [10 A/div], time [5 ms/div].
Table II lists the experimental data of –5A discharging. Compared with the conventional method, the proposed method reduces the evaluation time of the cost function by 35%. Moreover, the average switching frequency and the average switched voltage are reduced by 8% and 14%, respectively, which means a smaller switching loss.
EXPERIMENTAL DATA WITH DIFFERENT ACTIVE SWITCHING STATE SELECTION METHODS
PPT Slide
Lager Image
EXPERIMENTAL DATA WITH DIFFERENT ACTIVE SWITCHING STATE SELECTION METHODS
In Method 1, η is set to 1 in both the charging and discharging modes. Like Is * , the generated idc * changes in a stepped manner. The grid current presents a very large overshoot in the charging mode and undershoot in the discharging mode because of the drastic change in the DC current. On the other hand, because the actual efficiency is less than one, idc * is not the accurate value that results in the desired grid current. Therefore, the grid current has a steady state amplitude error, and this problem is even worse in the discharging mode because of lower efficiency.
In Method 2, kp = 0.4 and ki = 800. As expected, the grid current amplitude error is eliminated. However, idc * shows very different dynamic behaviors in the charging and discharging modes. idc * takes a long time to reach the steady state in the charging mode. On the contrary, it responses fast and even oscillates in the discharging mode. Therefore, designing PI parameters suitable for both charging and discharging modes is difficult. The PI control is very sensitive to input noise and ripple, thus resulting in a distorted idc * . The DC current is forced to follow this non-ideal reference, so the grid current quality is damaged. As shown in Table III , the grid current has a higher THD.
EXPERIMENTAL DATA WITH DIFFERENT DC CURRENT REFERENCE GENERATION METHODS
PPT Slide
Lager Image
EXPERIMENTAL DATA WITH DIFFERENT DC CURRENT REFERENCE GENERATION METHODS
- C. Comparison of DC Current Reference Generation Methods
Figs. 9 and 10 show the experimental waveforms with different DC current reference generation methods. Table III lists the grid current amplitude and THD. The grid current changes from 3A to 5A and –3A to –5A, respectively.
PPT Slide
Lager Image
Experimental waveforms using different DC current reference generation methods, charging. (a) Method 1. (b) Method 2. (c) Method 3. isa* and isa [3 A/div], idc* and idc [5 A/div], time [10 ms/div].
PPT Slide
Lager Image
Experimental waveforms using different DC current reference generation methods, discharging. (a) Method 1. (b) Method 2. (c) Method 3. isa* and isa [3 A/div], idc* and idc [5 A/div], time [10 ms/div].
In Method 3, because the inertia term considered the dynamic relationship of the two currents, the grid current tracks the reference very fast, but without the large overshoot, while the PI control ( kp = 0.1, ki = 200) removes the amplitude error. PI has little effect on the quality of idc * and eventually the grid current THD because of the small PI parameters. In fact, because both currents are provided the correct references, they can be controlled more precisely.
- D. Robustness under Grid Voltage Distortion
The actual grid usually contains 5 th and 7 th background harmonics, and this will be even worse for a weak grid. To demonstrate the robustness of the proposed control strategy on grid voltage distortion, simulations are carried out.
Fig. 11 shows the simulation waveforms. The grid voltage is set to have 5% and 3% of the 5 th and 7 th harmonics. Thanks to the high sampling frequency and fast dynamics of the Luenberger observer, the estimated input voltage tracks the actual value accurately despite of the distortion. The estimated value is one-step ahead because of the predictive nature of the observer. With the MPC scheme, the grid current is directly forced to follow a sinusoidal reference, regardless of the grid voltage distortion. Therefore, the grid current waveforms are almost sinusoidal and the 5 th and 7 th harmonic components are only 1.4% and 0.45%, respectively. Obtaining ideal waveforms on both sides when the grid voltage is abnormal is impossible because of the direct coupling of the input and output in matrix converters. However, acceptable performances can be achieved because the cost function involves grid current tracking and DC current tracking at the same time.
PPT Slide
Lager Image
Simulation waveforms under grid voltage distortion.
- E. Parameter Sensitivity
As shown in Equs. (8) and (11), the MPC requires information of the filter parameters to estimate and predict the future quantities. The mismatch of these parameters can degrade the control performance and affect the system stability. Therefore, investigating the control behavior under parameter changes is important. The inductor parasitic resistances Rf and Ro are very small and their modeling errors have negligible influence. The parameter variations of Lf , Cf , and Lo are analyzed through simulations. The actual filter parameters are nominal, while the values used in the predictive controller are considered to change from –30% to +30%.
Fig. 12 plots the grid current THD and the average switching frequency favg , and Fig. 13 shows the waveforms of the grid and DC currents. The system has similar behaviors under the Lf and Cf variations. The underestimation of the two parameters has very small effects on the grid current, whereas the overestimation leads to worse grid current because when Lf and Cf are lower than the actual values, the coefficient of input current term B d (1,1) in Equ. (8) becomes larger, then the prediction of grid current relies more on the input current that is directly associated with the selection of switching states. As a result, the switches operate more frequently for grid current reference tracking. Conversely, when Lf and Cf are higher than the actual values, the role of switching states on the grid current regulation is weakened. This can also be seen in the trends of the average switching frequency shown in Fig. 12 (b). favg reduces greatly when Lf and Cf are overestimated, thus the grid current performance deteriorates significantly. The analysis of Equ. (11) indicates that the mismatch of Lo affects the DC current performance in a similar way. When Lo is higher than its actual value, the DC current is controlled more precisely. However, because of the tradeoff between grid current control and DC current control, the grid current quality is sacrificed. This is verified by Fig. 13 (b), when Lo changes by –30%, the DC current ripple is smaller than that with the case when Lo changes by +30%, but the grid current contains more harmonics. Extreme parameter variations would severely degrade the control performance and even make the system unstable.
PPT Slide
Lager Image
Simulation results under model parameter variations. (a) Grid current THD. (b) Average switching frequency favg.
PPT Slide
Lager Image
Simulation waveforms under model parameter variations. (a) AC filter inductance Lf changes ±30%. (b) DC filter inductance Lo changes ±30%.
VII. CONCLUSION
The MPC scheme for AC/DC matrix converter applied in grid-connected BESS is discussed in this paper. To achieve good grid current regulation, the DC current is controlled as well. By using an improved DC current reference generation method, the grid current can track the reference quickly without a large overshoot or undershoot, and the steady state tracking error is reduced. The AC input voltage is estimated using a closed-loop observer, so the input voltage sensors are eliminated, which reduces the cost and enhances the system reliability. The number of active switching states for cost function evaluation is reduced by half, so the computation time is shortened, the switching loss is reduced, and large dv/dt at the DC terminal is avoided. The proposed MPC strategy is robust against grid voltage distortion. Parameter sensitivity is also investigated. Simulation and experimental results verify that the AC/DC matrix converter with the proposed strategy has good performance for battery charging and discharging control.
Acknowledgements
This work was supported by the National Basic Research Program of China under Award No. 2010CB227206.
BIO
Bo Feng was born in Hubei, China, in 1987. He received a B.S. degree in Electrical Engineering and Automation in 2009 from the Huazhong University of Science and Technology (HUST), Wuhan, China, where he is currently working toward a Ph.D. degree in Power Electronics and Electrical Drive. His research interests include matrix converter and battery energy storage system.
Hua Lin was born in Wuhan, Hubei, China, in 1963. She received a B.S. degree in Industrial Automation from the Wuhan University of Technology, Wuhan in 1984, an M.S. degree in Electrical Engineering from the Naval University of Engineering, Wuhan in 1987, and a Ph.D. degree in Electrical Engineering from the Huazhong University of Science and Technology (HUST), Wuhan in 2005. From 1987 to 1999, she was with the Department of Electrical Engineering, Naval University of Engineering, as a lecturer and associate professor. Since 1999, she has been with the College of Electrical and Electronic Engineering, HUST, where she became a full professor in 2005. In October 2010, she was a visiting scholar at the Center for Advanced Power Systems of Florida State University, Tallahassee. She has been engaged in research and teaching in the field of power electronics and electric drive. Her research interests include high-power high-performance AC motor drives, novel power converters, and their control. She has authored or coauthored more than 50 technical papers in journals and conferences. Dr. Lin was a recipient of the second-grade National Scientific and Technological Advance Prize of China in 1996 and 2003.
References
Miao Z. , Xu L. , Disfani V. R. , Fan L. 2014 “An SOC-based battery management system for microgrids,” IEEE Trans. Smart Grid 5 (2) 966 - 973    DOI : 10.1109/TSG.2013.2279638
Inoue S. , Akagi H. 2007 “A bidirectional DC-DC converter for an energy storage system with galvanic isolation,” IEEE Trans. Power Electron. 22 (6) 2299 - 2306    DOI : 10.1109/TPEL.2007.909248
Holmes D. G. , Lipo T. A. 1992 “Implementation of a controlled rectifier using AC-AC matrix converter theory,” IEEE Trans. Power Electron. 7 (1) 240 - 250    DOI : 10.1109/63.124596
You K. , Xiao D. , Rahman M. F. , Uddin M. N. 2014 “Applying reduced general direct space vector modulation approach of AC/AC matrix converter theory to achieve direct power factor controlled three-phase AC-DC matrix rectifier,” IEEE Trans. Ind. Appl. 50 (3) 2243 - 2257    DOI : 10.1109/TIA.2013.2285956
Julian A. L. , Oriti G. 2005 “A novel clamp circuit for a regenerative rectifier using AC/AC matrix converter theory,” IEEE Trans. Ind. Appl. 41 (1) 68 - 74    DOI : 10.1109/TIA.2004.840963
Ratanapanachote S. , Han J. C. , Enjeti P. N. 2006 “A digitally controlled switch mode power supply based on matrix converter,” IEEE Trans. Power Electron. 21 (1) 124 - 130    DOI : 10.1109/TPEL.2005.861197
Metidji R. , Metidji B. , Mendil B. 2013 “Design and implementation of a unity power factor fuzzy battery charger using an ultrasparse matrix rectifier,” IEEE Trans. Power Electron. 28 (5) 2269 - 2276    DOI : 10.1109/TPEL.2012.2211107
Su M. , Wang H. , Sun Y. , Yang J. , Xiong W. , Liu Y. 2013 “AC/DC matrix converter with an optimized modulation strategy for V2G applications,” IEEE Trans. Power Electron. 28 (12) 5736 - 5745    DOI : 10.1109/TPEL.2013.2250309
Wheeler P. W. , Rodriguez J. , Clare J. C. , Empringham L. , Weinstein A. 2002 “Matrix converters: a technology review,” IEEE Trans. Ind. Electron. 49 (2) 276 - 288    DOI : 10.1109/41.993260
Rodriguez J. , Rivera M. , Kolar J. W. , Wheeler P. W. 2012 “A review of control and modulation methods for matrix converters,” IEEE Trans. Ind. Electron. 59 (1) 58 - 70    DOI : 10.1109/TIE.2011.2165310
Empringham L. , Kolar J. W. , Rodriguez J. , Wheeler P. W. , Clare J. C. 2013 “Technological issues and industrial application of matrix converters: a review,” IEEE Trans. Ind. Electron. 60 (10) 4260 - 4271    DOI : 10.1109/TIE.2012.2216231
Liu X. , Zhang Q. , Hou D. 2014 “One-cycle control strategy with active damping for an AC-DC matrix converter,” Journal of Power Electronics 14 (4) 778 - 787    DOI : 10.6113/JPE.2014.14.4.778
Kouro S. , Cortes P. , Vargas R. , Ammann U. , Rodriguez J. 2009 “Model predictive control – A simple and powerful method to control power converters,” IEEE Trans. Ind. Electron. 56 (6) 1826 - 1838    DOI : 10.1109/TIE.2008.2008349
Rodriguez J. , Kazmierkowski M. P. , Espinoza J. R. , Zanchetta P. , Abu-Rub H. , Young H. A. , Rojas C. A. 2013 “State of the art of finite control set model predictive control in power electronics,” IEEE Trans. Ind. Informat. 9 (2) 1003 - 1016    DOI : 10.1109/TII.2012.2221469
Rivera M. , Rodriguez J. , Wheeler P. W. , Rojas C. A. , Wilson A. , Espinoza J. R. 2012 “Control of a matrix converter with imposed sinusoidal source currents,” IEEE Trans. Ind. Electron. 59 (4) 1939 - 1949    DOI : 10.1109/TIE.2011.2167114
Rivera M. , Rojas C. , Wilson A. , Rodriguez J. , Espinoza J. , Baier C. , Munoz J. 2014 “Review of predictive control methods to improve the input current of an indirect matrix converter,” IET Power Electron. 7 (4) 886 - 894    DOI : 10.1049/iet-pel.2013.0327
Zavala P. , Rivera M. , Kouro S. , Rodriguez J. , Wu B. , Yaramasu V. , Baier C. , Munoz J. , Espinoza J. , Melin P. “Predictive control of a current source rectifier with imposed sinusoidal input currents,” in Conf. IECON 2013 2013 5842 - 5847
Zhang Z. , Wang F. , Sun T. , Rodriguez J. , Kennel R. “FPGA based experimental investigation of a quasi-centralized DMPC scheme for a back-to-back converter,” IEEE Trans. Power Electron. to be published
Cortes P. , Ortiz G. , Yuz J. I. , Rodriguez J. , Vazquez S. , Franquelo L. G. 2009 “Model predictive control of an inverter with output LC filter for UPS applications,” IEEE Trans. Ind. Electron. 56 (6) 1875 - 1883    DOI : 10.1109/TIE.2009.2015750
Cortés P. , Wilson A. , Kouro S. , Rodriguez J. , Abu-Rub H. 2010 “Model predictive control of multilevel cascaded H-bridge inverters,” IEEE Trans. Ind. Electron. 57 (8) 2691 - 2699    DOI : 10.1109/TIE.2010.2041733
Xia C. , Liu T. , Shi T. , Song Z. 2014 “A simplified finite-control-set model-predictive control for power converters,” IEEE Trans. Ind. Inform. 10 (2) 991 - 1002    DOI : 10.1109/TII.2013.2284558
Hu J. , Zhu J. , Lei G. , Platt G. , Dorrell D. G. 2013 “Multi-objective model-predictive control for high-power converters,” IEEE Trans. Energy Convers. 28 (3) 652 - 663    DOI : 10.1109/TEC.2013.2270557
Shi T. , Zhang C. , Geng Q. , Xia C. 2014 “Improved model predictive control of three-level voltage source converter,” Electric Power Components and Systems 42 (10) 1029 - 1038    DOI : 10.1080/15325008.2014.913736
Vargas R. , Rodriguez J. , Rojas C. A. , Rivera M. 2014 “Predictive control of an induction machine fed by a matrix converter with increased efficiency and reduced common-mode voltage,” IEEE Trans. Energy Convers. 29 (2) 473 - 485    DOI : 10.1109/TEC.2014.2299594
You K. , Rahman M. F. 2009 “Analytical model of conduction and switching losses of matrix-Z-source converter,” Journal of Power Electronics 9 (2) 275 - 287
Larrinaga S. A. , Vidal M. A. R. , Oyarbide E. , Apraiz J. R. T. 2007 “Predictive control strategy for DC/AC converters based on direct power control,” IEEE Trans. Ind. Electron. 54 (3) 1261 - 1271    DOI : 10.1109/TIE.2007.893162
Wang X. , Lin H. , She H. , Feng B. 2012 “A research on space vector modulation strategy for matrix converter under abnormal input-voltage conditions,” IEEE Trans. Ind. Electron. 59 (1) 93 - 104    DOI : 10.1109/TIE.2011.2157288