Reduced Order Generalized Integrator Based Modular Multilevel Converter Loop Current Suppression Strategy under Unbalanced Conditions in Distribution Networks

: Under the condition of grid voltage imbalance, the circulation of the bridge arm inside the modular multilevel converter (MMC) increases significantly, which leads to the aggravation of the distortion of the bridge arm current, and, thus, increases the system loss and reduces the power quality. To address this problem, this paper analyzes the mechanism of circulating current generation and proposes a circulating current suppression strategy based on a reduced-order generalized integrator (ROGI), which firstly uses the ROGI system to separate the second-harmonic positive-and negative-sequence components in the circulating current from the DC, and then converts the rotating coordinates of the circulating current’s second octave component into the DC to be fed into the proportional–integral quasi-resonance (PIR) controller for suppression. A simulation model of a 23-level MMC inverter is built in MATLAB/Simulink, and the control strategy proposed in this paper is compared with the classical proportional–integral (PI) control in simulation experiments. The simulation results show that the amplitude of the circulating current fluctuation of the classical PI control is reduced from 90 A to 22 A, and the harmonic distortion rate of the bridge arm current is reduced from 32.56% to 5.57%; the amplitude of the circulating current fluctuation of the control strategy proposed in this paper is reduced from 90 A to 5.7 A, and the harmonic distortion rate of the bridge arm current is reduced from 20.2% to 1.13%, which verifies the effectiveness of the pro-posed control strategy.


Introduction
Due to the penetration of new renewable energy sources and the economics of highvoltage DC transmission, high-voltage DC transmission systems are becoming the primary means of handling and transmitting high-power electricity [1][2][3][4].Currently, the MMC is considered the most ideal voltage source converter solution for HVDC systems due to its modular and scalable architecture.
When unbalanced conditions occur in the power grid, the MMC, due to its unique distributed arrangement of energy storage capacitors, is prone to causing imbalance in the capacitor voltage of each submodule.The energy imbalance between the phases leads to the generation of circulating current inside the converter, which adds to the current distortion of the bridge arm, increases the energy loss of the bridge arm, affects the stable operation and service life of the submodule components, and increases the system cost.Therefore, it is necessary to carry out suppression measures for the circulating current [5,6].
To suppress the internal circulating current within the MMC, extensive research has been conducted by scholars both domestically and internationally.Currently, the adopted circulating current suppression strategies are primarily classified into hardware-based and software-based methods [7].Hardware methods mainly involve increasing the arm inductance or employing auxiliary circuits to suppress the circulating current.However, such approaches tend to increase the system's size and losses, leading to adverse effects [8].
In reference [9], focusing on the unique submodule structure of the MMC in high-voltage battery energy storage systems, it is proposed that increasing the arm inductance value under certain threshold conditions can help suppress the interphase circulating current of the system.Nonetheless, this method requires substantial calculations and complex design.Additionally, reference [10] suggests using a zero-sequence circulating current controller to effectively increase the equivalent arm resistance, thereby enhancing system damping and effectively suppressing the impedance resonance peak of the MMC under islanding control mode.Nevertheless, this method's control complexity poses inconveniences for practical applications.
The software-based approach primarily focuses on suppressing the harmonic components of the secondary current through designing an appropriate controller, thereby inhibiting the circulating current.This method, known for its low cost and ease of debugging, has become the mainstream choice for circulating current suppression [11].However, it still suffers from drawbacks such as slow dynamic response and sensitivity to controller parameters.In reference [12], the abc-dq coordinate transformation is proposed to decompose the three-phase circulating current into direct current components, upon which a circulating current suppressor is designed.Nevertheless, this method requires coordinate transformation and interphase decoupling, limiting the system's stability and circulating current suppression capability.Reference [13] introduces a method that employs coordinate transformation to decompose the three-phase circulating current into dc components in the dq coordinate system, which are then suppressed using a PI controller.This method boasts a simple structure and convenient operation but is prone to external influences on its control parameters, resulting in limited control accuracy.References [14,15] present a composite controller combining repetitive control with a proportional element, capable of simultaneously eliminating multiple harmonic circulating currents in power systems and exhibiting excellent tracking performance.Subsequently, reference [16] incorporates a fractional-order phase lead compensator into the aforementioned structure to further enhance dynamic performance, albeit at the cost of increased computational complexity.
In reference [17], after separating the positive-and negative-sequence components of the circulating current, synchronous rotating coordinate transformation control is adopted.However, the negative-sequence component separation process is overly complex, requiring multiple low-pass filters, which slows down the control speed.On the other hand, reference [18] transforms the three-phase circulating current into the αβ0 coordinate system and utilizes a proportional resonance (PR) controller to suppress the circulating current in the stationary coordinate system.Nevertheless, the control effect is not ideal when grid deviations occur.In reference [19], a PR controller-based circulating current suppression method is proposed for both single-phase and three-phase systems.However, this approach cannot effectively handle situations involving multiple harmonics simultaneously, and the PR controller suffers from a narrow bandwidth.To address the bandwidth limitation of traditional PR controllers, references [20][21][22] propose a QPR controller-based circulating current suppression strategy, which enhances system stability and response speed.Nevertheless, this method is only applicable to two-terminal MMC transmission systems and single-terminal grid voltage imbalance scenarios.Reference [23] introduces a circulating current suppression strategy utilizing PR and PI controllers in the abc reference frame, which not only reduces circulating currents but also minimizes the fluctuation amplitude of capacitor voltages.However, this method involves complex calculations.On the other hand, reference [24] presents a simpler approach using two PI controllers to suppress circulating currents, offering a straightforward and convenient solution but with limited control accuracy.In reference [25], a direct power control strategy for MMC-HVDC based on a reduced-order vector resonator is proposed.This strategy employs PI controllers to regulate the dc components of instantaneous power, while the reduced-order vector resonator is utilized to eliminate the second-order positive-sequence fluctuating components.This approach enables direct control of active and reactive power.However, the control strategy is complex and involves numerous control parameters.Reference [26] presents a circulating current suppressor based on PIR control.However, the system exhibits relatively weak stability and is therefore not suitable for unbalanced systems.
References [27,28] propose circulating current suppression strategies based on the second order generalized integrator (SOGI), where a QPR controller is designed to suppress the second-order harmonics, thereby reducing the computational burden of the system.Additionally, reference [29] introduces a circulating current suppression strategy based on the ROGI, which further minimizes the overall computational load of the system.However, both methods are discussed and analyzed under the assumption of balanced grid conditions and are, thus, not applicable to unbalanced grid systems.
Addressing the aforementioned issues, this paper first establishes the topological structure and mathematical analysis of circulating current in MMCs, elucidating the mechanism of circulating current generation.Subsequently, the second-order positive-and negative-sequence components of the circulating current are separated using a reducedorder generalized integrator.These positive-, negative-, and zero-sequence second-order components are then transformed into appropriate rotating coordinates and suppressed by a PIR controller.This method enables independent control of the positive-, negative-, and zero-sequence second-order components in the bridge arms without requiring auxiliary circuits or modifications to the topology.While preserving the impedance characteristics of the MMC, it stabilizes submodule capacitor voltage fluctuations.The system features a simple structure, fast response, and robust stability.

Topology and Working Principle of Active Distribution Network MMCs
The three-phase MMC inverter topology is shown in Figure 1.The MMC consists of six bridge arms in three phases, with each phase divided into two bridge arms.Each bridge arm consists of an inductor and several half-bridge submodules with the same structure.Each half-bridge submodule consists of two IGBTs and a set of capacitors.L 0 , R 0 are the equivalent reactance and equivalent loss, respectively; i pj , i nj (j = a,b,c) are the currents of the upper and lower bridge arms of phase j, respectively, and u sj is the AC side voltage of phase j.To simplify the analysis, the reactor and the ideal source are considered as a whole and the whole output is considered as an ideal source to the outside world.From Kirchhoff's law: To simplify the analysis, the reactor and the ideal source are considered as a whole and the whole output is considered as an ideal source to the outside world.From Kirchhoff's law: Adding and subtracting Equations ( 2) and (3), respectively, gives the AC and DC side characteristics of the MMC as follows: where U di f f j , i cirj and u comj denote the differential mode voltage, circulating current, and common mode voltage, respectively, and are defined as where U di f f j and U comj are regulated by the inner circulation controller and i cirj is obtained by the circulation suppression controller.

Circulation Analysis within the MMC under Unbalanced Conditions
Under normal conditions, the MMC internal circulation can usually be expressed by Equation (9): where θ is the second-harmonic initial phase of the ring current, the high-frequency harmonic components in the ring current are generated by the twofold frequency layer-by-layer excitation of the AC quantity, and the harmonic amplitude content decreases as the frequency increases, where Q 1 is the sum of the harmonics above the second harmonic, which is generally ignored due to the small content.
In the case of unbalanced grid voltage, the MMC phase-to-phase circulating currents are altered, generating positive-sequence, negative-sequence, and zero-sequence components of twofold frequency.The AC-side converter in the MMC structure usually adopts the Y-∆ connection, and there is no circulating loop for the zero-sequence current.Neglecting the zero-sequence component effect U acj , I sj can be decomposed into positive-and negative-sequence components, expressed as where ω 0 is the angular frequency of the grid, U p , U n are the positive-and negativesequence components of the AC side voltage, θ p , θ n are the positive-and negative-sequence voltage of the initial phase angle, I p , I n are the positive-and negative-sequence components of the AC side current, and φ p , φ n are the positive-and negative-sequence current of the initial phase angle.
In the three-phase asymmetrical operation, redefining the currents in the upper and lower bridge arms have (11) Energies 2024, 17, 4270 5 of 16 Simplification shows that the instantaneous power of the phase cell is expressed as p j = u pj i pj + u nj i nj = U acj I sj + U dc I dcj (12) Substituting the above Equation ( 12) and analyzing it with phase a as an example, we obtain P a = P a0 + P p a2 + P n a2 + P z a2 (13) where P a0 is the dc component, and P p a2 , P n a2 , P z a2 are the twofold positive-, negative-, and zero-order components, respectively, expressed as From the analysis of the above equations, it can be seen that the negative-sequence current and voltage produce positive-sequence two-times frequency loop current; positivesequence current and voltage produce negative-sequence two-times frequency loop current; positive-sequence voltage and negative-sequence current and negative-sequence voltage and positive-sequence current produce zero-sequence two-times frequency loop current.
Integrating Equation ( 13) yields the instantaneous energy of phase j: where W j_dc is the DC component of the instantaneous energy, W p j2 , W n j2 , and W z j2 are the twofold positive-, negative-, and zero-order components of the instantaneous energy, respectively, and p j0 contains the DC component and the fundamental frequency component, defining the fundamental frequency component to be p ′ j0 , expressed as Through the above analysis of the energy, it can be seen that harmonics will be generated in the bridge arm circulating current, the fundamental frequency voltage and current will generate positive-, negative-, and zero-sequence two times the frequency of the bridge arm voltage and current, and it can be deduced that the circulating current under the unbalanced condition of the grid is defined as shown in Equation (17).The instantaneous energy equivalent circuit diagram of the phase unit of the MMC in the case of the unbalanced grid is shown in Figure 2.
Energies 2024, 17, 4270 6 of 16 the unbalanced condition of the grid is defined as shown in Equation ( 17).The instantaneous energy equivalent circuit diagram of the phase unit of the MMC in the case of the unbalanced grid is shown in Figure 2.
The twofold component in the loop current originates from the voltage difference between the three-phase bridge arms.Although it does not directly affect the normal operation of the MMC, it will increase the burden of the capacitance of the bridge arm submodules and trigger additional energy loss.Therefore, in order to improve the system efficiency, the second-harmonic component in the MMC loop current must be effectively suppressed.

Phase Sequence Separation
Among the methods for extracting the diophantine components, various techniques have their advantages and disadvantages.The trap method is slow to respond and may not be suitable for real-time applications.The delay method is memory-intensive and sensitive to signal variations and noise.The second-order generalized integrator is widely used due to its simple structure and its ability to suppress harmonics.The transfer function of SOGI is where r  is the resonant angular frequency.As can be seen from Equation ( 18), the poles of the system are ± r j , indicating that the SOGI exhibits frequency selectivity but lacks polarity selectivity, making it unable to directly separate positive-and negative-sequence fundamental components.Therefore, by performing a reduction process on the SOGI, two reduced-order generalized integrators with the same resonant frequency but opposite phase sequences can be obtained.The transfer function of this ROGI is () From this transfer function, it can be seen that there is only one pole, so it can only control the component whose angular frequency is r  , and when the angular frequency The twofold component in the loop current originates from the voltage difference between the three-phase bridge arms.Although it does not directly affect the normal operation of the MMC, it will increase the burden of the capacitance of the bridge arm submodules and trigger additional energy loss.Therefore, in order to improve the system efficiency, the second-harmonic component in the MMC loop current must be effectively suppressed.

Phase Sequence Separation
Among the methods for extracting the diophantine components, various techniques have their advantages and disadvantages.The trap method is slow to respond and may not be suitable for real-time applications.The delay method is memory-intensive and sensitive to signal variations and noise.The second-order generalized integrator is widely used due to its simple structure and its ability to suppress harmonics.The transfer function of SOGI is where ω r is the resonant angular frequency.As can be seen from Equation (18), the poles of the system are ±jω r , indicating that the SOGI exhibits frequency selectivity but lacks polarity selectivity, making it unable to directly separate positive-and negative-sequence fundamental components.Therefore, by performing a reduction process on the SOGI, two reduced-order generalized integrators with the same resonant frequency but opposite phase sequences can be obtained.The transfer function of this ROGI is From this transfer function, it can be seen that there is only one pole, so it can only control the component whose angular frequency is ω r , and when the angular frequency is equal to the resonant angular frequency, it has infinite gain, so it can be used to directly extract the positive-and negative-sequence components of the power grid.Furthermore, since at the same resonant frequency, the ROGI can provide different gains for positiveand negative-sequence signals, while the SOGI can only provide the same gain, the total computational load when using ROGI for signal modulation is approximately half that of using SOGI, thereby enhancing the computational efficiency of the control process.
In separating the positive-and negative-sequence components of harmonics, the presence of the complex number j in the ROGI transfer function is more detrimental to the mathematical realization, so that j can be realized in a three-phase system using the αβ Energies 2024, 17, 4270 7 of 16 coordinate system variable i α = ji β .At this time, the structural block diagram of ROGI is shown in Figure 3, and the transfer function is where k is the system gain.
computational load when using ROGI for signal modulation is approximately half that of using SOGI, thereby enhancing the computational efficiency of the control process.
In separating the positive-and negative-sequence components of harmonics, the presence of the complex number j in the ROGI transfer function is more detrimental to the mathematical realization, so that j can be realized in a three-phase system using the αβ coordinate system variable j ii  = .At this time, the structural block diagram of ROGI is shown in Figure 3, and the transfer function is where k is the system gain.
When the grid is unbalanced, the MMC contains a two-times frequency component internally, which is derived from Equation ( 5), and the dynamic equations are obtained by Clark transformation: From Equations ( 19) and (20), it can be seen that the output of D(s) can be regarded as a combination of a bandpass filter and a low-pass filter, and the positive-and negativesequence current components generated when the grid is unbalanced can be quickly separated by ROGI, but the existence of the DC component will bring about a large error, so the regulator adds a DC compensator to solve the problem of ROGI's inability to filter out the DC component by injecting the DC component.The principle block diagram of positive-and negative-sequence separation is shown in Figure 4.When the grid is unbalanced, the MMC contains a two-times frequency component internally, which is derived from Equation ( 5), and the dynamic equations are obtained by Clark transformation: From Equations ( 19) and ( 20), it can be seen that the output of D(s) can be regarded as a combination of a bandpass filter and a low-pass filter, and the positive-and negativesequence current components generated when the grid is unbalanced can be quickly separated by ROGI, but the existence of the DC component will bring about a large error, so the regulator adds a DC compensator to solve the problem of ROGI's inability to filter out the DC component by injecting the DC component.The principle block diagram of positive-and negative-sequence separation is shown in Figure 4.

Circulation Suppression Strategies
Due to the asymmetry of three-phase loads and asymmetrical faults in the power system, the grid is prone to unbalanced conditions, which leads to an increase in the harmonic content of the two-frequency loop current in the grid-connected current.As shown

Circulation Suppression Strategies
Due to the asymmetry of three-phase loads and asymmetrical faults in the power system, the grid is prone to unbalanced conditions, which leads to an increase in the harmonic content of the two-frequency loop current in the grid-connected current.As shown in the previous section, the loop current i ′ cirj is the superposition of the DC component and the harmonic component.The PI controller can realize the static-free tracking of the DC input, which has the advantages of high reliability and robustness, but this regulator cannot realize the static-free tracking of the sinusoidal input.The ideal PR control has a relatively narrow bandwidth frequency, with a small increase at the resonant frequency point as the frequency varies.This results in a significant drop in open-loop gain when the system frequency changes.To broaden the system bandwidth, enhance stability, and facilitate digital implementation, this paper adopts a combination of the QPR controller and the traditional PI controller, forming the PIR controller, whose structure block diagram is shown in Figure 5, and the transfer function is expressed as Equation ( 23): where K p , K i are the proportional and integral coefficients, ω c is the cutoff frequency, and K r is the twofold resonant controller gain.
Energies 2024, 17, 4270 10 of 18 The circulating current within the bridge arm flows back and forth within the three phases, and a twofold negative/positive-sequence rotating coordinate transformation is used to convert the negative/positive-sequence twofold circulating current into a direct flow to analyze the negative-sequence component as an example there: Based on Equation ( 24), the three-phase time-varying circulating currents can be decomposed into two DC components through the negative-sequence rotating coordinate transformation.This decomposition facilitates the design of a circulating current suppressor using PI controllers.The remaining double-frequency components of the zero-sequence circulating current will flow through the DC lines, affecting the normal operation of other converter stations, and, therefore, must also be suppressed by the QPR controller.Based on the above control approach, the control block diagram for the internal arm circulating current of MMC under unbalanced grid conditions is obtained, as shown in Figure 6.The principle block diagram of this controller is shown in Figure 5.The principle of the current loop PI controller in parallel with the QPR controller is to reduce the THD of the bridge arm current, in which the QPR controller at the fundamental frequency of the gain is ∞.When the command value of the given current and the feedback current there are harmonics, the PIR controller amplifies the error through the negative feedback link to inhibit harmonic components of the current.
The circulating current within the bridge arm flows back and forth within the three phases, and a twofold negative/positive-sequence rotating coordinate transformation is used to convert the negative/positive-sequence twofold circulating current into a direct flow to analyze the negative-sequence component as an example there: Based on Equation ( 24), the three-phase time-varying circulating currents can be decomposed into two DC components through the negative-sequence rotating coordinate transformation.This decomposition facilitates the design of a circulating current suppressor using PI controllers.The remaining double-frequency components of the zero-sequence circulating current will flow through the DC lines, affecting the normal operation of other converter stations, and, therefore, must also be suppressed by the QPR controller.Based on the above control approach, the control block diagram for the internal arm circulating current of MMC under unbalanced grid conditions is obtained, as shown in Figure 6.
transformation.This decomposition facilitates the design of a circulating current suppressor using PI controllers.The remaining double-frequency components of the zero-sequence circulating current will flow through the DC lines, affecting the normal operation of other converter stations, and, therefore, must also be suppressed by the QPR controller.Based on the above control approach, the control block diagram for the internal arm circulating current of MMC under unbalanced grid conditions is obtained, as shown in Figure 6.

PIR Controller Parameterization
The setting of parameters for the PIR controller primarily involves adjusting the parameters of the proportional-integral controller and the quasi-resonant controller.By altering the values of various parameters within the controller, the control effect can be fine-tuned.
K p determines the baseline gain across the entire frequency range.An increase in the K p value enhances the suppression capability of the resonant controller for harmonics near the resonant frequency.However, if the K p value is too high, it may compromise system stability, leading to oscillations.K i primarily affects the gain in the low-frequency range and is used to eliminate steady-state errors in the system.Nevertheless, excessively high values of K p and K i can cause phase lag and stability issues.The frequency domain analysis results are shown in Figure 7a,b.
K r determines the resonant gain at the resonant frequency.As K r increases, the gain at the tracked frequency point improves, but an excessively high K r value can lead to phase advance.ω c dictates the bandwidth of the controller; as ω c increases, so does the controller's bandwidth.However, an excessively large ω c can increase interference to nonharmonic frequencies.The frequency domain analysis results are presented in Figure 7c,d.
In summary, the optimal parameters of K p , K i , K r and ω c are determined based on actual conditions, with K p = 50, K i = 20, K r = 400, and ω c = 5 rad/s.Figure 8 shows the zero-pole diagram of the circulating current suppressor with these parameters, where represents zeros and × represents poles.All characteristic roots of the system have negative real parts, indicating that the controller is stable.
ics near the resonant frequency.However, if the p K value is too high, it may compromise system stability, leading to oscillations.○ represents zeros and × represents poles.All characteristic roots of the system have negative real parts, indicating that the controller is stable.In this paper, the circulation suppression strategy first measures the bridge arm current of the MMC, calculates the bridge arm loop current value ' cirj i , extracts the AC two- times frequency component in the loop current by the ROGI harmonic extractor, and then inputs it to the PIR loop current controller to obtain the voltage reference value for eliminating the harmonic component diffj u .The circulation suppression strategy in this paper is shown in Figure 9.In this paper, the circulation suppression strategy first measures the bridge arm current of the MMC, calculates the bridge arm loop current value i ′ cirj , extracts the AC two-times frequency component in the loop current by the ROGI harmonic extractor, and then inputs it to the PIR loop current controller to obtain the voltage reference value for eliminating the harmonic component u di f f j .The circulation suppression strategy in this paper is shown in Figure 9.
rent of the MMC, calculates the bridge arm loop current value ' cirj i , extracts the AC two- times frequency component in the loop current by the ROGI harmonic extractor, and then inputs it to the PIR loop current controller to obtain the voltage reference value for eliminating the harmonic component diffj u .The circulation suppression strategy in this paper is shown in Figure 9.

Modulation Strategy
In the conventional nearest level approximation modulation strategy, the actual output voltage between the upper and lower bridge arms deviates from the reference voltage, which makes the actual AC energy provided by the AC side unequal to that of the submodule capacitor voltage fluctuation, thus generating a double-frequency loop current.When the actual output characteristics of the bridge arm voltages are not affected by the fluctuations of the submodule capacitor voltages, there is no twofold frequency circulating current between the phase units, and the power quality can be effectively improved.
In this paper, an approximate superposition modulation strategy is chosen: This strategy considers that the submodule capacitance voltage is changing all the time, and through the superposition of submodule capacitance voltages, the actual values of the input submodule voltages of the upper and lower bridge arms are always close to the reference value of the modulating waveform.Considering that when the actual submodule capacitance voltages are superimposed, the sum of the final input submodule capacitance voltages cannot be exactly the same as the target level; it is chosen to take the absolute value of the deviation from the target level as the criterion, so as to achieve the purpose that the sum of the input submodule capacitance voltages is closest to the target level.The specific steps of this control strategy are as follows: 1. Boost sequencing of the submodule capacitor voltages for each bridge arm.2. When the current charges the capacitor voltage of the submodule (current greater than 0), the capacitor voltages of the submodules are summed from front to back until the sum of the capacitor voltages of the submodules is closest to the given value,

Modulation Strategy
In the conventional nearest level approximation modulation strategy, the actual output voltage between the upper and lower bridge arms deviates from the reference voltage, which makes the actual AC energy provided by the AC side unequal to that of the submodule capacitor voltage fluctuation, thus generating a double-frequency loop current.When the actual output characteristics of the bridge arm voltages are not affected by the fluctuations of the submodule capacitor voltages, there is no twofold frequency circulating current between the phase units, and the power quality can be effectively improved.
In this paper, an approximate superposition modulation strategy is chosen: This strategy considers that the submodule capacitance voltage is changing all the time, and through the superposition of submodule capacitance voltages, the actual values of the input submodule voltages of the upper and lower bridge arms are always close to the reference value of the modulating waveform.Considering that when the actual submodule capacitance voltages are superimposed, the sum of the final input submodule capacitance voltages cannot be exactly the same as the target level; it is chosen to take the absolute value of the deviation from the target level as the criterion, so as to achieve the purpose that the sum of the input submodule capacitance voltages is closest to the target level.The specific steps of this control strategy are as follows: 1.
Boost sequencing of the submodule capacitor voltages for each bridge arm.

2.
When the current charges the capacitor voltage of the submodule (current greater than 0), the capacitor voltages of the submodules are summed from front to back until the sum of the capacitor voltages of the submodules is closest to the given value, thereby determining the number of submodules that should be engaged, i.e., engaging the first S submodules after sorting.

3.
When the current discharges the submodule capacitor voltage (with the current being less than 0), the submodule capacitor voltages are summed up from the latter to the former until the sum of the submodule capacitor voltages is the closest to the given value, thereby determining the number of submodules to be activated, i.e., the last N−S+1 submodules after sorting.
This control strategy casts cuts in each cycle strictly according to the submodule voltage magnitude, which can limit the submodule voltage unbalance to a fairly small range.The specific flow of the modulation strategy is shown in Figure 10: the former until the sum of the submodule capacitor voltages is the closest to the given value, thereby determining the number of submodules to be activated, i.e., the last N−S+1 submodules after sorting.This control strategy casts cuts in each cycle strictly according to the submodule voltage magnitude, which can limit the submodule voltage unbalance to a fairly small range.The specific flow of the modulation strategy is shown in Figure 10

Simulation and Results Analysis
In order to verify the effectiveness of the circulating current suppression method mentioned in this paper, a 23-level MMC simulation model is built in MATLAB/Simulink simulation software; the modulation method adopts the nearest level approximation modulation, and the model simulation parameters are shown in Table 1.In order to reflect the control effect of the circulating current suppression strategy proposed in this paper under the unbalanced condition of the grid, it is set that the grid voltage fluctuation occurs 0.3 s after the balanced state, the amplitude of the A-phase voltage falls to 60%, and the circulating current suppressor is activated at the moment of 0.5 s.

Simulation and Results Analysis
In order to verify the effectiveness of the circulating current suppression method mentioned in this paper, a 23-level MMC simulation model is built in MATLAB/Simulink simulation software; the modulation method adopts the nearest level approximation modulation, and the model simulation parameters are shown in Table 1.In order to reflect the control effect of the circulating current suppression strategy proposed in this paper under the unbalanced condition of the grid, it is set that the grid voltage fluctuation occurs 0.3 s after the balanced state, the amplitude of the A-phase voltage falls to 60%, and the circulating current suppressor is activated at the moment of 0.5 s.Under the unbalanced condition, the three-phase output voltage on the AC side is shown in Figure 11a, which is changed at t = 0.3 s and stabilized the rest of the time.The FFT analysis of this output voltage is shown in Figure 11b, and the THD of the three-phase output voltage is 2.94%, so it can be seen that the distortion rate of this output voltage is small and still available.
Under the unbalanced condition, the three-phase output voltage on the AC s shown in Figure 11a, which is changed at t = 0.3 s and stabilized the rest of the tim FFT analysis of this output voltage is shown in Figure 11b, and the THD of the threeoutput voltage is 2.94%, so it can be seen that the distortion rate of this output volt small and still available.Figure 12 shows the changes of the submodule capacitor voltage waveforms during the simulation process.The capacitor voltage waveforms are in good condition in steady state, the capacitor voltage distortion is serious after the fluctuation of the grid voltage in 0.3 s, and the distortion is obviously weakened after turning on the circulating current suppressor in 0.5 s.The fluctuation amplitude of the submodule capacitor voltages is also reduced, which is in good condition.Figures 13 and 14, respectively, depict the current waveform and the circulating current waveform of the a-phase bridge arm.After the voltage of the a-phase grid drops (at t = 0.3 s), the current waveform of the a-phase bridge arm experiences severe distortion, with the second-harmonic amplitude increasing from 49 A to 90 A and the harmonic distortion rate rising from 19.69% to 20.2%, adversely affecting the stable operation of the system.After 0.5 s of operation, the circulation suppressor designed in this paper is activated, resulting in a reduction in the distortion of the a-phase bridge arm current waveform.Specifically, the second-harmonic amplitude of the circulating current decreases from 90 A to 5.7 A, and the harmonic distortion rate drops to 1.13%, significantly improving the waveform quality and transforming it into a more ideal sinusoidal wave.Figures 13 and 14, respectively, depict the current waveform and the circulating current waveform of the a-phase bridge arm.After the voltage of the a-phase grid drops (at t = 0.3 s), the current waveform of the a-phase bridge arm experiences severe distortion, with the second-harmonic amplitude increasing from 49 A to 90 A and the harmonic distortion rate rising from 19.69% to 20.2%, adversely affecting the stable operation of the system.After 0.5 s of operation, the circulation suppressor designed in this paper is activated, resulting in a reduction in the distortion of the a-phase bridge arm current waveform.Specifically, the second-harmonic amplitude of the circulating current decreases from 90 A to 5.7 A, and the harmonic distortion rate drops to 1.13%, significantly improving the waveform quality and transforming it into a more ideal sinusoidal wave.
tortion rate rising from 19.69% to 20.2%, adversely affecting the stable operation of the system.After 0.5 s of operation, the circulation suppressor designed in this paper is activated, resulting in a reduction in the distortion of the a-phase bridge arm current waveform.Specifically, the second-harmonic amplitude of the circulating current decreases from 90 A to 5.7 A, and the harmonic distortion rate drops to 1.13%, significantly improving the waveform quality and transforming it into a more ideal sinusoidal wave.Figure 15 presents the simulation results of the circulating current suppression based on the PI control strategy.Specifically, during the grid voltage imbalance at 0.3 s, the harmonic oscillation and fluctuation of the circulating current in the a-phase bridge arm change from 48 A to 90 A, and the distortion rate increases from 19.7% to 32.56%.After activating the PI circulating current suppressor at 0.5 s, the fluctuation amplitude of the circulating current is suppressed to 22 A, and the harmonic distortion rate decreases to 5.57%.Although there is some improvement in waveform quality compared to before applying PI control, the waveform still exhibits a certain degree of distortion.Figure 15 presents the simulation results of the circulating current suppression based on the PI control strategy.Specifically, during the grid voltage imbalance at 0.3 s, the harmonic oscillation and fluctuation of the circulating current in the a-phase bridge arm change from 48 A to 90 A, and the distortion rate increases from 19.7% to 32.56%.After activating the PI circulating current suppressor at 0.5 s, the fluctuation amplitude of the circulating current is suppressed to 22 A, and the harmonic distortion rate decreases to 5.57%.Although there is some improvement in waveform quality compared to before applying PI control, the waveform still exhibits a certain degree of distortion.monic oscillation and fluctuation of the circulating current in the a-phase bridge arm change from 48 A to 90 A, and the distortion rate increases from 19.7% to 32.56%.After activating the PI circulating current suppressor at 0.5 s, the fluctuation amplitude of the circulating current is suppressed to 22 A, and the harmonic distortion rate decreases to 5.57%.Although there is some improvement in waveform quality compared to before applying PI control, the waveform still exhibits a certain degree of distortion.

Conclusions
In this paper, based on the mathematical analysis of the circulating current inside the MMC, a new circulating current suppression method is proposed for the unbalanced situation of the power grid.The following conclusions are obtained after simulation experiments by building a 23-level MMC inverter model: 1.When the grid is unbalanced, the circulating current asymmetry inside the MMC is aggravated, and the circulating current suppressor mentioned in this paper can significantly reduce the twofold frequency component in the circulating current compared with the traditional PI control, which can reduce the fluctuating range of the submodule capacitance voltage to a certain extent and improve the quality of the waveform of the bridge arm current.The structure avoids the use of traps, streamlines the calculation, and increases the stability of the system.2. When the grid is unbalanced, the circulating harmonic component inside the MMC increases significantly.There is almost no fluctuation in the system AC side voltage before and after the addition of the circulating current suppressor, which proves that the positive-and negative-sequence second-harmonic components only circulate between the internal three-phase bridge arms.

Conclusions
In this paper, based on the mathematical analysis of the circulating current inside the MMC, a new circulating current suppression method is proposed for the unbalanced situation of the power grid.The following conclusions are obtained after simulation experiments by building a 23-level MMC inverter model: When the grid is unbalanced, the circulating current asymmetry inside the MMC is aggravated, and the circulating current suppressor mentioned in this paper can significantly reduce the twofold frequency component in the circulating current compared with the traditional PI control, which can reduce the fluctuating range of the submodule capacitance voltage to a certain extent and improve the quality of the waveform of the bridge arm current.The structure avoids the use of traps, streamlines the calculation, and increases the stability of the system.2.
When the grid is unbalanced, the circulating harmonic component inside the MMC increases significantly.There is almost no fluctuation in the system AC side voltage before and after the addition of the circulating current suppressor, which proves that the positive-and negative-sequence second-harmonic components only circulate between the internal three-phase bridge arms.
Funding: This work was supported by the National Natural Science Foundation of China, grant number 52077120.

Figure 4 .
Figure 4. Block diagram of positive-and negative-sequence separation.

Figure 4 .
Figure 4. Block diagram of positive-and negative-sequence separation.

iKFigure 7 . 5 cFigure 7 .
Figure 7.The bode plot of the PIR controller with different parameters.(a) Bode plots as Kp varies; (b) Bode plots as Kp varies; (c) Bode plots as Kr varies; (d) Bode plots as wc varies.

Figure 8 .
Figure 8.The zero-pole diagram of the PIR circulating current suppressor.

Figure 8 .
Figure 8.The zero-pole diagram of the PIR circulating current suppressor.

Figure 11 .
Figure 11.Three-phase AC voltage output.(a) Three-phase AC voltage output; (b) FFT anal the three-phase AC voltage output.

Figure 12
Figure12shows the changes of the submodule capacitor voltage waveforms d the simulation process.The capacitor voltage waveforms are in good condition in s state, the capacitor voltage distortion is serious after the fluctuation of the grid volt 0.3 s, and the distortion is obviously weakened after turning on the circulating cu suppressor in 0.5 s.The fluctuation amplitude of the submodule capacitor voltages reduced, which is in good condition.

Figure 11 .
Figure 11.Three-phase AC voltage output.(a) Three-phase AC voltage output; (b) FFT analysis of the three-phase AC voltage output.

Figure 13 .
Figure 13.A-phase bridge arm current.(a) A-phase upper bridge arm current; (b) FFT analysis at t = 0.25 s; (c) FFT analysis at t = 0.65 s.

Figure 14 .
Figure 14.Simulation results of the control strategy circulation in this paper.

Figure 15 .
Figure 15.PI control strategy circulation simulation results.

Figure 15 .
Figure 15.PI control strategy circulation simulation results.

Table 1 .
Main parameters of the simulation model.