Next Article in Journal
Effect of Manufacturing Inaccuracies on the Wake Past Asymmetric Airfoil by PIV
Previous Article in Journal
A Quantitative Evaluation Method of Anti-Sloughing Drilling Fluid Inhibition for Deep Mudstone
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis and Design of the Energy Storage Requirement of Hybrid Modular Multilevel Converters Using Numerical Integration and Iterative Solution

1
Department of Electrical Engineering, Tsinghua University, Beijing 100084, China
2
Electric Power Research Institute of China Southern Power Grid, Guangzhou 510663, China
*
Author to whom correspondence should be addressed.
Energies 2022, 15(3), 1225; https://doi.org/10.3390/en15031225
Submission received: 29 December 2021 / Revised: 2 February 2022 / Accepted: 3 February 2022 / Published: 8 February 2022

Abstract

:
Increasing the modulation index by utilizing the negative voltage states of full-bridge submodules (FBSMs) can greatly reduce capacitor usage of modular multilevel converters (MMCs), thereby optimizing the cost and volume. The hybrid MMC is composed of half-bridge submodules (HBSMs) and FBSMs, and the capacitor voltages of the two types of submodules (SMs) have different shapes as long as negative voltage states exist. This condition greatly complicates the analysis and design of the energy storage requirement of the hybrid MMC, which utilizes the negative voltage states of FBSMs to boost the AC voltage. A numerical calculation method for solving the capacitor voltages and designing the capacitances of FBSMs and HBSMs is proposed in order to accurately determine the minimum energy storage requirement considering the difference between the energy variations in FBSMs and HBSMs. In the numerical calculation, the energy storage and voltage of the arm are decomposed into FBSM and HBSM parts. According to the physical switching process, the output voltages of FBSM and HBSM parts are determined separately. The one-cycle waveforms of the capacitor voltages are then obtained by numerical integration of the power flows in FBSM and HBSM parts. An iterative solution procedure and the termination criterion that can ensure the accuracy of the obtained one-cycle waveforms are also proposed. Using the numerical integration and iterative solution procedure as the kernel algorithm, the proposed method can accurately analyze the capacitor voltages of the FBSMs and HBSMs and determine the minimum energy storage requirement of the hybrid MMC. Furthermore, the proposed method is applicable for various operating working conditions and various proportions of FBSMs. The simulation results verify the feasibility and accuracy of the analysis and design method.

1. Introduction

Modular multilevel converters (MMCs) have been widely used in voltage-source-converter-based high-voltage direct-current (VSC-HVDC) transmission. Compared with the two- or three-level converters, MMCs have a modular structure and can achieve higher output voltages with lower harmonics and losses [1,2,3,4,5]. Each submodule (SM) has a large capacitor to maintain the capacitor voltages and limit the capacitor voltage ripple within an acceptable level. However, because a large energy variation in each arm of an MMC exists, the required SM capacitance is relatively large. In each SM, the capacitor may account for approximately 30% of the cost and 50% of the volume [6], resulting in a cost and volume problem for the MMC. Thus, reducing capacitor usage has become a research focus.
Existing studies generally consider two technical routes to reduce the SM capacitance. The first route reduces the energy variation in the arm through circulating current injection [7,8,9,10]. However, the capacitance reduction effect is limited because the injected circulating current also increases the root-mean-square (RMS) and peak values of the arm current, which reduces the safe operating margin [9,10]. The other route increases the modulation index of the MMC to a higher value. Previous studies have found that this method can greatly reduce both the capacitor usage and the rated AC current simultaneously [11,12,13,14,15,16,17,18]. For example, literature [11] found that the fundamental-frequency capacitor ripple can be eliminated by increasing the modulation index to 1.414 under a pure active power condition, greatly reducing the required capacitance. Literature [12,13] considered the influence of reactive power on the fundamental-frequency capacitor ripple and the double-frequency ripple. In addition, some studies combined the approach of increasing the modulation index with a double-frequency circulating current injection, third harmonic voltage injection, and other methods to further reduce the capacitor voltage ripple and capacitor usage [14,15,16,17,18].
The half-bridge submodule (HBSMs)-based MMC (HB-MMC) is the most commonly used topology. However, the maximum attainable modulation index of the HB-MMC is limited. To increase the modulation index to such a high value that the energy variation and energy storage requirement can be reduced, an MMC composed of FBSMs (FB-MMC) is used. The capacitor voltage ripple and the required SM capacitance can be obtained through the overall energy stored in the arm, which can be obtained by integrating the power flow in the arm. The analysis approach of the overall arm energy, which is based on the assumption that all of the SMs in each arm equally share the stored energy all the time, is only applicable to an FB-MMC that contains only one type of SM. However, the semiconductor usage of the FB-MMC is extremely high, and the resulting cost and volume of the MMC are even larger than those of a conventional MMC, even if the energy storage requirement is minimized by increasing the modulation index.
Another approach is adding several extra full-bridge SMs (FBSMs) in each arm (utilizing their negative voltage states) to the conventional HB-MMC to form a hybrid MMC [19], which can lower the increase in semiconductor usage. In this scenario, the rated DC voltage remains unchanged, and the rated AC voltage and modulation index can be increased linearly with the increase in FBSMs. However, exactly evaluating and designing the energy storage requirement of a hybrid MMC that utilizes negative voltage states is greatly challenging. The hybrid MMC contains FBSMs and HBSMs in each arm and the energy variation and the capacitor voltages of the two types of SMs are significantly different if the arm voltage has a period of negative states [20,21,22]. This condition is mainly because only the FBSMs are charged or discharged, whereas all of the HBSMs are bypassed during the negative-voltage period. Therefore, the shape of the capacitor voltage ripple and the required SM capacitance of the two types of SMs are relatively different in a hybrid MMC, and the analysis approach through the overall arm energy is inapplicable. Although the reduction effect of the cost and size of hybrid MMC using negative voltage states is analyzed in [13], the difference between the energy storage variation of the two types of SMs is not considered, and the analysis is still based on the overall arm energy—the accuracy of which can only be ensured to the MMCs, whose arms are composed of the same type of SMs. Some research efforts have been conducted on the analysis and design method of the energy storage requirement of hybrid MMCs [14,15,16]. The accuracy and applicability of these methods for various operating conditions should be further explored because they are based on some assumptions that are only true for some special conditions. The charging and discharging processes of the capacitors in HBSMs and FBSMs are analyzed physically using piecewise analysis in [17]; they are mainly focused on the working condition with a power factor of no less than 0.9 and a proportion of FBSMs equal to 2/3.
In this study, a numerical calculation method for solving the capacitor voltages and designing the capacitances of FBSMs and HBSMs is proposed in order to provide a method that can accurately determine the minimum energy storage requirement, considering the difference between the energy variations in FBSMs and HBSMs. First, the energy storage and voltage of the arm are decomposed into FBSM and HBSM parts, and the output voltage of these parts is determined separately according to the physical switching process. A numerical integration method is then adopted to obtain the one-cycle waveforms of the capacitor voltages. To solve the problem that the initial values of the capacitor voltages of the FBSMs and HBSMs cannot be exactly determined, an iterative solution procedure and the termination criterion that can ensure the accuracy of the obtained one-cycle waveforms are proposed. Then, based on the numerical integration and iterative solution procedure as the kernel algorithm, the method for finding the minimum energy storage requirement of the hybrid MMC and determining the capacitances of the FBSMs and HBSMs is proposed. The simulation results verify the feasibility and accuracy of the analysis and design method.

2. Hybrid MMC Utilizing Negative Voltage States of FBSMs

Previous studies have found that increasing the modulation index can significantly reduce the capacitor usage of MMCs. However, the arm of the conventional HB-MMC cannot output negative voltage states, so the modulation index is limited to 1.0. To achieve a high modulation index while lowering the increase in semiconductor usage, several extra FBSMs can be appended to each arm of the HB-MMC to form a hybrid MMC. As shown in Figure 1, an HB-MMC with N0 submodules is appended with F FBSMs to form a hybrid MMC. In Figure 1, phase A is considered for illustration, where uap and uan are the arm voltages, iap and ian are the arm currents, and Udc is the DC voltage.
Ignoring the losses, the equivalent circuit of the MMC connected to the AC grid is shown in Figure 2. The AC grid voltage phasor is denoted using Us = Us∠0, where Us is the rated RMS value of the phase voltage. The phasor of the AC port voltage of the MMC is denoted using Uac = Uacδ. According to the equivalent circuit shown in Figure 2, the phasor of the AC port voltage can be expressed as follows:
U ac = U s + j X I ac ,
where Iac = Iac∠-φ is the phasor of the output AC current, Iac is the RMS value of the AC current, and φ is the power factor angle. X is the equivalent reactance between the MMC and the AC grid, including the line reactance between the MMC and PCC, transformer leakage reactance, and arm reactance. When outputting the AC current Iac∠-φ, the RMS value and phase angle of the AC port voltage of the MMC can be expressed as follows:
U ac = ( U s + X I ac sin φ ) 2 + ( X I ac cos φ ) 2 U s ( 1 + X I ac sin φ U s ) ,
δ = arctan X I ac cos φ U s + X I ac sin φ .
Therefore, at the operating point when outputting the AC current Iac∠-φ, the corresponding modulation index Mac is expressed as follows:
M ac = 2 U ac U dc / 2 = 2 ( U S + X I ac sin φ ) U dc / 2 = M 0 ( 1 + X * I ac * sin φ ) ,
where X* is the per-unit value of the equivalent reactance, I*ac is the per-unit value of the ac-side current, and M0 is the base modulation index when the MMC outputs zero current, and is expressed as follows:
M 0 = 2 U S U dc / 2 .
The base modulation index M0 also represents the per-unit AC grid voltage normalized by Udc/2.
Figure 3 shows an example of the relation between the number of extra FBSMs and the MMC arm voltage. The rated voltage of the capacitor is denoted as Uc and can be expressed as follows:
U c = U dc N 0 .
If the number of the extra FBSMs in each arm is denoted using F, then the range of the arm voltage is [−FUc, N0 + FUc]. To satisfy the linear modulation constraint, the peak and valley values of the arm voltages at any operating point should be within the range [−FUc, N0 + FUc], which can be expressed as follows:
{ M ac + 1 2 U dc < ( N 0 + F ) U c , F U c > 1 M ac 2 U dc .
According to Equation (7), the number of extra FBSMs should satisfy the following condition:
F N 0 > M ac 1 2 .
Evidently, the maximum modulation index of all possible operating points determines the minimum required number of extra FBSMs. The operating range of MMC can be described using the PQ diagram, as shown in Figure 4. The main function of the MMC in the VSC-HVDC system is to transfer the active power. Existing studies have demonstrated that the capacitor voltage ripple under a high capacitive power condition is much greater than a pure active power condition. Outputting a reactive power of rated capacity in actual applications is not required, and a maximum reactive power requirement, which is denoted using ±Qmax, is usually set according to the requirement of the connected AC grid. Therefore, the operating region of the MMC is shown as the shaded area in Figure 4. According to Equation (4), the maximum modulation index within the given operating region can be expressed as follows:
M ac _ max M 0 ( 1 + X * Q max ) ,
where the maximum modulation index is related to M0 and Qmax given that the value of X* is certain. Then, combining Equations (8) and (9), the required minimum number of the extra FBSMs is determined by the base modulation index and maximum reactive power requirement, as follows:
F min = M ac _ max 1 2 N 0 .

3. Algorithm for Solving the One-Cycle Waveforms of the Capacitor Voltages of FBSMs and HBSMs

The peak values of the capacitor voltages for a certain operating point should be exactly calculated, because the capacitances of the SMs should be designed to satisfy the constraint that the peak values of the capacitor voltages be below the permitted upper limit under all possible operating conditions. A physical switching process-based numerical calculation method for solving the one-cycle waveforms of the capacitor voltage is proposed: The energy storage and voltage of the arm are decomposed into FBSM and HBSM parts, and the output voltages of these parts are determined separately according to the actual switching process; then, a numerical integration and iterative solution procedure are proposed to obtain the accurate one-cycle waveforms of the capacitor voltages of the FBSMs and FBSMs. The obtained accurate one-cycle waveforms can be used to determine the peak values of the capacitor voltages of the FBSMs and HBSMs and further provide a base for finding the minimum capacitor usage of the hybrid MMC.

3.1. Arm Energy and Voltage Decomposition for Hybrid MMC

For the hybrid MMC, the capacitor voltages of the HBSMs and FBSMs should be calculated separately. The capacitor voltage sorting algorithm is assumed to achieve a good balance among all the FMSMs and among all the HBSMs. Therefore, the FBSM capacitor voltages can be denoted using a single variable ucf(t), and the HBSM capacitor voltages can be denoted using a single variable uch(t). The capacitor voltages of ucf(t) and uch(t) may not be equal; thus, directly solving the capacitor voltage through the overall stored energy in each arm is impossible. However, the energy stored in each arm can be decomposed into two parts of FBSMs and HBSMs. For each part, the same relation between stored energy and the capacitor voltages still holds.
The required capacitance of the FBSMs and HBSMs is also different because the FBSMs and HBSMs in a hybrid MMC have different energy variations. The HBSM capacitance is denoted as Cd. The ratio of the capacitance of the FBSMs to that of the HBSMs is the capacitance ratio kf, and the FBSM capacitance is kfCd. Under the operating conditions when the modulation index is greater than 1.0, the capacitor voltages of the FBSMs fluctuate greater than those of the HBSMs, thereby increasing the total capacitor usage requirement. By appropriately increasing the capacitance ratio, the voltage ripple of the FBSMs can be reduced to decrease the overall capacitor usage of the MMC. The calculation method considers this factor.
The instantaneous energy stored in each arm can then also be decomposed into FBSM and HBSM parts, expressed as follows:
{ e f ( t ) = F 1 2 ( k f C d ) u ch 2 ( t ) , e h ( t ) = N 0 1 2 ( C d ) u cf 2 ( t ) .
The capacitor usage can also be represented using the energy storage requirement Enom, which is defined as the normalized energy that is stored in all of the SM capacitors when the SM capacitor voltage is its nominal value [11,12]. In hybrid MMCs, the total energy storage requirement can be decomposed into the energy storage requirements of FBSMs and HBSMs, which are denoted as Enomf and Enomh, respectively, as follows:
{ E nomf = 6 × 1 2 ( k f C d ) F U c 2 S N = k f F 1 + k f F E nom , E nomh = 6 × 1 2 C d N 0 U c 2 S N = 1 1 + k f F E nom ,
where SN is the rated capacity of the MMC, and Enom is the total energy storage requirement, which can be expressed as follows:
E nom = E nomf + E nomh .
To calculate the capacitor usage of the FBSM and HBSM parts, the arm voltages should also be decomposed into two parts, whose voltages can be determined according to the sorting algorithm. The commonly used sorting algorithm is shown in Figure 5. The output voltages of the FBSM and HBSM parts can be obtained according to the relation among arm voltage, arm current, capacitor voltages, and the number of inserted SMs.
Using the upper arm voltage of phase A as an example, the arm voltage and current can be expressed as follows:
u arm ( t ) = U dc 2 2 U ac sin ( ω t + δ ) ,
i arm ( t ) = 1 3 I dc + 2 2 I ac sin ( ω t φ ) ,
where Idc is the DC current of the MMC. According to the active power balance between the AC and DC sides of the MMC, the DC current of the MMC can be expressed by AC current as follows:
I dc = 3 2 2 M 0 I ac cos ( φ ) .
Figure 6 shows an example of the variation process of the instantaneous capacitor voltages in one fundamental-frequency cycle when the MMC operates under a pure active absorption condition. The negative zero-crossing point of the arm voltage is selected as the origin of the time axis. According to the zero-crossing point of the arm voltage and current, the energy variation of the arm in a fundamental-frequency cycle can be divided into four phases:
(1)
Phase 1 (charge)
In this phase, the arm voltage is negative, and only FBSMs participate in voltage output. Therefore, the voltages of the FBSM and HBSM parts can be expressed as follows:
{ u f ( t ) = u arm ( t ) , u h ( t ) = 0 .
According to nearest-level modulation (NLM), the actual number of the inserted SMs is consistent with the arm voltage. Figure 6d shows the insertion indices of FBSMs and HBSMs. In this phase, only FBSMs are inserted and HBSMs are bypassed. As shown in Figure 6a, the capacitor voltages of the FBSMs increasingly diverge from those of the HBSMs, and reach a peak at the end of this phase because only the FBSMs are charged in this phase.
(2)
Phase 2 (discharge)
The arm voltage becomes positive and the arm current is negative. As a result, the capacitors are discharged in this phase. According to Figure 6a, the capacitor voltages of FBSMs are greater than those of the HBSMs at the beginning of the second phase due to the first charging phase. Therefore, in this phase, the switching strategy initially inserts as many FBSMs as possible under the constraints that the insertion index of FBSMs does not exceed F to discharge their capacitor voltages. The HBSMs then output the remaining arm voltage correspondingly. The voltages of the FBSM and HBSM parts in this phase are expressed as follows:
{ u f ( t ) = min ( u arm ( t ) , F U c ) , u h ( t ) = u arm ( t ) u f ( t ) .
As shown in Figure 6d, the insertion index of FBSMs in the first half of the second phase gradually increases with the demand for arm voltage and finally saturates to F; the remaining arm voltage is generated by HBSMs. As shown in Figure 6a, the capacitor voltages of the FBSMs and HBSMs start to converge gradually in this phase.
(3)
Phase 3 (charge)
In this phase, the arm voltage and current are positive, and the capacitors are charged in this phase. As shown in Figure 6a, at the beginning of this phase, as many HBSMs as possible should be initially inserted when the capacitor voltages of FBSMs are still higher than those of the HBSMs. The voltages of the FBSM and HBSM parts in this phase are expressed as follows:
{ u h ( t ) = min ( u arm ( t ) , N 0 U c ) , u f ( t ) = u u h ( t ) .
Subsequently, the capacitor voltages of the FBSMs and HBSMs are equal. At this time, the arm-switching strategy balances the two capacitor voltages. To maintain equal capacitor voltages, the energy absorbed by the FBSMs and HBSMs should be proportional to their nominal energy storage. The output voltages of the FBSM and HBSM parts should be proportional to their nominal energy storages because the current flowing through the arm is equal, as follows:
{ u f ( t ) = E nomf E nom u arm ( t ) , u h ( t ) = E nomh E nom u arm ( t ) .
As shown in Figure 6d, in the first half of the third phase, as many HBSMs are inserted as possible. In the second half, the insertion indices of FBSMs and HBSMs are proportional.
(4)
Phase 4 (discharge)
In this phase, the arm voltage is positive and the current is negative, resulting in discharging the capacitors. The same relation of (20) exists because the capacitor voltages of the FBSMs and HBSMs are equal.
This analysis is based on the example of pure active absorption condition. However, the proposed method is not limited to this circumstance. In general, the voltages of FBSM and HBSM parts based on the switching process are summarized in Table 1.
After obtaining the output voltages of the FBSM and HBSM parts, the instantaneous power flowing into the two parts can be obtained by the product of the arm current with the respective voltages of the two parts, as follows:
{ p f ( t ) = u f ( t ) i arm ( t ) , p h ( t ) = u h ( t ) i arm ( t ) .
The instantaneous energy stored in the FBSM and HBSM parts can then also be expressed by integrating the corresponding instantaneous power, as follows:
{ e f = p f d t = u f ( t ) i arm ( t ) , e h = p h d t = u h ( t ) i arm ( t ) .
Referring to Equation (11), the capacitor voltages (p.u.) of the FBSM and HBSM parts can then be calculated as follows:
{ u cf * ( t ) = 6 e f ( t ) S N E nomf , u ch * ( t ) = 6 e h ( t ) S N E nomh .

3.2. Numerical Integration and Iterative Solution Procedure

Based on the decomposed arm energy and voltage, a numerical integration and iterative solution procedure is proposed. The numerical calculation aims to obtain the one-cycle waveforms of the capacitor voltages under a certain operating condition and thereby obtains the peak values of the capacitor voltages for capacitance design. The flowchart of the numerical integration and iterative solution procedure is shown in Figure 7. The detailed descriptions of the main steps of the numerical calculation process for obtaining the one-cycle waveforms are as follows:
(1)
Initial values of the numerical calculation.
The initial values should be set for the numerical calculation. As analyzed in Section 3.1, the numerical calculation starts at the negative zero-crossing point of the arm voltage. From the view of the overall energy, if the initial value of the overall energy is set as e(0), then the energy stored in the overall arm can be expressed as follows:
e ( t ) = e ( 0 ) + 0 t u arm ( τ ) i arm ( τ ) d τ .
Note that the Volterra operator is adopted, which has also gained applications for storage system modeling and control [23]. Under steady-state operation, the average value of the overall arm energy storage in one fundamental-frequency cycle should be the nominal energy storage, and the following equation can be obtained:
1 T 1 0 T 1 e ( t ) d t = e ( 0 ) + 1 T 1 0 T 1 ( 0 t u arm ( τ ) i arm ( τ ) d τ ) d t = S N E nom 6 ,
where T1 is the fundamental-frequency period. Therefore, the initial value of the overall arm energy can be determined as follows:
e ( 0 ) = S N E nom 6 1 T 1 0 T 1 ( 0 t u arm ( τ ) i arm ( τ ) d τ ) d t .
The initial value of the overall arm energy can be determined in advance of the numerical calculation because the arm voltage and current can be determined as shown in (14) and (15) in advance.
However, the numerical calculation needs the initial values of the capacitor voltages of the FBSMs and HBSMs. Considering that the ripple component is related to the power factor angle of the operating point, the values of the first time point of the one-cycle waveforms of the capacitor voltages are generally not equal to the DC rated value. Furthermore, the values of the first time point of the capacitor voltages of the HBSMs and FBSMs may be different. This condition brings great difficulties in determining the initial values of the capacitor voltages. To start the numerical calculation, the initial values of the capacitor voltages of FBSMs and HBSMs are assumed to be equal and are determined by the initial overall arm energy, as follows:
u cf * ( 0 ) = u ch * ( 0 ) = 6 e ( 0 ) S N E nom
This assumption only aims to initiate the numerical calculation, and it does not indicate that the values of the first time points of the capacitor voltages of FBSMs and HBSMs always appear equal in the numerical calculation. The values of the first time points of the one-cycle waveforms are corrected iteratively in step (3).
(2)
Numerical integration to obtain temporary one-cycle waveform.
For a certain operating point, the arm voltage and current can be calculated, as shown in Equations (14) and (15). According to the analysis in Section 3.2, the output voltages of the FBSM and HBSM parts can be determined referring to Table 1. The discretized instantaneous power of the two parts at the k-th time-step point can then be obtained as follows:
{ p f ( k ) = u f ( k ) i arm ( k ) , p h ( k ) = u h ( k ) i arm ( k ) .
The energy storages of the FBSM and HBSM parts at the next time-step point can be obtained through numerical integration of the respective instantaneous power. A left-hand rectangle integration formula is adopted for its simplicity and accuracy in small time steps. It satisfies
{ e f ( k + 1 ) = e f ( k ) + p f ( k ) Δ t , e h ( k + 1 ) = e h ( k ) + p h ( k ) Δ t ,
where ∆t is the time-step of the numerical calculation. According to the relationship between the energy storage and capacitor voltage, the capacitor voltages of the FBSM and HBSM parts at the next time-step point can be expressed as follows:
{ u cf ( k + 1 ) = 6 e f ( k + 1 ) S N E nom f , u ch ( k + 1 ) = 6 e h ( k + 1 ) S N E nom h .
The calculation for the next time step point is then performed until the calculation of a fundamental-frequency period is completed and the temporary one-cycle waveforms of the capacitor voltages are obtained.
(3)
Iterative solution procedure for obtaining the final one-cycle waveforms.
At step (1), the initial voltages of the FBSM and HBSM parts are assumed to be equal, which may be not true when the capacitance ratio of FBSMs and HBSMs is not 1. As a result, a single calculation using step (2) cannot guarantee the accuracy of the obtained one-cycle waveforms of the capacitor voltages. To guarantee accuracy, an iterative solution procedure is proposed to obtain the final accurate one-cycle waveforms. In the iterative solution procedure, the calculation for solving one-cycle waveforms in step (2) is performed iteratively, and the values of the final time-step point of the obtained one-cycle waveforms are used as the initial values of the next iterative calculation. Based on the physical process of capacitor voltage balancing which drives the capacitor voltages of FBSMs and HBSMs equal, the values of the first and final time-step points of the accurate one-cycle waveforms will converge because the waveforms of the capacitor voltages are periodic. Therefore, the termination criterion of the iterative solution procedure is that the relative error between the values of the first and final time-step points of the obtained one-cycle waveforms is less than a preset tolerance. Moreover, it must be mentioned that the convergence of capacitor voltages of the hybrid MMC is related to the hybrid ratio (F/N0), modulation index, and AC power factor [20]. As the main focus of this paper is the calculation method of capacitor voltages and energy storage design of hybrid MMCs, detailed convergence analysis is omitted.
Note that a very low tolerance will increase iteration time, as the capacitor voltage needs many cycles to converge under this tolerance. On the contrary, a high tolerance will cause inaccurate results, as the capacitor voltages may not reach steady states. Thus, a moderate tolerance of 0.1% is set in order to balance the accuracy and speed in this paper. The 0.1% tolerance is relatively accurate, as it is much lower than the 10% capacitor voltage ripple to be calculated. On the other hand, the iteration process only takes a few cycles under the 0.1% tolerance in our calculation, which improves the calculation’s efficiency. Using the iterative solution procedure, the final accurate one-cycle waveforms of the capacitor voltages can be obtained after a few iterations.

4. Method for Finding the Minimum Required Capacitances

The capacitances of the FBSMs and HBSMs should be selected to ensure that the peak values of the capacitor voltages of the FBSMs and HBSMs are below the permitted upper limit under all possible points. To minimize the cost and volume, the minimum capacitor usage (energy storage requirement) and the respective capacitances of the FBSMs and HBSMs that can satisfy this constraint should be determined. For a certain operating point and a certain set of parameters, the peak value of the capacitor voltages can be determined according to the obtained one-cycle waveforms. Therefore, using the numerical integration and iterated solution procedure presented in Section 3 as the kernel algorithm, a method for finding the minimum capacitor usage and the respective capacitances of the FBSMs and FBSMs is proposed.
Prior to finding the minimum required capacitor usage, the base modulation index M0, the P/Q range, and the permitted upper limit of the capacitor voltage εlim, which are the key factors that affect the energy variation and capacitor usage of the MMC, should be specified initially. The maximum reactive power Qmax can be specified, and the P/Q range is depicted as the shaded area in Figure 4.
One of the complexities of the hybrid MMC is that the capacitances of the FBSMs and HBSMs may not be equal. Therefore, the total capacitor usage Enom and the capacitance ratio kf should be considered the search variables. Figure 8 shows the flowchart to find the minimum capacitor usage for a certain configuration (M0 and P/Q range). The searching algorithm includes outer and inner loops. The outer-loop variable is the capacitance ratio kf. A reasonable range of kf, such as [1.0, 4.0], can be set for the outer-loop search. For each value of kf, the maximum value of Enom_i that allows all operating points within the given P/Q range to satisfy the constraint of peak values of the capacitor voltages is found using the dichotomy method. Finally, all the obtained values of Enom_kf are compared, and the minimum value is selected as the final capacitor usage Enom; the corresponding capacitance ratio kf is also determined. According to Equations (12) and (13), the energy storage requirements and capacitances of the FBSMs and HBSMs are also determined.
The curves of the energy storage requirement with the change in M0 and Qmax can also be obtained by varying the values of M0 and Qmax.

5. Case Study and Simulation

To verify the accuracy of the proposed numerical calculation method, an electromagnetic simulation model of a hybrid MMC was built in MATLAB/Simulink, where all parts of the simulation model were kept consistent with the real application, e.g., main circuit topology and parameters, converter control, pulse generation, capacitor voltage balancing, etc. As shown in Figure 9, the AC side of the hybrid MMC is connected with an AC grid, and the DC side is connected with a DC current source. The hybrid MMC operates under the Udc/Q control mode. The DC voltage is controlled as the rated value, and the active power can be adjusted by the DC current. The reactive power is controlled by the MMC. Simulations have been conducted under different scenarios and show high accuracy with the proposed numerical calculation method. Due to the limitation of paper length, only one typical case is considered as an example.

5.1. Verification of the Theoretical Calculation Algorithm under Certain Values of M0 and Qmax

First, the case of M0 = 1.2 and Qmax = 1.0 p.u. is considered an example. The permitted upper limit of the capacitor voltage is set to 1.1 p.u. Table 2 lists the obtained parameters using the design method presented in Section 3 and Section 4.
The algorithm presented in Section 3 can find the one-cycle waveforms of the capacitor voltages of FBSMs and HBSMs for any given operating points. Figure 10 illustrates the obtained one-cycle waveforms of the capacitor voltages for several typical operating points, such as pure active power input, pure active power output, pure inductive reactive power, and pure capacitive reactive power operating conditions. Consistent with the theoretical analysis, the waveforms of the capacitor voltages of the FBSMs and HBSMs have different shapes as long as negative voltage states exist. The capacitor voltages of the FBSMs and HBSMs can be accurately calculated based on the proposed arm voltage decomposition method and numerical integration and iterative solution procedure. The simulation results are also illustrated in Figure 10. The figure shows that the simulation and theoretical calculation results fit effectively, thereby verifying the accuracy of the proposed calculation method.
The peak values of the capacitor voltages can be determined according to the obtained one-cycle waveforms. Under the constraint that the peak values of the capacitor voltages of the FBSMs and HBSMs are below the permitted upper limit under all possible points, the minimum capacitor usage and the corresponding capacitances of the FBSMs and FBSMs can be found using the method presented in Section 4. In this case study, when M0 = 1.2, Qmax = 1.0 p.u. and εlim = 1.1 p.u., the search result of the minimum energy storage requirement is 35.7 kJ/MVA. Correspondingly, the capacitances of the FBSMs and FBSMs are 18.2 and 14 mF, respectively, and the capacitance ratio kf is 1.3. In this example, the capacitor voltage ripple under the rated capacitive reactive power condition (φ = π/2) is the highest, because Qmax is set to 1.0 p.u. and it is relatively high. As illustrated in Figure 10d, using the designed capacitances, the peak values of the capacitor voltages of the FBSMs and HBSMs are relatively close to 1.1 p.u., which also verifies the feasibility of the method for finding the minimum capacitor usage and the capacitances of the FBSMs and HBSMs.

5.2. Analysis of the Impacts of the Base Modulation Index and Maximum Reactive Power on Capacitor Usage and the Number of Extra FBSMs

The base modulation index M0, which also represents the normalized rated AC grid voltage, is a key parameter that determines capacitor usage and semiconductor usage. In the existing studies, in which a FB-MMC that only contains FBSMs is considered, the capacitor usage can reach the minimum when the base modulation index is increased to approximately 1.414. However, this conclusion may not be applicable to a hybrid MMC that contains two types of SMs, because determination of the capacitor usage is relatively different from the FB-MMC. The maximum reactive power requirement Qmax is another key parameter that affects capacitor usage. Generally, capacitor usage increases with the increase in the maximum reactive power requirement.
For certain values of M0 and Qmax, the minimum required capacitor usage can be found using the method presented in Section 4; then, the impacts of the base modulation index and maximum reactive power on capacitor usage can be observed. Figure 11 shows the curve of the capacitor usage Enom, following the change in the base modulation index M0. To illustrate the impacts of the maximum reactive power on capacitor usage, three typical cases, namely, Qmax = 1.0 p.u., Qmax = 0.5 p.u., and Qmax = 0 p.u., are shown as examples. In addition, the curves of the capacitor usage of the FB-MMC are shown in Figure 11 for comparison.
Figure 11 shows that increasing M0 can decrease capacitor usage for the FB-MMC and the hybrid MMC. In addition, decreasing the maximum reactive power requirement can reduce the capacitor usage.
However, Figure 11 also reveals the difference in capacitor usage between the hybrid MMC and the FB-MMC when M0 is increased further. This concept is overlooked by previous studies when evaluating the cost and volume. As shown in Figure 11, when the value of Qmax is large, capacitor usage between the hybrid MMC and FB-MMC is slightly different. However, when the value of Qmax is low, capacitor usage of the hybrid MMC starts to increase when M0 is increased above an inflection point (approximately 1.0), which is much lower than the inflection point (approximately 1.414) of FB-MMC. This finding is caused by the low Qmax, where the capacitor usage is determined by the rated active power operating point, and the hybrid MMC has the highest capacitor voltage imbalance under higher M0. Thus, the capacitor usage of the hybrid MMC is significantly higher than for the FB-MMC. However, when Qmax is higher, the capacitor usage is determined by the rated capacitive reactive power operating point because of its inherent large capacitor voltage ripple, even if the capacitor voltage imbalance in the active power operating point is considered. Thus, the capacitor usage in this case differs slightly from the FB-MMC.

6. Discussion on Advantages of the Proposed Method and Conclusions

Hybrid modular multilevel converters are a promising choice for future flexible high-voltage DC transmission. In applications of boosted AC voltage or reduced DC voltage, the capacitor voltage of FBSMs and HBSMs differs because of the negative voltage states of FBSMs. It significantly increases the complexity of capacitor voltage ripple calculation and energy storage design. Several methods have been proposed to solve this problem, and a numerical integration and iterative solution method is proposed in this paper.
Table 3 is listed below to evaluate the above methods qualitatively in three dimensions: accuracy, completeness, and low complexity. Accuracy reflects the simplification degree of the model, i.e., totally neglecting the capacitor voltage difference between FBSMs and FBSMs results in the simplest model but the least accurate; Completeness reflects the application scope of the model, i.e., some calculation methods only consider a limited active power region, a certain hybrid ratio, or a low modulation index. Thus, the model is only applicable in a certain scenario; Low complexity reflects the analytical complexity of the model, i.e., the piecewise calculation method involves considerable discussion and calculation, which is quite complicated in real applications. The comparison table clearly shows the advantages of the proposed method to the previous methods.
For the proposed method, the energy storage and voltage of the arm are decomposed into FBSM and HBSM parts in the numerical calculation. According to the physical switching process, the output voltage of the FBSM and HBSM parts are described separately. The one-cycle waveforms of the capacitor voltages are then obtained by numerical integration of the power flows in the FBSM and HBSM parts. Furthermore, using the termination criterion, where the relative error between the values of the first and final time-step points of the obtained one-cycle waveforms is less than a preset tolerance, the proposed iterative solution procedure can ensure the accuracy of the final obtained one-cycle waveforms of the capacitor voltages; then, using the numerical integration and iterative solution procedure as the kernel algorithm, the method for finding the minimum capacitor usage and the capacitances of the FBSMs and FBSMs is also proposed. The feasibility and accuracy of the analysis and design method of the energy storage requirement of hybrid MMC are verified using the simulation results.

Author Contributions

Conceptualization, K.W.; methodology, K.W.; software, K.W.; validation, K.W.; formal analysis, K.W.; investigation, K.W.; resources, K.W.; data curation, K.W.; writing—original draft preparation, K.W.; writing—review and editing, Q.S.; visualization, K.W.; supervision, Q.S.; project administration, Q.S. and S.X.; funding acquisition, Q.S. and S.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 51977119.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

N0Base number of HBSM.
FNumber of FBSM appended.
uapArm voltage of phase A upper arm.
uanArm voltage of phase A lower arm.
iapArm current of phase A upper arm.
ianArm current of phase A lower arm.
UdcDC voltage.
IdcDC current.
UsAC grid RMS phase voltage.
UacAC port RMS phase voltage of MMC.
δAC port phase voltage phase of MMC.
IacAC RMS current of MMC.
φAC current phase of MMC.
XEquivalent reactance between MMC and AC grid.
MacModulation index.
M0Base modulation index.
UcBase capacitor voltage.
QmaxMaximum reactive power of MMC.
ucfCapacitor voltage of FBSM.
uchCapacitor voltage of HBSM.
CdCapacitance of HBSM.
kfCapacitance ratio of FBSM to HBSM.
efArm capacitor energy of FBSM.
ehArm capacitor energy of HBSM.
EnomEnergy storage requirement of MMC.
EnomfEnergy storage requirement of FBSM.
EnomhEnergy storage requirement of HBSM.
SNCapacity of MMC.
uarmArm voltage.
iarmArm current.
ufArm voltage of FBSM.
uhArm voltage of HBSM.
pfArm power of FBSM.
phArm power of HBSM.
T1Fundamental-frequency period.
εlimUpper limit of capacitor voltage.
Enom_iEnergy requirement at operating point i.
Enom_kfEnergy requirement at capacitance ratio kf.
Superscript *Per unit value.
Functional spaceFundamental-frequency period: [0, T1].

References

  1. Lesnicar, A.; Marquardt, R. An innovative modular multilevel converter topology suitable for a wide power range. In Proceedings of the IEEE Bologna Power Tech Conference Proceedings, Bologna, Italy, 23–26 June 2003; Volume 3, p. 6. [Google Scholar] [CrossRef]
  2. Hagiwara, M.; Akagi, H. Control and Experiment of Pulsewidth-Modulated Modular Multilevel Converters. IEEE Trans. Power Electron. 2009, 24, 1737–1746. [Google Scholar] [CrossRef]
  3. Tu, Q.; Xu, Z.; Xu, L. Reduced Switching-Frequency Modulation and Circulating Current Suppression for Modular Multilevel Converters. IEEE Trans. Power Deliv. 2011, 26, 2009–2017. [Google Scholar]
  4. 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]
  5. Han, X.; Sima, W.; Yang, M.; Li, L.; Yuan, T.; Si, Y. Transient Characteristics under Ground and Short-Circuit Faults in a ±500 kV MMC-Based HVDC System with Hybrid DC Circuit Breakers. IEEE Trans. Power Deliv. 2018, 33, 1378–1387. [Google Scholar] [CrossRef]
  6. Davidson, C.; Trainer, D. Innovative concepts for hybrid multi-level converters for HVDC power transmission. In Proceedings of the 9th IET International Conference on AC and DC Power Transmission (ACDC 2010), London, UK, 19–21 October 2010; pp. 1–5. [Google Scholar]
  7. Picas, R.; Pou, J.; Ceballos, S.; Agelidis, V.G.; Saeedifard, M. Minimization of the capacitor voltage fluctuations of a modular multilevel converter by circulating current control. In Proceedings of the IECON 2012—38th Annual Conference on IEEE Industrial Electronics Society, Montreal, QC, Canada, 25–28 October 2012; pp. 4985–4991. [Google Scholar] [CrossRef]
  8. Picas, R.; Pou, J.; Ceballos, S.; Zaragoza, J.; Konstantinou, G.; Agelidis, V.G. Optimal injection of harmonics in circulating currents of modular multilevel converters for capacitor voltage ripple minimization. In Proceedings of the 2013 IEEE ECCE Asia Downunder, Melbourne, Australia, 3–6 June 2013; pp. 318–324. [Google Scholar]
  9. Pou, J.; Ceballos, S.; Konstantinou, G.; Agelidis, V.G.; Picas, R.; Zaragoza, J. Circulating current injection methods based on instantaneous information for the modular multilevel converter. IEEE Trans. Ind. Electron. 2015, 62, 777–788. [Google Scholar] [CrossRef] [Green Version]
  10. Li, X.; Song, Q.; Liu, W.; Xu, S.; Zhu, Z.; Li, X. Performance Analysis and Optimization of Circulating Current Control for Modular Multilevel Converter. IEEE Trans. Ind. Electron. 2016, 63, 716–727. [Google Scholar] [CrossRef]
  11. Ilves, K.; Norrga, S.; Nee, H.-P. On energy variations in modular multilevel converters with full-bridge submodules for Ac-Dc and Ac-Ac applications. In Proceedings of the 2013 15th European Conference on Power Electronics and Applications (EPE), Lille, France, 2–6 September 2013; pp. 1–10. [Google Scholar]
  12. Ilves, K.; Norrga, S.; Harnefors, L.; Nee, H.-P. On Energy Storage Requirements in Modular Multilevel Converters. IEEE Trans. Power Electron. 2014, 29, 77–88. [Google Scholar] [CrossRef]
  13. Song, Q.; Yang, W.; Zhao, B.; Xu, S.; Rao, H.; Zhu, Z. Energy Storage Requirement Reduction Using Negative-Voltage States of a Full-Bridge Modular Multilevel Converter. IEEE Trans. Power Electron. 2019, 34, 5243–5255. [Google Scholar] [CrossRef]
  14. Judge, P.D.; Chaffey, G.; Merlin, M.M.C.; Clemow, P.R.; Green, T.C. Dimensioning and Modulation Index Selection for the Hybrid Modular Multilevel Converter. IEEE Trans. Power Electron. 2017, 33, 3837–3851. [Google Scholar] [CrossRef] [Green Version]
  15. Dong, P.; Lyu, J.; Cai, X. Optimized Design and Control for Hybrid MMC With Reduced Capacitance Requirements. IEEE Access 2018, 6, 51069–51083. [Google Scholar] [CrossRef]
  16. Zhao, C.; Lei, M.; Hu, Y.; Li, Z.; Gao, F.; Wang, P.; Li, Y. Energy Storage Requirement Optimization of Hybrid Modular Multilevel Converter with Circulating Current Injection. IEEE Trans. Ind. Electron. 2019, 66, 6637–6648. [Google Scholar] [CrossRef]
  17. Lin, L.; Lin, Y.; Xu, C.; Chen, Y. Comprehensive Analysis of Capacitor Voltage Fluctuation and Capacitance Design for Submodules in Hybrid Modular Multilevel Converter with Boosted Modulation Index. IEEE J. Emerg. Sel. Top. Power Electron. 2018, 7, 2369–2383. [Google Scholar] [CrossRef]
  18. Hu, J.; Xiang, M.; Lin, L.; Lu, M.; Zhu, J.; He, Z. Improved Design and Control of FBSM MMC with Boosted AC Voltage and Reduced DC Capacitance. IEEE Trans. Ind. Electron. 2018, 65, 1919–1930. [Google Scholar] [CrossRef]
  19. Zeng, R.; Xu, L.; Yao, L.; Williams, B.W. Design and Operation of a Hybrid Modular Multilevel Converter. IEEE Trans. Power Electron. 2015, 30, 1137–1146. [Google Scholar] [CrossRef] [Green Version]
  20. Lu, M.; Hu, J.; Zeng, R.; Li, W.; Lin, L. Imbalance Mechanism and Balanced Control of Capacitor Voltage for a Hybrid Modular Multilevel Converter. IEEE Trans. Power Electron. 2018, 33, 5686–5696. [Google Scholar] [CrossRef]
  21. Lu, M.; Hu, J.; Zeng, R.; He, Z. Fundamental-Frequency Reactive Circulating Current Injection for Capacitor Voltage Balancing in Hybrid-MMC HVDC Systems During Riding Through PTG Faults. IEEE Trans. Power Deliv. 2018, 33, 1348–1357. [Google Scholar] [CrossRef]
  22. Lee, J.-H.; Jung, J.-J.; Sul, S.-K. Balancing of Submodule Capacitor Voltage of Hybrid Modular Multilevel Converter under DC-Bus Voltage Variation of HVDC System. IEEE Trans. Power Electron. 2019, 34, 10458–10470. [Google Scholar] [CrossRef]
  23. Sidorov, D.; Panasetsky, D.; Tomin, N.; Karamov, D.; Zhukov, A.; Muftahov, I.; Dreglea, A.; Liu, F.; Li, Y. Toward zero-emission hybrid AC/DC power systems with renewable energy sources and storages: A case study from lake Baikal region. Energies 2020, 13, 1226. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Illustration of the hybrid MMC.
Figure 1. Illustration of the hybrid MMC.
Energies 15 01225 g001
Figure 2. Equivalent circuit of the MMC connected to the AC grid.
Figure 2. Equivalent circuit of the MMC connected to the AC grid.
Energies 15 01225 g002
Figure 3. Relation between the number of FBSMs and the arm voltage: (a) HB-MMC; (b) hybrid MMC.
Figure 3. Relation between the number of FBSMs and the arm voltage: (a) HB-MMC; (b) hybrid MMC.
Energies 15 01225 g003
Figure 4. P/Q operating region of MMC considering maximum reactive power requirement.
Figure 4. P/Q operating region of MMC considering maximum reactive power requirement.
Energies 15 01225 g004
Figure 5. Flowchart of the sorting algorithm.
Figure 5. Flowchart of the sorting algorithm.
Energies 15 01225 g005
Figure 6. Illustration of capacitor voltage variation process: (a) capacitor voltages; (b) arm voltage; (c) arm current; (d) input submodule number.
Figure 6. Illustration of capacitor voltage variation process: (a) capacitor voltages; (b) arm voltage; (c) arm current; (d) input submodule number.
Energies 15 01225 g006
Figure 7. Computing process of capacitor voltages.
Figure 7. Computing process of capacitor voltages.
Energies 15 01225 g007
Figure 8. Computing process of capacitor usage and FBSM quantity.
Figure 8. Computing process of capacitor usage and FBSM quantity.
Energies 15 01225 g008
Figure 9. Schematic of the main circuit of the case study.
Figure 9. Schematic of the main circuit of the case study.
Energies 15 01225 g009
Figure 10. Simulation and calculation waveform of capacitor voltages in different operating points: (a) pure active output (inverting) (φ = 0); (b) pure active input (rectifying) (φ = π); (c) pure inductive reactive power (φ = −π/2); (d) pure capacitive reactive power (φ = π/2).
Figure 10. Simulation and calculation waveform of capacitor voltages in different operating points: (a) pure active output (inverting) (φ = 0); (b) pure active input (rectifying) (φ = π); (c) pure inductive reactive power (φ = −π/2); (d) pure capacitive reactive power (φ = π/2).
Energies 15 01225 g010
Figure 11. Capacitor usage of hybrid MMC varying with M0 and Qmax.
Figure 11. Capacitor usage of hybrid MMC varying with M0 and Qmax.
Energies 15 01225 g011
Table 1. Relation between arm voltages of FBSMs and HBSMs.
Table 1. Relation between arm voltages of FBSMs and HBSMs.
u arm < 0   u arm > 0
i arm > 0 i arm < 0
Only FBSM input:
u f = u arm
u h = 0
u c f > u c h Maximum HBSM input:
u h = min ( u arm ,   N 0 U c )
u f = u arm u h
Maximum FBSM input.
u c f = u c h Equally input:
u f = E nomf E nom u arm
u h = E nomh E nom u arm
Equally input.
u c f < u c h Maximum FBSM input:
u f = min ( u arm , F U c )
u h = u arm u f
Maximum HBSM input.
Table 2. Main parameters of the case study.
Table 2. Main parameters of the case study.
ParametersValues
Rated capacity1250 MVA
Rated DC voltage400 kV
Rated AC voltage (line)294 kV (M0 = 1.2)
Maximum reactive power1 p.u.
Rated capacitor DC voltage2 kV
Equivalent AC reactance44 mH (0.25 p.u.)
HBSM quantity200 (1 p.u.)
FBSM quantity50 (0.25 p.u.)
Capacitor usage35.7 kJ/MVA
Capacitance ratio1.3
Capacitance of FB/HBSM18.2/14 mF
Table 3. Comparison table of different methods.
Table 3. Comparison table of different methods.
AccuracyCompletenessLow Complexity
Arm power calculation [13,16]
Peak voltage balancing [14]
Piecewise calculation [19,20,21,22]
Detailed piecewise calculation [17]
Proposed method
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, K.; Song, Q.; Xu, S. Analysis and Design of the Energy Storage Requirement of Hybrid Modular Multilevel Converters Using Numerical Integration and Iterative Solution. Energies 2022, 15, 1225. https://doi.org/10.3390/en15031225

AMA Style

Wang K, Song Q, Xu S. Analysis and Design of the Energy Storage Requirement of Hybrid Modular Multilevel Converters Using Numerical Integration and Iterative Solution. Energies. 2022; 15(3):1225. https://doi.org/10.3390/en15031225

Chicago/Turabian Style

Wang, Kailun, Qiang Song, and Shukai Xu. 2022. "Analysis and Design of the Energy Storage Requirement of Hybrid Modular Multilevel Converters Using Numerical Integration and Iterative Solution" Energies 15, no. 3: 1225. https://doi.org/10.3390/en15031225

APA Style

Wang, K., Song, Q., & Xu, S. (2022). Analysis and Design of the Energy Storage Requirement of Hybrid Modular Multilevel Converters Using Numerical Integration and Iterative Solution. Energies, 15(3), 1225. https://doi.org/10.3390/en15031225

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