Journal of Modern Power Systems and Clean Energy

ISSN 2196-5625 CN 32-1884/TK

网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

Enhanced AC Fault Ride-through Control for MMC-integrated System Based on Active PCC Voltage Drop  PDF

  • Haihan Ye 1
  • Wu Chen 1
  • Heng Wu 2
  • Wu Cao 1
  • Guoqing He 3
  • Guanghui Li 3
1. Center for Advanced Power Conversion Technology and Equipment, School of Electrical Engineering, Southeast University, Nanjing 210096, China; 2. Department of Energy Technology, Aalborg University, 9220 Aalborg, Denmark; 3. China Electric Power Research Institute, Beijing 100192, China

Updated:2023-07-24

DOI:10.35833/MPCE.2022.000285

  • Full Text
  • Figs & Tabs
  • References
  • Authors
  • About
CITE
OUTLINE

Abstract

When a renewable energy station (RES) connects to the rectifier station (RS) of a modular multilevel converter-based high-voltage direct current (MMC-HVDC) system, the voltage at the point of common coupling (PCC) is determined by RS control methods. For example, RS control may become saturated under fault, and causes the RS to change from an equivalent voltage source to an equivalent current source, making fault analysis more complicated. In addition, the grid code of the fault ride-through (FRT) requires the RES to output current according to its terminal voltage. This changes the fault point voltage and leads to RES voltage regulation and current redistribution, resulting in fault response interactions. To address these issues, this study describes how an MMC-integrated system has five operation modes and three common characteristics under the duration of the fault. The study also reveals several instances of RS performance degradation such as AC voltage loop saturation, and shows that RS power reversal can be significantly improved. An enhanced AC FRT control method is proposed to achieve controllable PCC voltage and continuous power transmission by actively reducing the PCC voltage amplitude. The robustness of the method is theoretically proven under parameter variation and operation mode switching. Finally, the feasibility of the proposed method is verified through MATLAB/Simulink results.

I. INTRODUCTION

WITH the increase in the number of installations of large-scale renewable energy stations (RESs), their stable control and transmission are urgently demanded. A modular multilevel converter based high-voltage direct current (MMC-HVDC) can not only avoid commutation failure and reactive power compensation, but also build grid-connected voltage at the point of common coupling (PCC), thereby avoiding the need for large-capacity thermal power plant construction or AC grid expansion in remote areas. Thus, MMC-HVDC is a mainstream option for islanded renewable energy transmission.

Each RES interfaces the rectifier station (RS) of the MMC-HVDC through the PCC, which poses a high risk for an AC fault due to the large number of interconnected transmission lines. With a focus on fault analysis and fault ride-through (FRT), numerous studies have been conducted. However, they have mainly focused on grid-connected systems. References [

1]-[3] considered the existence of equilibrium points of permanent magnetic synchronous generator (PMSG) systems under low-voltage ride-through (LVRT), and eigenvalue trajectory was adopted to analyze system stability during the fault steady state. Both symmetric and asymmetric faults of modular multilevel converter (MMC) connected to a power grid were analyzed in [4], and q-axis current reduction was recommended to suppress the overcurrent. Directed arm current control was proposed to suppress DC current oscillation in [5], and the superiority of vertical balance control was demonstrated in comparison with energy-based control under an asymmetric fault in [6].

These previous studies focused on grid-connected systems, where the power grid could be Thevenin equivalent to an ideal voltage source in series with line impedance. Unlike the classic grid-connected system, the RS external response in an MMC-integrated system is determined by the RS AC voltage loop, which may become saturated under the fault, causing the RS to change from an equivalent voltage source to an equivalent current source [

7]. Under these conditions, the voltage at the fault point is affected by multiple factors that produce greater challenges for fault analysis. The AC fault of an MMC-integrated system was analyzed in [8], and DC voltage oscillation was determined to have been caused by unbalanced MMC internal energy. However, analysis was limited to the DC side, thus mostly ignoring the AC side of the FRT. The PCC overvoltage caused by the saturation of the RS AC voltage loop was reported in [9], and several strategies were then proposed to clear the integrator of the proportional-integral (PI) controller. This overvoltage can be suppressed by the anti-windup in the PI controller. This is a commonly used method but has little to do with fault response. PCC overvoltage suppression under an asymmetric fault without introducing negative sequence (NS) current was studied in [10]. However, the theoretical analysis in the study was insufficient, and numerical calculations were not conducted, thus preventing the veracity of the analysis from being verified. Unlike NS current suppression, NS voltage suppression was recommended in [11] because NS current suppression control maintains NS PCC voltage following single-phase-to-ground (P-G) fault clearance, which hinders system restoration from the fault and results in FRT failure. However, NS voltage suppression does not generate a response under a symmetric fault and thus does not additionally contribute to symmetric AC FRT. Fault current injection in an RS q-axis current reference was adopted in [12] to trigger overcurrent protection. This type of method realizes AC FRT through reasonable fault detection and circuit breaker configuration. However, considerable transmitted power is lost, and recovery after fault clearance is required, which is time-consuming. In addition, the aforementioned studies mainly focused on the solid fault, where the resistance of the fault point was 0. The validation under the solid fault only proves the validity of the proposed method under specific case, but it does not ensure its effectiveness under uncertainties.

Based on existing studies, the fault response under RES and RS control switching must be further examined, and the FRT method for an MMC-integrated system under multiple uncertainties must be improved. Accordingly, symmetric and asymmetric faults of an MMC-integrated system are fully analyzed in this study, and a high-performance FRT method is proposed, with its effectiveness demonstrated under parameter variations. The contributions of this study are described as follows.

1) The responses of symmetric and asymmetric faults under an MMC-integrated system are demonstrated. Under a symmetric fault, the MMC-integrated system has five operation modes and three common characteristics, and controllable PCC voltage and continuous power transmission are possible by optimizing PCC voltage. Under an asymmetric fault, NS voltage suppression rather than NS current suppression is adopted, and RS overcurrent and saturated modulation can be observed.

2) A high-performance AC FRT method is proposed. The proposed method can maintain continuous power transmission under fault, which is very helpful in blocking the unbalanced power that flows from the DC grid to the fault point. In addition, RS overcurrent and saturated modulation under an asymmetric fault can be effectively suppressed.

The remainder of this paper is organized as follows. Section II introduces the system description, where a basic block of the system model is presented to describe the research scope of the study. Section III analyzes the symmetric AC faults, including the RES and RS operation modes and solutions to the RS equilibrium points under fault resistance variation. Section IV analyzes asymmetric AC faults based on the findings from the Section III. Section V presents the active PCC voltage drop control, where the RS continuous power transmission area is identified and the maximum power transmission curve is derived. Moreover, the effective range calculation is conducted, and the robustness of the method is verified. Section VI presents simulation verification, where the veracity of the theoretical analysis and of the proposed method are proven through MATLAB/Simulink results. Finally, conclusions are given in Section VII.

II. System Description

The block of an RES connected to the MMC-HVDC system is shown in Fig. 1, where ur and ir are the RS terminal voltage and input current, respectively; and Rf is the fault resistance.

Fig. 1  Block of an RES connected to MMC-HVDC system.

RESs are commonly used as models of lumped controllable current sources with choppers and grid-side inverters [

13], [14], where the current references of all RESs must be determined by FRT grid code [15], [16]. An RES control algorithm is composed of a DC voltage loop, current loop, NS current suppression, and phase locked loop (PLL) based on the all-pass filter [17], [18].

The AC transmission line can be characterized by lumped resistor-inductor (RL) model [

9], [19] because the results between the aggregated parameter model and distributed line model under fault detection and analysis have proven to be very similar [20], [21]. In this study, measured data from the Zhangbei DC project of China was adopted for calculating equivalent line impedance, where 0.175 p.u. of reactive power is required when 1 p.u. of active power is transmitted.

The AC fault, marked red in Fig. 1, is generally modeled by fault resistance Rf [

22]-[24]. Both symmetric and asymmetric faults are studied, and the NS voltage suppression recommended by [11] is adopted in the RS. Analysis can be conducted in a classic two-terminal system, a multi-terminal system, or even a DC grid, as AC and DC fault characteristics are decoupled by MMC internal energy under unsaturated modulation [13], [25].

The RS adopts a bipolar configuration, where both positive and negative converters are half-bridge MMCs with identical parameters. The MMC in the RS adopts a detailed equivalent model based on double diodes [

26], and its control algorithms are composed of classic rated AC voltage control, NS voltage suppression, current loop, circulating current suppression [27], and third harmonic injection [28]. Wherein classic AC voltage control and NS voltage suppression form a centralized control, but the other control methods are distributed in both poles.

III. Symmetric AC Fault Analysis

A. Normal Operation Mode with RS Positive Power Flow

If the power loss at the fault point is less than the RES output power, the RS receives the residual power and continuously transmits it to the DC bus. In this case, the absolute value of RS current reference irdref is less than its limitation value Irmax, and thus, the d-axis PI controller of the RS AC voltage loop is not saturated, and ur can be controlled to be the rated value Urn.

The critical Rf value can be calculated by the constraint that the RS received power is greater than 0:

Pr=Ps-1.5Ur2/Rf0Ur=Urn (1)

where Ps and Pr are the active power output by the RES and that transmitted to the RS, respectively; and Ur and Urn are the measured and rated amplitudes of ur, respectively.

The main parameters of simulation are listed in Table I, where Pbase is the base of the power; Psn and udcn are the rated values of the RES output power and DC bus voltage, respectively; ωs is the rated angular frequency; ks and kr are the turn ratios of the RES- and RS-side transformers, respectively; and Lac is the equivalent collector line inductance. Figure 2 presents the simulation results when Rf=30 Ω, where usa is the a-phase voltage waveform of the equivalent RES terminal voltage us; and isq is the q-axis current of the equivalent RES terminal current is. The fault occurs at 0.6 s and clears at 1 s.

TABLE I  Main Parameters of Simulation
SymbolValueSymbolValueSymbolValue
Pbase 3000 MW Psn 2250 MW Urn 230 kV
udcn 500 kV ωs 100π ks 35/230
Lac 13.08 mH Irmax 1.133 p.u. kr 230/290

Fig. 2  Simulation results when Rf=30 Ω.

In Fig. 2, usa and isq remain stable during the fault, indicating that PCC voltage can be controlled to be the rated value, and the RES continues to operate in the rated mode. Pr decreases but is greater than 0, and therefore the continuous power transmission holds. irdref does not reach Irmax, and thus the RS voltage loop is valid in maintaining rated PCC voltage.

B. Normal Operation Mode with RS Negative Power Flow

When the RES output power is less than the power loss of the fault point, the RS active power flow is reversed. If irdref<Irmax, the PI controller of the RS AC voltage loop is not saturated, and ur is still controlled to be the rated value under the AC fault:

Ps-1.5Ur2/Rf-PrmaxUr=Urn (2)

where Prmax=1.5UrIrmax is the RS output power limitation.

Figure 3 presents the simulation results when Rf=15 Ω. The waveforms of usa and isq are similar to those shown in Fig. 2, and the RES operates under rated mode during the fault. However, Pr decreases to less than 0 during the period of 0.6-1 s. Thus, the RS power flow is reversed. In addition, irdref is increased but does not reach Irmax, which means the RS voltage loop remains effective in this case.

Fig. 3  Simulation results when Rf=15 Ω.

C. RES Normal Operation Mode Under RS Current Limitation Mode

If irdref reaches Irmax, the PI controller of the RS AC voltage loop is saturated, and RS becomes a constant current source. In this case, the PCC voltage is out-of-control and its amplitude is determined by the RES output current, Irmax, and Rf.

If the amplitude of us remains higher than 0.9 p.u., the RES does not enter LVRT mode:

0.9UsnksjωsLacksisd+jisq+UrUr=Rfksisd+Irmax (3)

where Usn is the rated value of the amplitude of us; and isd is the d-axis current of is, which can be regarded as a rated value because the RES remains in normal operation mode.

Figure 4 shows the simulation results when Rf=8.5 Ω. The amplitude of usa drops from 1.04 to 0.96 p.u., and isq remains stable during the AC fault. Thus, RES does not enter LVRT mode. Because irdref=Irmax, the PI controller of the voltage loop is saturated, and the PCC voltage is out-of-control. Pr decreases to less than 0, and thus the power transmission is reversed in this case.

Fig. 4  Simulation results when Rf=8.5 Ω.

Because of the saturation of the voltage loop, the output value of the PI integrator keeps increasing and causes overvoltage of usa [

9]. Thus, isq is reversed and RES enters high-voltage ride-through (HVRT) after fault clearance, as shown in Fig. 4, where the anti-windup PI controller should be adopted to terminate the output value of the PI integrator and suppress PCC overvoltage.

D. RES LVRT Mode with RS Current Limitation Mode

When Rf is further decreased, the RES and RS output power can no longer maintain normal PCC voltage, and thus the RES begins to enter LVRT mode.

In this case, the fault point voltage is:

Ur=RfksIsn2-isq2+irmax (4)

where Isn is the rated value of is.

Based on the LVRT standard [

15], [16], the RES output reactive current is expressed as:

isq=-1.5×0.9-Us /UsnIsn (5)

According to (4) and (5), the relationship between Ur and Us must first be solved to obtain the critical value of Rf in this case.

The voltage equation given in Fig. 1 under the αβ frame is:

usαβks=Lacksdisαβdt+urαβ (6)

where subscript αβ denotes the αβ frame.

Equation (6) is transformed into the PLL frame:

usdq=ks2Lac(s+jωPLL)isdq+ksurdq (7)

where subscript dq denotes the PLL frame; and ωPLL is the angular frequency measured by the RES PLL.

Combined with the final value theorem, the fault steady- state equation of (7) becomes:

usdq=jωskbb2Lacisdq+ksurdq (8)

As us is oriented to the d-axis by the RES PLL, (8) can be simplified as:

Us=-ωsks2Lacisq+ksurd0=ωsks2Lacisd+ksurq (9)

where urd and urq are d- and q-axis components of ur, respectively.

Equation (9) can be revised as:

Us=ksUr2-ωsksLacisd2-ωsks2Lacisq (10)

As isd is calculated by Isn in the AC fault, we can obtain:

Us=ksUr2-ωs2ks2Lac2Isn2-isq2-ωsks2Lacisq (11)

Equation (11) can then be solved as:

isq=ks2Ur2-Isn2ωs2Lac2ks4-Us22ωsLacUsks2 (12)

When (5) is combined with (12), PCC voltage amplitude Ur under the fault steady state can be expressed as:

Ur=1ksUs2+Isn2ωs2Lac2ks4-3ωsLacUsIsnks20.9-UsUsn (13)

When (5) and (13) are substituted into (4), the relationship between Rf and Us can be solved.

Figure 5 shows the simulation results when Rf=5 Ω. The amplitude of usa is close to 0.59 p.u., which can be obtained by (13). Irdref reaches Irmax, and thus the RS AC voltage loop is saturated, and the PCC voltage is out-of-control in this case. isq increases in the opposite direction, indicating that the RES enters LVRT mode. The RS becomes the constant current source with the value of Irmax, and transmission power is reversed. Overvoltage is suppressed (compare this with Fig. 4) because anti-windup is added to the PI controller of the RS AC voltage loop.

Fig. 5  Simulation results when Rf=5 Ω.

E. Current Limitation Mode of RES and RS

When Rf further decreases, the RES output current reaches its limitation, and therefore, both the RES and RS become constant current sources. Substituting these constraints into (10) and (4) yields:

Us=ksUr-ωsLacisqks2Ur=RfIrmax (14)

Simulation results when Rf=0.01 Ω are shown in Fig. 6. Both isq and ird reach their limitations, and thus RES and RS can be treated as constant current sources with limitation values. Substituting these limitations into (14) obtains Us=0.184 p.u., which is consistent with the amplitude of usa shown in Fig. 6. As Rf0, PCC voltage is very close to 0 and thus Pr=0, as shown in Fig. 6.

Fig. 6  Simulation results when Rf=0.01 Ω.

The system operation mode under Rf variation is shown in Fig. 7. Based on (1)-(4), (13), (14), and Table I, the RS equilibrium curve under a symmetric fault is shown in Fig. 8. In Figs. 7 and 8, the green line indicates the area of controllable PCC voltage and continuous power transmission; the red line denotes the area of the controllable PCC voltage but interrupted power transmission; and the purple line indicates the area of the out-of-control PCC voltage and interrupted power transmission.

Fig. 7  System operation mode under Rf variation.

Fig. 8  RS equilibrium curve under a symmetric fault.

Figures 7 and 8 show that the system has five operation modes and three fault characteristics. Thus, the state space equations are discontinuous, making fault calculation and optimized design rather challenging. In addition, out-of-control PCC voltage and RS power transmission interruption occur under multiple modes.

IV. Asymmetric AC Fault Analysis

Asymmetric fault characteristics are studied based on the results presented in Fig. 8, and the P-G fault is adopted as an example for analysis.

Three-phase voltages, ua, ub, and uc are expressed as:

uaubuc=U+cosωst+θ+U+cosωst+θ+-2π/3U+cosωst+θ++2π/3+U-cosωst+θ-U-cosωst+θ-+2π/3U-cosωst+θ--2π/3 (15)

where the superscripts + and - denote the positive sequence (PS) and NS components, respectively; and U and θ are the magnitude and phase of the voltage, respectively.

The Clarke transformation is applied to (15), and the result can be rewritten in an exponential form:

uαβ=U+ejωst+θ++U-e-jωst+θ- (16)

Transforming (6) into the global PS and NS rotating frame yields:

udq=U+ejθ++U-e-jθ-=udq++u¯dq- (17)

where ()¯ means the conjugate. Equations (15)-(17) show that a conjugate operation is required to transform the NS component from the three-phase frame to the rotating frame. Based on the aforementioned transformations, the sequential network is built in the rotating frame to connect the control method of the RS and RES and the response of the fault point.

Reference [

11] recommends NS voltage suppression, where the sequence network under the P-G fault is shown in Fig. 9, in which the superscript z denotes the zero sequence (ZS) component; and Rsg and Rrg are the transformer grounding resistance at the RES and RS sides, respectively.

Fig. 9  Sequence network under a P-G fault when RS adopts NS voltage suppression.

Figure 9 shows that the current operation mode of the RES and RS must be judged in real time, and thus an iterative calculation is needed for fault analysis.

According to the superposition theorem, we can obtain:

us+=ks2is+Zac+ksur+u¯s-=ks2i¯dq-Zacir+=ksis+-ur+/Zdemi¯r-=ksi¯s--ur+/ZdemZdem=3Rf+3RrgZac+3Rsg3Rsg+Zac+3Rrg (18)

where Zac denotes the line impedance.

The absolute value of the NS voltage equation in (18) is:

Us-=ks2Rac2+ωsLac2i¯s- (19)

where Us- is the amplitude of u¯s-; and Rac is the resistance component of Zac.

The NS current specified by the grid code is:

isd-=0isq-=kiq-IsnUs-/Usn (20)

where kiq- is the NS reactive current coefficient.

Combining (19) and (20) yields:

Us-=ks2Rac2+ωsLac2kiq-IsnUs-/Usn (21)

Obviously, the solution to (21) is Us-=0. From (20), isq- equals 0. Therefore, no NS component exists in the RES during a P-G fault. However, an NS current is generated, which may cause RS overcurrent and saturated modulation during AC FRT.

The MMC valve-side voltage equation and modulation ratio Mr in Fig. 1 are expressed as:

ur/krejπ6=jωsLreqkrirp+umMr=2Um/(1.15udcn) (22)

where Lreq is the equivalent MMC inductor; um is the MMC valve-side voltage; Um is the amplitude of um; and irp is the current of the RS positive pole, which equals (ir++ir-)ejπ/6/2.

Solving (18) and (22) yields Mr and the amplitude of irp (Irp), as shown in Fig. 10 and verified by Fig. 11. In Fig. 11, is exhibits a pure sinusoidal shape, and thus the NS component is not contained. The amplitude of ir is 1.56 p.u., and RS modulation mr is saturated, which must be considered under FRT design.

Fig. 10  Irp and Mr under a P-G fault.

Fig. 11  Simulation results under a P-G fault (Rf=1 Ω).

V. Active PCC Voltage Drop Control

A. Control Design

Figure 7 shows that although RES has the modes of rated operation, LVRT, and current limitation under an AC fault, the response of the LVRT mode should strictly meet the FRT grid code, and the rated operation and current limitation modes cannot be easily optimized. Thus, the RES leaves little degree of freedom for improvement. For this reason, this study optimizes the FRT curve shown in Fig. 8 by changing the RS fault response.

From the analysis in Section III, if Urn is adopted as the reference of RS AC voltage control, the fault response is uniquely determined, as shown in Figs. 7 and 8, and thus the freedom of RS optimization remains lost. Therefore, the reference of RS AC voltage control is changed as an independent variable to provide the freedom for FRT optimization, after which an optimal voltage reference is determined by analyzing the influence of the RES on the RS during the AC fault. This is necessary when designing an enhanced FRT control and optimizing the equilibrium curve.

If the RS AC voltage loop is unsaturated, the general expression of transmitted power is:

Pr=1.5UsIsn2-isq2-1.5Ur2/Rf (23)

Combined with RES LVRT mode analysis, (5) and (13) are substituted into (23), and the result of Pr0 is plotted in Fig. 12.

Fig. 12  RS active power under LVRT mode.

In Fig. 12, the three-dimensional (3D) surface is the solution to (23) larger than 0, the solid grid plane is Pr=0, and the 3D curve equals the equilibrium curve shown in Fig. 8, where the pink, red, and green segments in the 3D curve correspond to the pink, red, and green modes in the left of Fig. 7, respectively. The solid and dashed segments in the 3D curve are below the solid grid plane, which means the RS input power is reversed under this part. However, the corresponding 3D surface remains, and thus continuous power transmission during the fault period is theoretically feasible. Therefore, continuous power transmission during an AC fault can be realized by moving the RS equilibrium point from the 3D curve to the 3D surface, where the reversed power shown in Figs. 3-5 and the red and pink modes presented in Figs. 7 and 8 can be effectively avoided. In addition, the optimal selection of the equilibrium point on the 3D surface can realize the maximum power transmission under a fault steady state.

Both the colored surface and green areas of the 3D curve in Fig. 12 satisfy Pr0, where the green areas of the 3D curve indicate that the RES and RS are in rated operation modes. Solving (23), i.e., PrPsn-1.5Urn 2/Rf to obtain optimal transmitted power yields the result shown in Fig. 13. The 3D curve is the maximum power transmission curve of the 3D surface and yields the desired equilibrium points for obtaining the maximum transmission under different Rf scenarios. The specific derivation of this curve can be found in Appendix A. To reduce the requirement for fault detection, the relationship between Pr and Ur on the solid line in Fig. 13 is fitted to Fig. 14.

Fig. 13  RS active power optimization.

Fig. 14  Proposed fault ride-through function.

In addition to the maximum power transmission, active PCC voltage drop can not only help reduce voltage error, desaturate the PI controller, and suppress the PCC overvoltage after fault clearance, but also decrease the current loop feedforward to avoid over-modulation, as shown in Fig. 11. When combined with (18), the NS current can also be suppressed after PCC voltage is reduced. Thus, the RS overcurrent shown in Fig. 10 can be significantly improved. Therefore, the aforementioned method is more advanced than the widely used current limitation mode of all sources under an AC fault.

To avoid algorithm switching during a fault, the fitting result shown in Fig. 14 is extended to Ur=1:

Ur=2.2Pr+0.11Urn    Pr[0.01,0.4]p.u. (24)

According to (24), when an AC fault occurs, Pr decreases immediately, and thus the proposed method actively reduces Ur to move the RS equilibrium point to the maximum power transmission curve shown in Fig. 13. When the fault is cleared, Pr increases immediately, and the proposed method actively increases Ur. From Fig. 14, Pr is further increased, and the positive feedback is constituted until Ur reaches Urn. Therefore, fault detection is not required. To avoid steady-state undervoltage under low wind speed, a trigger signal kt is added:

kt=1    Ur<0.9Urn  and  Pr<0.4Prn0    otherwise (25)

where 0.9 p.u. is the lower limitation of normal PCC voltage; and 0.4 p.u. is the maximum valne of Pr in (24).

Combining (24) and (25) enables us to generate the proposed method, given as Fig. 15, where kp is the coefficient of the minimum power constraint; and 0.13Urn is derived from the lower limitation of (24).

Fig. 15  Block of proposed method.

When Pr<0.01 p.u., the RS equilibrium point violates the constraint in (24), kp=0, and the proposed method is equal to the rated AC voltage control. Similarly, when Ur and Pr violate (25), kt=0 and the proposed method also becomes the rated AC voltage control to avoid steady-state PCC undervoltage. When kp=kt=1, the output of (24) is activated to output the reference of the d-axis RS voltage loop.

The design guidelines of the proposed method are presented in Fig. 16. It can be observed that the proposed method does not require either fault detection, Rf measurement, operation mode identification, or parameter tuning, and the method is very easy to implement.

Fig. 16  Design guidelines of proposed method.

B. Effective Range Calculation and Contribution Analysis

The effective range of the proposed method is shown in the green area of Fig. 7, and its effect on the continuous power supply is presented in Fig. 17.

Fig. 17  Effect of proposed method on continuous power supply.

The dotted-dashed and dotted lines indicate the RS equilibrium points under rated AC voltage control, where the dotted-dashed line indicates the continuous power transmission area. The green area represents the continuous power transmission area of the proposed method, which is a projection of the surface shown in Fig. 12. The solid line in this green area is the optimized RS equilibrium curve, which is a projection of the maximum power transmission curve shown in Fig. 13.

As Fig. 17 shows, the RS continuous power transmission area can be significantly expanded by the proposed method. Combined with Figs. 3-6, the RS equilibrium points in the aforementioned four cases are all located in the power transmission interruption area of Fig. 17, which can be proven based on the reversed Pr from the respective simulation under fault. To achieve continuous power transmission during the fault, the aforementioned equilibrium points must be dropped to the green area, or most specifically to the solid line in this green area, to maintain the maximum power transmission under fault. It can be concluded that excessive PCC voltage amplitude is the main cause of RS voltage loop saturation and power transmission interruption and can be optimized by the proposed method.

C. Robustness Analysis

With the results in Figs. 9-11 calculated by (5), (13), and (23), Lac is the most likely parameter to introduce a calculation error. Therefore, the robustness of the proposed method under Lac variation is verified in Fig. 18, where the pink 3D curve is equal to that shown in Fig. 13. Combined with the green area in Fig. 7, the effective area decreases from Rf3.83 Ω to only Rf7.02 Ω under 20% Lac variation. Therefore, the proposed method is effective under severe measurement distortion.

Fig. 18  Proof of robustness under Lac variation.

In addition, if the topology in Fig. 1 changes and an additional 0.25 p.u. RES is connected to the PCC, the original maximum power transmission curve as shown in Fig. 19 exhibits a small failure interval of 1.66 ΩRf2.83 Ω, which can be easily compensated by updating the calculation data.

Fig. 19  Proof of robustness under a newly connected RES.

As the wind speed usually changes in practice, the proposed method should be tested under a steady-state low-power operation. The results are shown in Fig. 20, where Psmax is the RES maximum output active power. When the maximum output active power decreases from 0.75 p.u. to 0.3 p.u. due to the decrease in wind speed, the effective range of the proposed method remains unchanged, but the transmitted power reduces after 15.2 Ω.

Fig. 20  Proof of robustness under a steady-state low-power operation.

Therefore, the proposed method is robust under inaccurate parameter measurements, unplanned RES expansion, and wide variations in wind speed.

VI. Simulation verification

To prove the veracity of the analysis presented in Figs. 7, 8, 15, and 16, the proposed method shown in Fig. 14 is applied to all cases shown in Figs. 2-5. In addition, the performance of the proposed method under an asymmetric fault is also tested and analyzed. In each simulation case, the dotted and solid lines in the figures denote the results of rated AC voltage control (original method) and of the proposed method, respectively. The system as shown in Fig. 1 enters into steady state at 0.4 s, and an AC fault occurs at 0.6 s and is cleared at 1 s. In Figs. 21-27, idiff is the MCC circulating current, usm is the voltage of submodule capacitors, and udc is the DC bus voltage.

Fig. 21  Simulation comparisons when Rf=30 Ω. (a) Fault response. (b) MMC internal states under proposed method.

Fig. 22  Simulation comparisons when Rf=15 Ω. (a) Fault response. (b) MMC internal states under proposed method.

Fig. 23  Simulation comparisons when Rf=8.5 Ω. (a) Fault response. (b) MMC internal states under proposed method.

Fig. 24  Simulation comparisons when Rf=5 Ω. (a) Fault response. (b) MMC internal states under proposed method.

Fig. 25  Verification under a P-G fault (Rf=1 Ω). (a) Fault response. (b) MMC internal states under proposed method.

Fig. 26  Verification under a 2P-G fault (Rf=5 Ω). (a) Fault response. (b) MMC internal states under proposed method.

Fig. 27  Verification under a newly connected RES (Rf=5 Ω).

A. Comparisons with Case Shown in Fig. 2

The dotted lines in Fig. 21(a) indicate the results that are equivalent to those shown in Fig. 2. Figure 21(a) shows that usa is reduced by the proposed method during an AC fault. According to (13), Ur equals 0.63 p.u. in this case, which corresponds to the black point on the solid curve shown in Fig. 17. Here, isq increases in the reverse direction, which means that the proposed method drives the RES into LVRT mode. Under the proposed method, the transmitted power during the fault increases from 0.15 to 0.24 p.u., which is consistent with the solid line shown in Fig. 20. From Fig. 21(b), MMC internal states are stable under the duration of the fault, and the DC bus voltage can be maintained in the normal range. Thus, the DC grid has little effect on the AC bus.

B. Comparisons with Case Shown in Fig. 3

The dotted lines in Fig. 22(a) indicate the results that are equivalent to those shown in Fig. 3. According to Fig. 18(a) and (b), the amplitude of usa is lower than 0.9 p.u. and isq increases in reverse direction during the fault. Thus, the RES operates in LVRT mode. According to (13), Ur is 0.33 p.u. in this case, which corresponds to the black point on the solid curve shown in Fig. 17. Under the proposed method, the obtained active power of the RS during the AC fault increases from -0.43 to 0.1 p.u., which is consistent with the curve shown in Fig. 14. The reversed irdref is suppressed by the proposed method to prevent the power from flowing into the fault point.

C. Comparisons with Case Shown in Fig. 4

The dotted lines in Fig. 23(a) indicate the results that are equivalent to those shown in Fig. 4. As Fig. 23(a) shows, under the proposed method, the amplitude of usa is reduced during the fault and restored without overvoltage at fault clearance. According to (13), Ur equals 0.21 p.u. in this case, which corresponds to the black point on the solid curve shown in Fig. 17. Here, isq increases in a reverse direction during the fault, and thus the RES operates in LVRT mode. From the power waveforms, the proposed method can not only avoid RS power reversal but can also receive 0.042 p.u. of active power during the fault, which is consistent with the solid line in Fig. 20. In addition, under the proposed method, the RS AC voltage loop is no longer saturated, and thus PCC voltage is controllable during the fault period. Combined with the analysis under (3), irdref is no longer greater than 0 at fault clearance after voltage loop desaturation. Thus, the steep rise in modulation amplitude stops, and PCC overvoltage is clearly suppressed by the proposed method.

D. Comparisons with Case Shown in Fig. 5

The dotted lines in Fig. 24(a) indicate results that are equivalent to those shown in Fig. 5. It can be observed that the amplitude of usa is further reduced and isq is further increased in the reverse direction under the proposed method, which can optimize the RES LVRT mode. According to (13), Ur equals 0.14 p.u. in this case, which is consistent with the data point on the solid curve shown in Fig. 17. The power reversal is optimized by the proposed method, and the RS receives 0.013 p.u. of active power during the AC fault, which is consistent with the curve shown in Fig. 14. In addition, under the proposed method, the RS voltage loop is no longer saturated, and thus PCC voltage is controllable during the fault, and PCC overvoltage is obviously suppressed.

E. Verification Under a P-G Fault

The proposed method was verified under a P-G fault, and the results are shown in Fig. 25. Since the RS adopts NS voltage suppression as shown in Fig. 1, NS voltage is suppressed but ZS voltage remains, and thus, PCC voltage is not sinusoidal. Under the proposed method, AC bus overcurrent is effectively suppressed and the transmitted power increases under the fault, thus showing that the proposed method is effective under a P-G fault. RES terminal voltage is sinusoidal because ZS voltage is blocked by the YD11 transformer. The terminal current is also sinusoidal due to NS negative current suppression. Figure 25(b) shows that circulating current waveforms are stable but do not overlap, which is consistent with the results presented in [

27].

F. Verification Under a Phase-to-phase-to-ground (2P-G) Fault

The proposed method was verified under a 2P-G fault, and the results are shown in Fig. 26. In this case, PCC PS voltage is equal to NS voltage [

11], but PS voltage is controlled to be the rated value, whereas NS voltage is controlled to be 0 by the RS. Thus, a control conflict is initiated and RS modulation is saturated. When the proposed method is used, overcurrent is suppressed and modulation is desaturated, and thus RS control capability is restored. Reverse power flow is suppressed and continuous power transmission is maintained. The proposed method is effective under a 2P-G fault.

Although the ZS voltage and current are 0, the responses of a phase-to-phase (P-P) fault are similar to those of a 2P-G fault [

11]. Accordingly, the same analysis and results are not repeated here.

G. Verification Under a Newly Connected RES

Under a newly connected RES, the topology in Fig. 1 changes and an additional 0.25 p.u. RES is connected to the PCC. Results are presented in Fig. 27, where subscript i indicates the ith RES. Figure 27(a) and (b) shows that both RESs exhibit similar responses under the fault because they must fulfill the FRT grid code. Here, Ps is greater than 0 under the duration of the fault, which is consistent with Fig. 19. DC bus voltage can be controlled under a normal range. Thus, the AC fault is decoupled with a DC-side response, which can be simplified by an ideal voltage source under fault.

H. Comparisons with Existing Method

In this section, the latest version of the voltage-optimized FRT method [

11] is introduced and compared with the proposed method. A control block is presented in Appendix A, where Ur- and Urz are the NS and ZS amplitudes of ur, respectively.

As the NS and ZS values are 0 under a symmetric fault, the method presented in Appendix A shows the similar response as the rated voltage control and its current loop. Thus, the response of this method is equal to that shown in Figs. 2-6 under fault resistance variation, and the superiority of the proposed method is well demonstrated in Figs. 21-24.

Figure 28 presents the results of the method in [

11] under a P-G fault when the RES adopts NS current suppression. Figure 28(a) shows that RES PLL output frequency fpll is stable and RES terminal waveforms are sinusoidal under the fault. The d-axis NS RES voltage usd- is not zero, but is- maintains zero output because the RES adopts NS current suppression. As Fig. 28(b) shows, Pr+=0.25 p.u., which is 5% greater than the result presented in Fig. 25(a). The d-axis NS current of RS ird- reaches its limited value (0.25 p.u.), and thus the negative voltage loop presented in Fig. A6 is saturated and the d-axis NS voltage urd- is not zero.

Fig. 28  Results of method in [

11] under a P-G fault (Rf=1 Ω) when RES adopts NS current suppression. (a) RES response. (b) RS response.

When FRT grid code is used in the RES, (20) requires that the RES output an NS current according to its NS terminal voltage, which changes the fault point voltage and in turn generates RES voltage regulation and current redistribution, resulting in a fault response interaction. From (19)-(21), the NS RES current is maintained at 0 in all cases when the RS adopts NS voltage suppression. Thus, the interaction is blocked. However, both the RES and RS become NS current sources when the method in [

11] is used, which exacerbates the fault interactions and results in divergence. This can be verified by Fig. 29, where the RES and RS adopt FRT code and the method in [11], respectively, under a P-G fault.

Fig. 29  Results of method in [

11] under a P-G fault (Rf=1 Ω) when RES adopts FRT code.

Thus, although the method in [

11] can improve the power transmission by 5% when the RES adopts NS current suppression, it cannot achieve performance optimization under a symmetric fault; and in fact, it leads to divergence when the RES adopts FRT code under an asymmetric fault.

VII. Conclusion

This study analyzed the AC fault responses of MMC-integrated wind farms. It showed that the integrated system has five operation modes and three common fault characteristics under various fault resistances at the fault point. Wherein out-of-control PCC voltage and RS power reversal occur under multiple operation modes. This study also determined the RS continuous power transmission area during a fault and revealed that excessive PCC voltage amplitude is the primary cause of these instances of performance degradation. In addition, the active PCC voltage drop can greatly expand the continuous power transmission area under an AC fault. A novel FRT method was then proposed to realize the maximum power transmission and controllable PCC voltage under wide fault resistance variation. The method was shown to be robust under inaccurate line parameter measurements, unplanned RES expansion, and wide variations in wind speed.

Appendix

APPENDIX A

A. RES Control Method

The blocks of RES control methods in PS and NS are shown in Figs. A1 and A2, respectively.

Fig. A1  Block of RES control method in PS.

Fig. A2  Block of RES control method in NS.

B. MMC Control Methods

Fig. A3  Block of PS voltage control and its current loop.

Fig. A4  Block of NS voltage suppression control and its current loop.

Fig. A5  Block of circulating current suppression control.

C. Derivation of Maximum Power Transmission Curve

The mathematical expression of Pr surface is:

Pr=-1.5Ur2Rf+0.752.7a+bUsn+3a/UsnIsn21-2.250.9-1.35a+0.5b/Usn2+3a2a=Isnks2LacUsnωsb=7.29a2+4Usn+3a/UsnUsn3Ur2-a2/Usn (A1)

Based on Pr=0, the continuous power transmission area in the (Rf, Ur) plane can be obtained.

Solving Pr/Ur=0 and substituting the result into Pr, we can obtain the maximum curve in the Rf direction. A detailed expression is given by (A2).

Pr=AB+0.752.7a+bUsn+3a/UsnIsn21-2.25×0.9-(1.35a+0.5b)/(Usn2+3a)]2A=1.46×1015Isn2Usn2Ur2Usn30.18+Ur2+aUsn-0.73+3Ur2-3a3/Usn3-4.1×10-7a/Usnc+Usn-1.18a2/Usn2-4.1×10-7cB=Usn2+6a+9a2/Usn2cIsn2Usn3-1.43×1016-3.9×1016Ur2+aUsn9.22×1015-1.17×1017Ur2+d/UsnUsn+3a/Usn2c=Usn1.08×1013Usn3Ur2+3.25×1013aUsnUr2+8.9×1012a2/Usn-3.25×1013a3/Usn3d=1.17×1017a3/Usn3+3.2×1010a/Usnc+Usn5.28×1016a2/Usn2+2.13×1010c (A2)

D. Control Block of Method in [

11]

Fig. A6  Control block of method in [

11].

References

1

H. Yuan, H. Xin, L. Huang et al., “Stability analysis and enhancement of type-4 wind turbines connected to very weak grids under severe voltage sags,” IEEE Transactions on Energy Conversion, vol. 34, no. 2, pp. 838-848, Jun. 2019. [Baidu Scholar] 

2

M. Talha, S. R. S. Raihan, and N. A. Rahim, “A grid-tied PV inverter with sag-severity-independent low-voltage ride through, reactive power support, and islanding protection,” Journal of Modern Power Systems and Clean Energy, vol. 9, no. 6, pp. 1300-1311, Nov. 2021. [Baidu Scholar] 

3

Y. Qu, L. Gao, G. Ma et al., “Conjoint analysis of crowbar switching scheme and its resistance based on statistical sampling for LVRT of DFIG,” Journal of Modern Power Systems and Clean Energy, vol. 7, no. 3, pp. 558-567, May 2019. [Baidu Scholar] 

4

J. Yang, S. Mao, X. Li et al., “Influence mechanism and suppression control of the MMC on short-circuit current under AC faults,” IEEE Access, vol. 8, pp. 138307-138317, Jul. 2020. [Baidu Scholar] 

5

J. Moon, J. Park, D. Kang et al., “A control method of HVDC-modular multilevel converter based on arm current under the unbalanced voltage condition,” IEEE Transactions on Power Delivery, vol. 30, no. 2, pp. 529-536, Apr. 2015. [Baidu Scholar] 

6

S. Wang, G. P. Adam, A. M. Massoud et al., “Analysis and assessment of modular multilevel converter internal control schemes,” IEEE Journal of Emerging and Selected Topics in Power Electronics, vol. 8, no. 1, pp. 697-719, Mar. 2020. [Baidu Scholar] 

7

Y. Li, X. Wang, J. Guo et al., “PLL synchronization stability analysis of mmc-connected wind farms under high-impedance AC faults,” IEEE Transactions on Power Systems, vol. 36, no. 3, pp. 2251-2261, May 2021. [Baidu Scholar] 

8

H. Ye, W. Cao, W. Chen et al., “An AC fault ride through method for MMC-HVDC system in offshore applications including DC current-limiting inductors,” IEEE Transactions on Power Delivery, vol. 37, no. 4, pp. 2818-2830, Aug. 2022. [Baidu Scholar] 

9

Y. Li, J. Guo, X. Zhang et al., “Over-voltage suppression methods for the MMC-VSC-HVDC wind farm integration system,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 67, no. 2, pp. 355-359, Feb. 2020. [Baidu Scholar] 

10

X. Wang, L. Yang, B. Lin et al., “Mechanism and restraining strategy of the sending-end grid overvoltage in offshore wind farm-flexible HVDC transmission system under grid faults,” High Voltage Engineering, vol. 47, no. 8, pp. 2688-2699, Aug. 2021. [Baidu Scholar] 

11

L. Shi, G. P. Adam, R. Li et al., “Control of offshore MMC during asymmetric offshore AC faults for wind power transmission,” IEEE Journal of Emerging and Selected Topics in Power Electronics, vol. 8, no. 2, pp. 1074-1083, Jun. 2020. [Baidu Scholar] 

12

Y. Jing, R. Li, L. Xu et al., “Enhanced AC voltage and frequency control of offshore MMC station for wind farm connection,” IET Renewable Power Generation, vol. 12, no. 15, pp. 1771-1777, Sept. 2018. [Baidu Scholar] 

13

J. Freytes, S. Akkari, P. Rault et al., “Dynamic analysis of MMC-based MTDC grids: use of MMC energy to improve voltage behavior,” IEEE Transactions on Power Delivery, vol. 34, no. 1, pp. 137-148, Feb. 2019. [Baidu Scholar] 

14

L. Shi, G. P. Adam, R. Li et al., “Enhanced control of offshore wind farms connected to MTDC network using partially selective DC fault protection,” IEEE Journal of Emerging and Selected Topics in Power Electronics, vol. 9, no. 3, pp. 2926-2935, Jun. 2021. [Baidu Scholar] 

15

Technical Rule for Connecting Wind Farm to Power System, GB/T 19963-2011, China, 2011. [Baidu Scholar] 

16

Technical Requirements for Connecting Photovoltaic Power System to Power System, GB/T 19964-2012, 2012. [Baidu Scholar] 

17

X. Sun, Q. Zhang, B. Ren et al., “A novel digital phase-locked loop for distributed grid-connected power generation systems under unbalanced and distorted utility conditions,” in Proceedings of the 2nd International Symposium on Power Electronics for Distributed Generation Systems, Hefei, China, Aug. 2010, pp. 813-817. [Baidu Scholar] 

18

W. Zhou, “The algorithm research about three-level APF virtual flux models predict SVPWM and FPGA implementation,” M.S. thesis, Department of Electronic Engineering, Lanzhou University of Technology, Lanzhou, China, 2018. [Baidu Scholar] 

19

J. Freytes, J. Li, G. de Préville et al., “Grid-forming control with current limitation for MMC under unbalanced fault ride-through,” IEEE Transactions on Power Delivery, vol. 36, no. 3, pp. 1914-1916, Jun. 2021. [Baidu Scholar] 

20

H. Zeng, F. Lin, and X. Kang. “Distributed parameter based distance protection for transmission lines and identification of short-circuit faults under low-, medium- and high-resistance at fault point,” Power System Technology, vol. 38, no. 1, pp. 262-268, Jan. 2014. [Baidu Scholar] 

21

M. Langwasser, G. de Carne, M. Liserre et al., “Fault current estimation in multi-terminal HVDC grids considering MMC control,” IEEE Transactions on Power Systems, vol. 34, no. 3, pp. 2179-2189, May 2019. [Baidu Scholar] 

22

X. Long, Y. Li, W. Xu et al., “A new technique to detect faults in de-energized distribution feeders ‒ Part II: symmetrical fault detection,” IEEE Transactions on Power Delivery, vol. 26, no. 3, pp. 1902-1910, Jul. 2011. [Baidu Scholar] 

23

S. Chatterjee and B. K. S. Roy. “Fast identification of symmetrical or asymmetrical faults during power swings with dual use line relays,” CSEE Journal of Power and Energy Systems, vol. 6, no. 1, pp. 184-192, Mar. 2020. [Baidu Scholar] 

24

M. G. Taul, X. Wang, P. Davari et al., “An overview of assessment methods for synchronization stability of grid-connected converters under severe symmetrical grid faults,” IEEE Transactions on Power Electronics, vol. 34, no. 10, pp. 9655-9670, Oct. 2019. [Baidu Scholar] 

25

G. Bergna-Diaz, J. A. Suul, E. Berne et al., “Optimal shaping of the mmc circulating currents for preventing AC-side power oscillations from propagating into HVDC grids,” IEEE Journal of Emerging and Selected Topics in Power Electronics, vol. 7, no. 2, pp. 1015-1030, Jun. 2019. [Baidu Scholar] 

26

W. Li and J. Bélanger, “An equivalent circuit method for modelling and simulation of modular multilevel converters in real-time HIL test bench,” IEEE Transactions on Power Delivery, vol. 31, no. 5, pp. 2401-2409, Oct. 2016. [Baidu Scholar] 

27

J. Wang, J. Liang, C. Wang et al., “Circulating current suppression for MMC-HVDC under unbalanced grid conditions,” IEEE Transactions on Industry Applications, vol. 53, no. 4, pp. 3250-3259, Mar. 2017. [Baidu Scholar] 

28

R. Li, J. E. Fletcher, and B. W. Williams. “Influence of third harmonic injection on modular multilevel converter-based high-voltage direct current transmission systems,” IET Generation, Transmission & Distribution, vol. 10, no. 11, pp. 2764-2770, Jun. 2016. [Baidu Scholar]