Microsoft Word - 10-3343_s_ETASR_V10_N2_pp5396-5401 Engineering, Technology & Applied Science Research Vol. 10, No. 2, 2020, 5396-5401 5396 www.etasr.com Lam et al.: Simulation Models for Three-phase Grid-connected PV Inverters Enabling Current … Simulation Models for Three-phase Grid-connected PV Inverters Enabling Current Limitation under Unbalanced Faults Le Hong Lam Faculty of Electrical Engineering The University of Danang—University of Science and Technology Da Nang, Vietnam lhlam@dut.udn.vn Tran Dai Hoang Phuc Faculty of Electrical Engineering The University of Danang—University of Science and Technology Da Nang, Vietnam 105150107@sv.dut.edu.vn Nguyen Huu Hieu Faculty of Electrical Engineering The University of Danang—University of Science and Technology Da Nang, Vietnam nhhieu@dut.udn.vn Abstract—Normally unbalanced grid voltage dips may lead to unbalanced non-sinusoidal current injections, DC-link voltage oscillations, and active and/or reactive power oscillations with twice the grid fundamental frequency in three-phase grid- connected Photovoltaic (PV) systems. Double grid frequency oscillations at the DC-link of the conventional two-stage PV inverters can further deteriorate the DC-link capacitor, which is one of the most important limiting components in the system. Proper control of these converters may efficiently address this problem. In such solutions, Current Reference Calculation (CRC) is one of the most important issues that should be coped with the reliable operation of grid-connected converters under unbalanced grid faults. Therefore, this paper proposes and simulates CRC methods and presents the results in order to improve the quality of the grid-connected PV system under unbalanced grid voltage fault. Keywords–Photovoltaic (PV) systems; unbalance voltage; two– stage converters; power oscillation; DC–link voltage oscillation I. INTRODUCTION With the fast increase of grid-connected PV generation [1, 2], PV systems should contribute to the grid stability by providing ancillary services, beyond the basic power delivery. The new grid requirements demand the grid-connected PV systems, single- or three-phase, to have the capability to operate in power factors other than unity. Also, based on the recently revised grid codes, PV inverters are preferred to stay connected during grid voltage faults [3-6]. When a fault happens, the converter has to detect the incident and respond quickly to the disturbance to mitigate the adverse effects on the inverter, the equipment connected to the grid, and the upstream system. Indeed, the revised grid codes require the PV systems to inject a certain amount of reactive power in case of low voltage fault [7, 8]. These issues are now gaining more consideration in PV systems, as the power capacity of an individual PV system is also increasing. Detection of voltage sags, current limitation, current reference calculation [9], active and reactive power oscillation, and DC-link voltage oscillation are such important issues. Besides, they are the key issues to the proper operation of grid-connected PV converters under faults. Among them, CRC plays the most important role to satisfy the grid requirements, especially under unbalanced grid faults. In current researched methods, d-q methodology is commonly implemented under grid voltage faults [3-6, 10] although it is quite complex since it requires building blocks to convert the signals and calculate the frequency as the Phase Lock Loop (PLL), Dual Second Order Generalized Integrator Based Frequency-Locked Loop (DSOGI–FLL). The aforementioned active power oscillation can have a negative impact on the reliable operation of the grid-connected PV converters. In two-stage PV converters, where a DC-DC converter operates as Maximum Power Point Tracking (MPPT) [11], it is common that a PI controller determines the active power reference. Thus, in case that the injected active power starts fluctuating, the PI controller cannot follow its sinusoidal variations because the PV power injected to the DC-link is constant. As a result, the DC-link voltage will fluctuate with the same frequency of the injected active power [12]. Notably, due to the high failure rates of the electrolytic capacitors of the two-stage PV converters, the system reliability is challenged. This is worsened by DC-link voltage ripples. In this paper there are two main contributions: (i) DC link ripples during unbalanced faults are reduced with proper control of DC-DC converter and (ii) the impact of the PV systems on the distribution grid is analyzed and a control strategy is proposed for the PV using calculation methods referenced in the static α- β reference frame to reduce the complexity of the control structure and more efficient instead of using moving d-q reference frames. Moreover, the deprecating power of the main grid can be made by changing the power factor of wind generators [13-15]. This control strategy may offer a new solution to solve the voltage sag issue on a power system due to the PV plants connected to the grid by streaming their reactive power to the grid to enhance its power quality. II. SYSTEM OPERATION This section analyzes the inverter operation under normal and abnormal conditions for a three-wire three-phase PV system. The two-stage three-phase system shown in Figure 1 Corresponding author: Le Hong Lam Engineering, Technology & Applied Science Research Vol. 10, No. 2, 2020, 5396-5401 5397 www.etasr.com Lam et al.: Simulation Models for Three-phase Grid-connected PV Inverters Enabling Current … includes a boost converter and a full-bridge inverter interconnected through the DC-link capacitor. PV ARRAY BOOST CONVERTER V dc 2 V dc 2 C C dc dc S4 S3 S6 S2 S5S1 a c b LCL FILTER Va V b Vc Fig. 1. Two-stage three-phase grid-connected PV system. The formulation is performed in the Stationary Reference Frame (SRF). The conversion from the three-phase system into the SRF is as: ��� = ������ = �� 1 − �� − ��0 √ � − √ � � � ������� (1) where ��, �� are the voltages in the SRF from ��, ��, ��. The apparent power S is written as: � = �. �∗ = � + �� (2) Since under normal conditions the grid voltages and loads are balanced, there will not be any oscillatory components in the active and reactive components of the power and the injected current is completely sinusoidal. However, under unbalanced conditions, the NS components will appear in both current and voltage vectors. Thus, the apparent power is re- written as: � = ���. ���∗ = ���! + ���" #. ���! + ���" #∗ = ���! . ���! ∗ +���! . ���" ∗ + ���" .���! ∗ + ���" . ���" ∗ (3) in which ���! , ���" are derived from: ���! = �� $1 − %% 1& ��� (4) ���" = �� $1 %−% 1& ��� (5) where % = '"()/� is a 90°-lagging phase-shifting operator applied to the time domain [16]. Similarly, ���! , ���" are achieved following (4) and (5): ���! = �� $1 − %% 1& ��� (6) ���" = �� $1 %−% 1& ��� (7) In (3), there are four terms in the apparent power formulation. ���! . ���! ∗ = ��! + ���!#. ��! + ���!#∗ = ��!��! + ��!��! +� ��!��! − ��!��!# = �� + ��� (8) ���! . ���" ∗ = ��! + ���!#. ��" + ���"#∗ = ��!��" + ��!��" +� ��!��" − ��!��"# = �� + ��� (9) ���" . ���! ∗ = ��" + ���"#. ��! + ���!#∗ = ��"��! + ��"��! +� ��"��! − ��"��!# = � + �� (10) ���" . ���" ∗ = ��" + ���"#. ��" + ���"#∗ = ��"��" + ��"��" +� ��"��" − ��"��"# = �+ + ��+ (11) The constant and oscillating parts of the total active and reactive power are written as: � = �, + �- (12) �, = �� + �+ = ��!��! + ��!��! + ��"��" + ��"��" (13) �- = �� + � = ��!��" + ��!��" + ��"��! + ��"��! (14) � = �, + �- (15) �, = �� + �+ = ��!��! − ��!��! + ��"��" − ��"��" (16) �- = �� + � = ��!��" − ��!��" + ��"��! − ��"��! (17) where � and � are the total active and reactive power, and �,, �,, �-, �- are the constant and oscillating parts in the active and reactive power. Under balanced voltage sag faults, there is no NS in the voltages and currents, thus there are no oscillatory components in the active and reactive power. However, during unbalanced faults, the NS components appear in the voltages and currents. From (12) and (15), it is concluded that the active and reactive power have a constant part named �0 and �0. Also, there are two oscillating parts in the active and reactive power, denoted as �. and �.. Fundamentally, in the PV power systems, all the active power generated by the PV panels is delivered to the DC-link. This active power is continuously processed by the inverter and injected into the grid. If the active power generated by the inverter is less than the power injected to the DC link from the PV source, the DC link voltage will increase. Proper control is needed to synchronize the power flow from the PV source to the grid by regulating the DC link voltage. Accordingly, in the case that the injected active power has double grid frequency oscillations, the DC-link voltage will inevitably oscillate with the same frequency. Double grid frequency oscillations of the DC-link voltage have a negative impact on the life cycle of the capacitive DC-link. III. PROPOSED CONTROL STRATEGY In this section, a new control strategy to overcome the unbalance grid voltage based on the Current Reference Generation (CRG) is presented and a current limitation scheme is used to warranty the overcurrent issue. A. Current Reference Generation In order to eliminate active power oscillations, (14) has to be zero. In (14) and (17), �, and �, are equal to the average active power ( �/01 which is the output of the DC-link regulator) and reactive power (�/01 which is calculated during grid faults) references. These values are continuously calculated by the control algorithm. The control goal is to eliminate the oscillatory components from the active power, while allowing reactive power to oscillate with the double grid Engineering, Technology & Applied Science Research Vol. 10, No. 2, 2020, 5396-5401 5398 www.etasr.com Lam et al.: Simulation Models for Three-phase Grid-connected PV Inverters Enabling Current … frequency. Hence, �- is considered equal to 2 ��!��" − ��!��"#. Accordingly, (17) is rewritten as: �- = ��!��" − ��!��" + ��"��! − ��"��! = 2 ��!��" − ��!��"# (18) ��!��" − ��!��" − ��"��! + ��"��! = 0 (19) Furthermore, (13), (14), (16), (19) are written as: 34 44 5 ��! ��! ��" ��"��" ��" ��! ��!��! − ��! ��" − ��"−��" ��" ��! − ��! 67 77 8 344 45��!��!��"��"677 78 = 9�/010�/010 : (20) A formulation for generating sinusoidal currents to deliver a certain amount of active and reactive power [3, 17, 18] is obtained as: ��! = ;<=>?">@ �/01 − ;?!>@ �/01 (21) ��" = − ;?">@ �/01 − ;?!>@ �/01 (22) ��! = ;C=>?">@ �/01 − ;CA=>?!>@ �/01 (23) ��" = − ;CB>?">@ �/01 − ;CAB>?!>@ �/01 (24) ���D��D� = �0 − 11 0� ������ (25) EF = ��!� + ��!� (26) EG = ��"� + ��"� (27) where ��D, ��D are the orthogonal voltages of the SRF voltage vectors. The SRF currents are driven from the average value of the active and reactive power. These references determine the peak-peak value of the oscillations on the reactive power. Accordingly, a general formulation is obtained as: ��H = ;<=";?"N>@>OPQR � (32) where � is the apparent power or the nominal power of the power converter, E��S0 is the base voltage, which is equal to the Root-Mean-Square (RMS) value of the line-line grid voltage. On the other hand, according to the voltage sag depth, the reactive power can be calculated as: T�/01 = 0 EFU > 0.9�/01 = � × 1.5 × 0.9 − EFU# 0.2 < EFU < 0.9�/01 = 1.05 × � EFU < 0.2 (33) With EFU being calculated as: EFU = �;OPQR (34) Given the NNP and reactive power of Q, the maximum allowed active power �[�\ for the inverter to inject to the grid while avoiding overcurrent can be achieved as: �[�\ = NMM�� − �� (35) For operation of the converter under very deep voltage sags, NNP will have a low value, since NEF − NEG becomes small. Therefore, under a deep voltage sag, the condition is: If � > MM� → � = MM� and �[�\ = 0 (36) Vpu Pmax P* Vpu<0.9 No Yes Fault Signal =1 Pmax < P* Enable Signal = 0 Compare Signal = 1 No Yes Calculate Pmax, NNP, Qref Pmax Non - MPPT MPPT Enable Signal = 1 Fig. 2. Flowchart of the proposed control algorithm The DC-DC converter operates as the MPPT–P&O, the DC-DC converter should switch to the Non-MPPT mode in case a grid fault occurred and the inverter could not inject the maximum PV power. The flowchart in Figure 2 summarizes and clarifies the control system. If EFU < 0.9 the voltage sag detection block will generate a fault signal activating the MM�, �[�\, �/01 calculator block. Then, for a comparison between �[�\ and �∗, a comparator signal will be generated There is an AND block, in which, if the comparator signal and Engineering, Technology & Applied Science Research Vol. 10, No. 2, 2020, 5396-5401 5399 www.etasr.com Lam et al.: Simulation Models for Three-phase Grid-connected PV Inverters Enabling Current … fault signal are equal to 1, the DC-DC converter switches to the Non-MPPT. In this case, the fault signal is 1, while the comparator signal remains zero. MPPT may continue working under abnormal operation when the fault exists in the grid. IV. CONTROL STRUCTURE DESIGN A control structure was implemented in MATLAB/Simulink [19]. It should be noted that the proposed strategy control is mainly in Voltage Source Inverter (VSI) and current controller. Since, these modules play an important role in creating the reference current to balance voltage after faults. A. Voltage Source Inverter In connected grid mode, the VSI controller [20] power flow comes from the PV system to respond to the requirements of the grid, and control the voltage on the capacitor abc of DC– link. Active power is required from the grid reference given �/01 from Ebc voltage regulator. Meanwhile, reactive power �/01 is calculated according to the current limit to stabilize the voltage on the AC bus. Figure 3 shows the schematic diagram of the Ebc voltage control loop, the quantity to be controlled because it is proportional to the energy stored in the capacitor. DC voltage on the capacitor is stable thanks to the energy balance through (ignoring losses in VSI): �H> = �bc + �d (37) where �H> is the power from the PV, �bc is the power on the capacitor, and �d is the power inject to the grid. Vdc Vdc_ref (.) 2 (.) 2 + - Kv(s) + - PPV P ref Fig. 3. Diagram control Eef voltage When the voltage on the DC capacitor is stable, the value of the �bc is considered as zero: �H> = �d (38) Equation (38) indicates that the controlling �d will be controlled through the �d to achieve the desired value. Then, the referenced currents g�/01 and g�/01 are calculated through the reference power value from (39) and (40) and then they are taken to the current control. ��/01 = ;<=";