#### Abstract

This study presents a new active and reactive power control scheme for a single-stage three-phase grid-connected photovoltaic (PV) system during grid faults. The presented PV system utilizes a single-stage three-phase current-controlled voltage-source inverter to achieve the maximum power point tracking (MPPT) control of the PV panel with the function of low voltage ride through (LVRT). Moreover, a formula based on positive sequence voltage for evaluating the percentage of voltage sag is derived to determine the ratio of the injected reactive current to satisfy the LVRT regulations. To reduce the risk of overcurrent during LVRT operation, a current limit is predefined for the injection of reactive current. Furthermore, the control of active and reactive power is designed using a two-dimensional recurrent fuzzy cerebellar model articulation neural network (2D-RFCMANN). In addition, the online learning laws of 2D-RFCMANN are derived according to gradient descent method with varied learning-rate coefficients for network parameters to assure the convergence of the tracking error. Finally, some experimental tests are realized to validate the effectiveness of the proposed control scheme.

#### 1. Introduction

The cerebellar model articulation neural network (CMANN) is a lookup-table based neural network which can be thought of as a learning mechanism that simulates mammalian cerebellum. Basically, the input space of the CMANN is divided into discrete states, which are called elements. Moreover, several elements construct a block and an overlapped area formed by quantized blocks is called hypercube. Different hypercubes can be obtained by shifting each block a small interval. Furthermore, the CMANN assigns each hypercube to a memory cell for storing learning information. To obtain the output of the CMANN for given input data, only the hypercubes covering the input data will be activated and effectively contribute to the corresponding output. The basic concept of the CMANN is to store learned data into hypercubes and the data can easily be recalled with less occupied storage space. Consequently, the behavior of storing the learned information in the CMANN is similar to the action of the cerebellar in mammalian. Advantages of the CMANN include fast learning, good localization capability, and simplicity of computation [1]. However, in the general CMANN, there are only two kinds of states in the mapped area of the state space, that is, 1 and 0. This directly affects the generalization capability of CMANN and may result in discontinuous output [2, 3]. In order to acquire the derivative information of input and output variables, a CMANN with a differentiable Gaussian receptive field basis function has been proposed in [4]. Due to the lack of time-delayed feedback in its inherent feedforward neural network structure, the application domain of CMANN is limited to static mapping problems. In addition, a recurrent fuzzy CMANN (RFCMANN) with Gaussian receptive field basis function was proposed in [5] to overcome the aforementioned disadvantage of CMANN. In this study, a two-dimensional RFCMANN (2D-RFCMANN) with Gaussian basis function control scheme is proposed for the active and reactive power control of the single-stage three-phase grid-connected photovoltaic (PV) system. The first dimension of the 2D-RFCMANN is for the tracking error input and the second dimension of the 2D-RFCMANN is for the derivative of the tracking error input.

The penetration of renewable energy sources has increased significantly because of raising awareness of the environmental issues and the deregulation of the electric power distribution industry. Among these renewable energy sources, PV system is one of the main candidates to play an important role in the future of power distribution. However, the connection of a large number of distributed energy sources (DESs) to the electrical power grid can cause instability when electrical disturbances appear in the grid. One of the most challenging disturbances is the voltage sags [6, 7]. The sources of the voltage sags are various such as short circuits between phases or phases and ground, lightning, and overloads. To minimize the negative effects of a large number of DESs on the reliability of the power grids, different countries have defined various low voltage ride through (LVRT) regulations for DESs. E.ON has defined the voltage profile of the low voltage ride through (LVRT) at the grid connection for type-2 generating plants [8]. Grid regulations require that the DESs have to inject reactive current to support grid stability during grid faults [8], in which 1 pu of reactive current is required as the depth of the sag reaches 50%. In this study, the E.ON LVRT regulation is adopted to determine the ratio of the required reactive current during grid faults.

The asymmetrical grid voltage sags will worsen the performance of the grid-connected inverter (GCI) especially the degradation of the power quality caused by the power ripple and the increase of current harmonic distortion. In order to overcome these disadvantages, some control methods based on symmetric components have been proposed to deal with unbalanced grid faults [9–13]. A control scheme which can adjust the characteristics of power quality by using two adjustable control parameters was developed in [9]. In [10], a selective active power quality control based on various current reference generating methods was proposed. Reference [11] has proposed a control method to minimize the peak value of the injection current during the voltage sag. In [12], a positive and negative sequences current injection method was proposed to provide the active power and reactive current without exceeding the current limit of inverter. The flexible reactive current injection scheme introduced in [13] provides voltage support during voltage sag by injecting reactive current via both positive and negative sequences. However, the aforementioned control methods principally focus on the power control of grid-connected inverter, whereas the fulfillment of the LVRT requirement and the control of the power flow between the dc-link bus and ac grid during grid faults are not well explored. In this study, the 2D-RFCMANN control scheme is proposed to control not only the output voltage of the PV panel via the maximum power point tracking (MPPT) in normal condition, but also the reactive power based on LVRT requirement with considering the power balance between dc-link bus and ac grid during grid faults.

In this study, the analysis of the instantaneous active/reactive output power of the three-phase GCI during grid faults is carried out by using symmetric component method. The MPPT control of PV system is also considered to ensure the power balance during grid faults. Subsequently, the dual second-order generalized integrator phase-locked loop (DSOGI-PLL) method [10, 14, 15] is adopted for synchronizing the GCI to the grid. Moreover, a new formula which is based on the positive sequence voltage is derived to evaluate the percentage of voltage sag and to determine the injected currents during the grid faults. Furthermore, the network structure and online training algorithms of network parameters of the proposed 2D-RFCMANN for the active and reactive power control are discussed in detail. Furthermore, two types of grid faults are tested and discussed to evaluate the performance of the proposed control scheme. The paper is organized as follows. In Section 2, the grid-connected PV system is introduced. In Section 3, the instantaneous active/reactive power formulation, DSOGI-PLL synchronization mechanism, and the control strategy of the grid-connected PV system during grid faults are discussed. The 2D-RFCMANN is presented in Section 4. In Section 5, the features of the proposed control scheme are examined by some experimental results. Finally, Section 6 draws the conclusions.

#### 2. System Description

Figure 1 shows the block diagrams of the proposed single-stage three-phase grid-connected PV system. The behavior of the PV panel is emulated by a Chroma 62100H-600S PV panel simulator operating at 203 Vdc output voltage with maximum output power 1 kW. The power grid was emulated by three programmable 2 kVA single-phase KIKUSUI PCR2000LE AC power supplies feeding a Y-connected 100 /phase resistive load and connected to the GCI with a 3 kVA Y- transformer. The 2 kVA three-phase GCI, with output 110 Vrms line to line voltage, is connected to the PV panel simulator through a dc-link capacitor and connected to the Y- transformer with three 10 mH/phase coupling inductors. The switching frequency of the three-phase GCI is set at 10 kHz. Table 1 shows the parameters of the experimental setup. Moreover, the control system is installed on a desktop personal computer (PC) which includes an A/D converter card (PCI-1716), a PC-based control system, and a D/A converter card (MRC-6810). In the PC-based control system, the SIMULINK control package is adopted for the implementation of the proposed algorithms of the MPPT and LVRT control. The analog input signals include output voltage and current of the PV panel simulator, the output line voltages , , and , and line currents , , and of the GCI. The control output of the PC-based control system is the three-phase reference currents of GCI , , and . Furthermore, the emulated faults occur at the point of common coupling (PCC).

In Figure 1, the output currents of the GCI are controlled by the maximum power point (MPP) voltage in the normal condition. The adopted MPPT method is the incremental conductance (IC) method with a fixed incremental step size set for the voltage [16]. Moreover, if a grid fault occurs at the PCC, the fault voltages at the GCI terminals appear differently from the fault voltages at the PCC [6, 7]. Furthermore, in the normal condition of the grid, the GCI regulates the active power injecting into the grid and maintains the voltage so as to track the MPP for maximizing the energy extraction. On the other hand, in the fault conditions of the grid, the reactive power control can mitigate the effects of voltage sag by injecting additional reactive current to the grid to ride through the perturbation and support the grid voltage. In addition, the injecting current has to meet the LVRT regulations, without going beyond the GCI current limit simultaneously.

#### 3. Reactive Power Control Strategy during Grid Faults

##### 3.1. Instantaneous Power Formulation under Grid Fault Conditions

For a three-phase GCI, the instantaneous active and reactive power outputs at the grid connection point are described as follows [17]:
where , , and represent the three phase voltages and currents of the GCI, respectively; and represent the instantaneous active and reactive power outputs of the GCI, respectively. In (1) and (2), and are orthogonal only when the three voltages of are balanced. The asymmetric three-phase voltage and current vector are represented as
where and represent the positive and negative sequences of** v**, respectively; and represent the positive and negative sequences of** i**, respectively. Note that the zero sequence of** v **and** i** does not exist due to the three-phase three-wire AC system configuration of the GCI. As a result, and in (1) and (2) can be rewritten as
where the corresponding orthogonal vectors in (5) can be derived by using the matrix transformation in (2). Thus, the instantaneous active and reactive output power of the GCI can be obtained by using (4) and (5) under both normal and abnormal conditions of the grid.

##### 3.2. Grid Synchronization of GCI

The active and reactive power control of the GCI is based on a synchronous reference frame (SRF) with proportional-integral (PI) or 2D-RFCMANN controllers. Moreover, it is important to detect the phase angle of the voltage at PCC in order to synchronize with the grid during grid faults. At present, the most popular PLL structure is the synchronous reference frame PLL (SRF-PLL) [14]. The input three-phase voltages of the SRF-PLL can be given by where and are the amplitude of the positive and negative sequence of the three-phase voltage; , and is the utility angular frequency. In the SRF-PLL, first, the Clarke transformation is applied to transfer the three-phase voltages from the natural reference frame to the stationary reference frame [14, 15]. and are the voltages of , axis which can be expressed as Then, (7) can be expressed in the synchronous reference frame by using Park transformation in the following [15]: where , are the voltages of the , axis; , and is the angular frequency of the synchronous reference frame. Suppose that is close to , then is approximately equal to , and is approximately equal to . The negative sequence component voltage of can be filtered by using a properly designed low pass filter. The difference between and will approach zero due to the closed-loop control of that makes approach zero. Therefore, with in steady state, the phase angle detected by the SRF-PLL will be in phase with the input voltage vector. However, the dynamic behavior of the SRF-PLL will become poor under asymmetric grid faults [15]. In order to overcome this disadvantage of SRF-PLL, a dual second-order generalized integrator PLL (DSOGI-PLL) based on the extraction of the positive sequence has been proposed in [10, 14, 15]. Define as the positive sequence voltage vector on the reference frame and it can be calculated as [10] where is a phase-shift operator in the time domain; and are the filtered voltages of and , respectively. The schematic diagram of DSOGI-PLL is shown in Figure 2, where the dual second-order generalized integrator (SOGI) in Figure 2 is used to generate , , and in-quadrature voltages and of and . Furthermore, and , which are depicted in Figure 2, are the positive sequence voltages on the synchronous reference frame. In addition, a voltage-controlled oscillator (VCO) is employed to convert the input voltage into . Finally, is successfully obtained by DSOGI-PLL by using the positive sequence voltages.

According to [15], the instantaneous active and reactive power outputs of the GCI can be calculated by using the voltage and current in the synchronous reference frame as follows: where and are the currents of the , axis. As mentioned above, an adequate tracking of the grid angular position of the SRF-PLL will make approach zero. Thus, the active and reactive power outputs of the GCI can be represented as Accordingly, and can be regulated by controlling and .

##### 3.3. Reactive and Active Power Control during Grid Faults

The positive sequence voltages of the three-phase line voltages , , and can be expressed as follows: where is the rms value of the magnitude of the positive sequence voltage of the three-phase line voltages and the operator is equal to . Since , , and are symmetric and balanced, can be expressed as Note that , , and are equal to the magnitude of the three-phase line voltages when the line voltages are balanced. According to (12), if the line voltages are unbalanced, , , and will be changed with the degree of unbalance. Therefore is adopted to evaluate the depth of voltage sag during grid faults in this study. The required percentage of compensation reactive current during grid faults for LVRT can be expressed as a function of [8]: and is defined as follows: where is the base value of voltage which equals 110 Vrms in this study. The grid regulation by E.ON demands reactive current injection during fault condition to support the grid stability when is greater than 0.1 pu. Moreover, the apparent output power of GCI will be changed due to the voltage sag at PCC under grid faults. Therefore, it is necessary to calculate the maximum allowable apparent power in order to decide the maximum injection reactive and active power. The maximum apparent power can be expressed as where is the rms value of the GCI current limit and , , and are the rms values of , , and , respectively. Thus, the required injection reactive power and the allowable maximum output active power of GCI as follows can be derived from (14) and (16):

A PV system should inject a certain reactive power in order to support the voltage sag at the PCC according to the LVRT requirements during grid faults. If is greater than , then the active power reference is set to equal the value of before the grid faults. Therefore, all the power generated by the PV panel can be totally delivered into the grid during grid faults and the MPPT control can be operated to keep tracking the MPP. On the other hand, if is smaller than , then is set to . Thus the output active power will be controlled to follow , and the GCI gives up tracking the MPP. Consequently, the power balance between and can be held during grid faults.

#### 4. D-RFCMANN Control Scheme

##### 4.1. Description of 2D-RFCMANN

The structure of 2D-RFCMANN with 2 inputs and 1 output is depicted in Figure 3(a), which consists of five spaces: input space, fuzzification space, receptive-field space, weight space, and out space. Moreover, one of the fuzzy if-then rules of the proposed 2D-RFCMANN is realized as follows.

**(a) Structure of 2D-RFCMANN**

(b) 2D-RFCMANN with , , and and its configuration of receptive-field space activated by state |

: if is and is , then where is the number of the layers for each input dimension; is the number of blocks for each layer; is the fuzzy set for the first input, th layer, and th block; is the fuzzy set for the second input, th layer, and th block; is the corresponding output weight in the consequent part. Furthermore, the signal propagation and the basic function in each space of the 2D-RFCMANN are introduced as follows.

*Space 1 (input space). *In this space, the input variables and are both divided into regions (which are called “elements”). The number of represents the input resolution. For example, a 2D-RFMANN with is shown as in Figure 3(b).

*Space 2 (fuzzification space). *According to the concept of CMANN, several elements can be accumulated as a block. In this space, each block performs a Gaussian receptive-field basis function. The fuzzification of the corresponding membership function in the 2D-RFCMANN is realized by the Gaussian basis function which can be expressed as
where is the output of the th block of the th layer of the th input variable; and are the mean and standard deviation of the Gaussian function, respectively; represents the recurrent weight of the th block of the th layer for the th input variable internal feedback unit; denotes the time step; and represents the value of through time delay. Note that the memory term stores the past information of the 2D-RFCMANN network and performs the dynamic property. Considering the example shown in Figure 4(b), the 2D-RFCMANN has two input variables with three layers in the 2D-RFCMANN and four blocks in each layer.

*Space 3 (receptive-field space). *In this space, each block has three adjustable parameters, , , and . The multidimensional receptive-field function is defined as
where is associated with the th layer and the th block for input and the th block for input according to the fuzzy rules in (12). The product is computed by using -norm product in the antecedent part. An organization diagram of the 2D-RFCMANN is depicted in Figure 4(b). Areas formed by blocks, referred to as , , and so forth, are called receptive fields (or hypercube). Moreover, the number of receptive fields is equal to in each layer.

*Space 4 (weight space). *Each receptive-field area corresponds to a particular adjustable value in the weight space. The output of the weight memory can be described as
where is the output of the weight memory corresponding to and denotes the connecting weight value of the output associated with the th block of and the th block of in the* i*th layer.

*Space 5 (output space). *The output of the 2D-RFCMANN is the summation of the entire activated weighted field:
where is the product of the output of the weight memory and the output of receptive-field. Moreover, for the active power control and for the reactive power control in the GCI.

An example is demonstrated to show the relation between 2D-RFCMANN and its fuzzy inference system. As shown in Figure 3(b), if the inputs of 2D-RFCMANN fall within the state (3, 5), three reception fields are activated which are , , and . Each hypercube of the reception-field space is subject to the respective weights , , and in the weight space.

##### 4.2. Online Learning Algorithm of 2D-RFCMANN

Since the online learning algorithms of the active and reactive power controls are the same, only the detailed derivation of the reactive power control is discussed in this study. To describe the online learning algorithm of the 2D-RFCMANN for the input command tracking using supervised gradient decent method, first the energy function is defined as where represents the tracking error in the learning process of the 2D-RFCMANN for each discrete time ; and represent the desired output and the current output of the system. Then the backpropagation learning algorithm is described as follows.

*Space 5. *The error term to be propagated is given by
The weights are updated by the amount:
where the factor is the learning rate of memory weights. is updated according to the following equation:

*Space 2. *Applying the chain rule, the update laws for , , and are
where the factors , , and are the learning rates. , , and are updated according to the following equation:
Moreover, the exact calculation of the Jacobian of the system, , is difficult to be determined due to the unknown dynamics of the single-stage three-phase grid-connected PV system. To overcome this problem, a delta adaptation law is adopted as follows [18]:
where and represent the first derivatives of and , respectively. Furthermore, the selection of the values of the learning rate parameters of the 2D-RFCMANN has a significant effect on the network performance. Therefore, in order to train the 2D-RFCMANN effectively, the varied learning rates, which guarantee the convergence of the tracking errors based on the analysis of a discrete-type Lyapunov function, are derived as follows:
where is a positive constant. The derivation is similar to the one in [18] and is omitted in this study.

#### 5. Experimental Results

The schematic diagram of the PC-based control system is depicted in Figure 4. The is obtained by DSOGI-PLL block. The and are obtained by the instantaneous active and reactive power calculation block shown in (4) and (5). The calculation block determines the value of and the grid fault signal is activated, while is greater than 0.1. The and are obtained by the active and reactive power references calculation block. The first PI controller (PI1) or the first 2D-RFCMANN (2D-RFCMANN1) is adopted for the control loop of or to regulate the control output . The control loop of introduces the second PI controller (PI2) or the second 2D-RFCMANN controller (2D-RFCMANN2) to regulate the control output . Moreover, since the single-stage three-phase grid-connected PV system is a nonlinear system with uncertainties, the plant model is difficult to obtain by mathematical analysis or experimentation. Therefore, the parameters of PI controllers are tuned by trial and error in this study to achieve the best control performance considering the stability requirement. The transfer functions of these PI controllers which are adopted in this paper can be generally expressed as , , where is the proportional gain and is the integral gain. The gains of these PI controllers , , , and are tuned to be 1.2, 1, 1, and 1.2. Furthermore, both PI and 2D-RFCMANN controllers are implemented in the experimentation for the comparison of the control performance. In order to measure the control performance of the mentioned controllers, the average tracking error , the maximum tracking error , and the standard deviation of the tracking error for the reference tracking are defined as follows [18]: where is the th value of the reference () and is the th value of the response (). The comparison of the control performance can be demonstrated using the maximum tracking error and the average tracking error. The oscillation of the reference tracking can be measured by the standard deviation of the tracking error.

The emulated irradiance varying from 600 W/m^{2} to 300 W/m^{2} has been programmed in the PV panel simulator to evaluate the MPPT performance at normal state. The measured , , , , and one of the line currents of GCI are shown in Figure 5. In the beginning, , , , and are 610 W, 546 W, 203.8 V, and 3.06 A, respectively, and the amplitude of is 4.1 A. At s, the irradiance changes abruptly to 300 W/m^{2}, and drops to 311 W. Then, drops to 267 W; drops to 201.6 V; decreases to 1.6 A; the amplitude of reduces to 2.3 A. From the experimental results, when the 2D-RFCMANN1 controller is used, the settling time of is about 0.3 s. Meanwhile, the total harmonic distortion (THD) of the current of the GSI is about 4.5%. The THD of the GCI is smaller than the THD regulation of IEEE 1547 standard. Therefore, the performances of the IC MPPT method satisfy the requirements of fast response and low THD for the proposed single-stage three-phase grid-connected PV system operating under normal condition.

Two cases of emulated grid faults with different voltage sags have been programmed in the ac power supplies to evaluate the behavior of the proposed PV system. The following grid faults at PCC are selected in this study: Case 1 two-phase ground fault with 0.2 pu voltage sag (sag type E) and Case 2 two-phase ground fault with 0.8 pu voltage sag (sag type E). In the experimentation, first, the grid line voltages are balanced and set at 1.0 pu. Then, at time s, the sag happens and the reactive and active power controls are enabled.

*Case 1. *In this case, the measured , , , , and and the measured three-phase currents are shown in Figure 6(a) for the PI controller and Figure 6(b) for the 2D-RFCMANN controller, respectively. In the beginning, is 608 W and is 538 W due to the power losses in GCI. Moreover, the amplitude of the inverter output currents is 4 A. At s, a two-phase ground fault occurs with 0.2 pu voltage sag on grid side and is 0.864 pu (95 Vrms) according to (13), where 1.0 pu is equal to 110 Vrms. and are equal to 0.136 pu and 27.2% according to (15) and (14), respectively. By using (16) and (17), and are determined to be 794 W and 224.6 var, respectively. Since is greater than , is set to and eventually both and are unchanged. Thus, rises to 224 var during the grid faults. Furthermore, the injection currents are approximately balanced and limited under 7.1 A peak value. From the experimental results, when the PI controllers are used, the settling time of is about 1.2 s. On the other hand, the settling time of is reduced to 0.3 s by using the 2D-RFCMANN controllers.

**(a)**

**(b)**

*Case 2. *In this case, the measured , , , , , and and the measured three-phase currents are shown in Figure 7(a) for the PI controllers and Figure 7(b) for the 2D-RFCMANN controllers, respectively. At s, two-phase ground fault occurs with 0.8 pu voltage sag on grid side. Thus, and are equal to 0.462 pu and 100%, respectively, and is set to zero and the GCI gives up tracking the MPP. As a result, and drop to 92 W and 2 W, respectively. Moreover, the GCI injects 493 var reactive power into the grid due to the 100% reactive current compensation. When either PI controller or 2D-RFCMANN controller is used, rises from 202 V to 230 V and drops from 3.1 A to 0.37 A due to the injected active power being set to zero. From the experimental results, when the PI controllers are used, the settling time of and is about 1.0 s and 0.6 s, respectively. On the other hand, the settling time of and is about 0.4s and 0.1 s by using the 2D-RFCMANN controllers. Furthermore, the maximum injection currents are still limited within the maximum current limit.

**(a)**

**(b)**

From the experimental results of the previous two cases, it can be stated that the proposed single-stage three-phase grid-connected PV system has clearly demonstrated its capabilities. These include the ability of dealing with voltage sags with satisfying the LVRT regulation and the ability of avoiding over current of the GCI. The performance measuring of the PI controller and 2D-RFCMANN controller for the command tracking is shown in Figure 8, using (31). From the resulted performance measuring of PI and 2D-RFCMANN controllers, the performances of 2D-RFCMANN controllers are superior to PI controllers.

**(a)**

**(b)**

#### 6. Conclusions

This study has proposed an intelligent active and reactive power control scheme using 2D-RFCMANN for single-stage three-phase grid-connected PV system operating during grid faults. The mathematical analysis of the injected reactive power during the grid faults has been introduced. A new formula for evaluating the percentage of voltage sag has been proposed. The proposed formula, which is based on positive sequence voltage, is used to determine the ratio of the injected reactive current during grid faults for committing to the LVRT regulation. Moreover, the network structure, online learning algorithms, and learning rates of the proposed 2D-RFCMANN have been discussed in detail. The analyses of the experimental results indicate that the proposed scheme is capable of providing reactive current support for a grid-connected PV system during grid faults. Furthermore, the control performance of the intelligent control scheme using 2D-RFCMANN is better than traditional control scheme using PI due to its online training of network parameters and the capability of parallel processing.

The contributions of this study are (1) the development of a single-stage three-phase grid-connected PV system with the compensation of reactive current during grid faults; (2) the design of the 2D-RFCMANN control scheme and its online learning algorithms with convergence analysis; and (3) the development of the 2D-RFCMANN controller in the reactive power and PV panel voltage control of a single-stage three-phase grid-connected PV system according to the LVRT regulation with the limitation of the maximum current output of GCI during grid faults.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

The authors would like to acknowledge the financial support of the National Science Council of Taiwan, through Grant no. NSC 103-3113-P-008-001.