C:\Users\raoh\Desktop\Paper 3.xps The Journal of Engineering Research Vol. 12, No. 1 (2015) 32-43 Investigations on Hysteresis-Based Current Control Techniques for Grid Connected Photovoltaic Systems V. Rajini Department of Electrical and Electronics Engineering, Sri Sivasubramaniya Nadar (SSN) College of Engineering, Kalavakkam, India Received 21 July 2014; accepted 16 December 2014 Abstract: This paper focuses on delivering an effective performance qualification (PQ) control strategy for a grid connected photovoltaic (PV) system. A PV-based distributed generation system with a capacity of 10 kilo volt amps (kVA) was connected to the utility side of a power grid operating at 415 V. The power grid was modelled with a capacity of 100 kVA, 11 kV, and 100 km transmission line length. The PQ control strategy was implemented with three hysteresis-based current control techniques. To overcome the drawbacks of conventional hysteresis-control techniques, adaptive tuning of the hysteresis band was carried out, and a vector-based hysteresis current control is proposed to improve the results. Keywords: PQ control, PV-based distribution grid system, Vector-based hysteresis current control, Adaptive hysteresis current control, Conventional hysteresis current control. PQ 10 41510011 100PQ PQ Corresponding author’s e-mail: rajiniv@ssn.edu.in V. Rajini 33 1. Introduction The growing demand for power and an impending energy crisis have revealed the acute need for the expansion of clean energy to better meet the world’s heightened power needs. Integration of solar power into the distribution grid (DG) is one of the greatest concerns in engineering. Efficient harnessing of this energy by a photovoltaic (PV) system and directly feeding it into the grid has always been the key solution to this concern. But such a system’s low efficiency and poor controllability have remained the major drawbacks of a PV-based DG system (Blaabjerg et al. 2004). For a grid-connected PV system, both a voltage source pulse-width modulation (PWM) inverter and a system of closed loop power control are required. Three-phase PWM voltage source inverters have been used in high performance ac motor drives, high power factor converters, active filters, etc. (Rahim and Mekhilef 2002). In most of these applications, the magnitude of alternating current (AC) voltage is a function of the voltage across the direct current (DC) link capacitor engaged in the input side of the inverter (Kwon et al. 1998; Mohseni and Islam 2010; Vahedi et al. 2011). Regulated PWM voltage source inverters where the inner current loop determines the control over the injected power are widely employed in grid-connected DG systems in order to achieve independent control over the active and reactive power, and improve load dynamics and power quality (Chitti Babu et al. 2008; Vahedi et al. 2011). The application of dynamic control for power injection is examined in the current study because of the varying demand for power in the utility grid due to the erratic variation in load requirements (Vahedi et al. 2011; Babu et al. 2012). Employing a performance qualification (PQ) control strategy can aid with the varying power requirements of the utility grid, and hence generate reference vectors for the inner current control loop of the voltage source PWM inverter. Effective control of the grid parameters to determine the real and reactive power references is presented and the simulation results analysed. A vector-based hysteresis control with adaptive tuning is proposed in this paper. The basic model of a grid-connected PV system employs three different stages, including the utility grid, the DG source side control, and the three-phase pulse-width modulation voltage- source invertor (PWM-VSI). The utility grid includes a three-phase power source, the distribution line, and a three-phase resistive inductive (RL) load which varies periodically. Figure 1 is a block diagram of a grid-connected PV system. The PV source—the solar array—is modelled mathematically through the basic parameters of a PV cell (Rahim and Mekhilef 2002). An incremental conductance-based maximum power point tracking (MPPT) algorithm was employed to obtain maximum power from the source at all stages. The output of the solar array was interfaced through a DC-DC boost converter and connected through a dc-link to the inverter. The three-phase PWM-VSI was connected to the utility grid in order to inject a three-phase current which was controlled by a PQ controller. Figure 2 shows the inverter topology connected to the grid with a filter. The solar photovoltaic (SPV) system is modelled using the five-parameter model with equations 1–4. A detailed modelling of the SPV system can be found in references 10 and 11. I exp ∗ 1 1 I I exp (2) 0 3 100 4 2. Performance Qualification (PQ) Strategy The paper aims to deliver a generic PQ Figure 1. Block diagram of a grid-connected photovoltaic (PV) system. Figure 2. Inverter Topology. control technique for any DG source. The time variables of three- phase current and voltage vectors are transformed into space vectors in a rotating reference frame (d-q) (Fig. 3). Hence, the real and reactive powers (P and Q) are controlled in the synchronous reference frame (d-q). In order to determine the phase angle (ɵ), Clark’s transformation was performed to obtain the voltage and current vectors in the stationary frame (α-β). This can be given by the Clark’s transformation matrix for a three-phase system. where CT is the Clark’s transformation matrix. The current and voltage vectors are hence given by Eqns. (7) and (8). Vdc R+jXL R+jXL R+jXL S1 S3 S5 S6S4S2 Va Vb Vc Ia Ib Ic 34 V. Rajini 35 Figure 3. Stationary and rotating frames. III 23 1 12 120 √32 √321√2 1√2 1√2 5 III 6 Park’s transformation was employed to convert the stationary two-phase (iα, iβ) system into a rotating (id, iq) system using the unit vector angle (ɵ) (Brod and Novotny 1985) as follows: cosθ sinθcosθ sinθ 7 Using the same transformation, the voltage vectors (Vd) and Vq were obtained using Eqn. (8). cosθ sinθcosθ sinθ (8) The real and reactive powers are given by Eqns. (9), (10) and (11). 9 (10) (11) The d-axis of the synchronous frame is aligned to rotate with the voltage vector (Vd) and, correspondingly, the value of Vq becomes zero. Hence, the power Eqns. (9) and (10) can be rewritten as Eqns. (12) and (13). (12) 13 The value of Vd corresponds to the voltage magnitude, which is a constant. Hence, precise control of the real and reactive powers to be injected can be achieved by controlling the magnitudes of Id and Iq. Two different constants, KP and KQ, are considered such that the reference current vectors along the direct and quadrature axes are given by Eqns. (14) and (15), ∗ 14 ∗ 15 where Kp and Kq are the real and reactive power coefficients of the PV source, to be tuned in accordance with the real and reactive power demands, respectively. By performing suitable inverse transformations, the above current reference variables are transformed to the a, b, c reference frame (Iaref, Ibref, Icref), and are used in the inner current control loop in order to generate the PWM pulses for the VSI. The PQ control system can be so defined in order to deliver preferred ratios or desired values of real power (Pi) and reactive power (Qi). In the case of a PV-based DG source of maximum capacity (Ppv), the values of Kp and Kq are selected so that the total power conservation stays inviolate. This can be ensured by following Eqn. (16). 16 3. Current Control The control of the inverter current by using the generated three-phase reference wave is employed for the calculated injection of power into the utility grid. I d βq d θ I β I q α θ I s I α Investigations on Hysteresis-Based Current Control Techniques for Grid Connected Photovoltaic Systems 36 3.1 Conventional hysteresis Current Contro- ller A conventional HCC technique is preferred for its dynamic performance and ease in implementation. In this technique, an error signal is generated by direct comparison of the inverter current with the reference current value (Fig. 4). The switching pulses are so generated that the current error is kept within a predefined hysteresis band. This method aids the control of switches in an inverter by asynchronously varying the current through the load so that it tracks a reference current signal. Figure 4. Gate Pulse Generation But one major drawback of this technique is the occurrence of a high switching frequency at a lower modulation index due to a lack of coordination between three-phase hysteresis comparators (Esram and Chapman 2007). Hysteresis band violations of up to twice the bandwidth are also prone to occur. One of the much squared solutions for this problem is the use of a dual hysteresis band wherein there are two bandwidths: δ and , such that < δ. This method provides a modest constructive effect in case of transients and peak overshoots. Figure 5 shows the error band region. 3.2 Adaptive bandwidth hysteresis current control (HCC) Since the switching frequency in the HCC technique depends on the rate of change of the line current, it varies with the inverter current waveform (Ic). This current also varies with the line inductance (L) and the dc link voltage (Vdc). The bandwidth of the hysteresis band determines the tolerance of current error. Thus, Figure 5. The error band region for conven- tional hysteresis current control (HCC). the bandwidth determines the rate at which the inverter current gets switched within the hysteresis band and, hence, the switching frequency (fs). It can be concluded that the relation between the bandwidth and switching frequency is: δα 1 f 17 Even though the increase in switching frequency (fs) can improve the inverter current (Ic), it augments the switching losses in the system, which is undesirable. Therefore, the bandwidth (δ) is selected to balance the switching losses, switching frequency (fs), and inverter current (Ic). In order to adaptively tune the bandwidth (δ), the fact that the switching frequency (fs) is a function of the rate of change of the line current ( ) must be considered. Thus, in a switching condition where the switch on time is ton and the switch off time is toff, the time period (T) and switching frequency (fs) are given by: 1 18 The switching of voltage and the line current (Il) for one of the phases at the switching frequency mentioned in Eqn. (18) is illustrated below: Iref Iact Ie Gate Pulse Iref Iβ Iα Ie δ δ δ V. Rajini 37 Figure 6. Switching of inverter line current within the hysteresis band. Figure 6 is analysed in order to determine the adaptive tuning of the hysteresis bandwidth (Vahedi et al. 2011). Considering the differential rise and fall in line currents, the line inductance (L), phase voltage (Vph) and dc input voltage (Vdc) are as follows: 2 (19) 2 (20) By approximating values in Fig. 7, the equations of the bandwidth are as follows: 2 21 2 22 Thus, by subtracting Eqn. (21) from (22), Eqn. (23) can be derived: 1 23 By adding Eqns. (21) and (22), the following is derived: 4 24 By analysing Eqns. (19), (20), (23), and (24), bandwidth (δ) can be represented as 1 12 4 251 2 3.3 Vector-based Hysteresis Current Control (HCC) The selection of optimum bandwidth (δ) was carried out in the case of the adaptive HCC technique. But this technique carries a prime disadvantage of exhibiting high switching frequencies at the zero crossing points of the reference current. Even though this approach yields an almost fixed switching frequency and minimizes inverter current oscillations, this improved performance is usually obtained at the expense of extra signal processing and control complexity requirements, which compromises the simplicity of conventional HCCs (Esram and Chapman 2007). Moreover, we may observe stability problems and limited transient performance in these adaptive tuning techniques. As an alternate to achieving reduced switching frequencies, vector-based HCC techniques can be used. Optimal vector-based methods systematically select zero and non-zero voltage vectors to follow the reference current vector, thereby significantly reducing the inverter switching frequency. A vector-based HCC that employs two levels of hysteresis comparators of bandwidths δ and δ+ δ, to generate an optimal switching pattern that restricts the current error vector within a hexagonal tolerance region, is discussed. Vdc/ 2 -Vdc/2 2 Investigations on Hysteresis-Based Current Control Techniques for Grid Connected Photovoltaic Systems 38 Table 1: Switching table for vector HCC. Sa1 Sb1 Sc1 Sa2 Sb2 Sc2 000 001 010 100 011 101 110 111 100 V01 V01 V01 V1 V01 V6 V2 V01 110 V00 V00 V3 V1 V00 V00 V2 V00 010 V01 V01 V3 V01 V4 V01 V2 V01 011 V00 V5 V3 V00 V4 V00 V00 V00 001 V01 V5 V01 V01 V4 V6 V01 V01 101 V00 V5 V00 V1 V00 V6 V00 V00 This technique retains the simplicity and ease of implementation showcased by the conventional HCC. Furthermore, it offers a considerably reduced switching frequency and inverter current oscillations, thereby exhibiting improved transient and steady state performance. Figure 7. Error band region for vector hysteresis current control (HCC). In order to understand the theory behind vector-based HCC schemes, an understanding of the basic voltage equations of the three phases in case of a standard inverter topology is necessary (Fig. 2). 13 2 13 2 13 2 (26) Figure 8. Tracking current error. Iref Iact Iβ Iα δ δ +∆δ Ie V2 V3 V4 V5 V6 Vopt V0,V7 , V1 V. Rajini 39 The vector form of the above equations based on Clark’s transformation can be given as 1 27 where Vn is the output voltage vector and eo is the vector notation of back-electromagnetic field (emf) (Kwon et al. 1998; Mohseni and Islam, 2010). The current error can be noted as the difference between the reference and the actual values of current. (28) Using Eqns. (27) and (28), we can denote the derivative of error current as follows: 29 By eliminating the R term in equation (29), the equation for the voltage vector corresponding to the optimal switching state, such that the error current vector Ie is zero, results. 30 Equations (29) and (30) denote the derivative of the current error vector (Fig. 9). This derivative can be written as (31) In the case of steady state conditions, vector- based HCC techniques usually follow control strategies to apply voltage vectors such that the tracking error is kept minimal. In Figure 8, where the optimal voltage vector lies in sector 1, of the six non-zero and two zero vectors (Vectors V0–V7) (Fig. 9). the switching states corresponding to the two zero vectors and the adjacent non zero vectors (V1 and V2) are the best to be applied. In contrast, during transient conditions it may be necessary to switch back the actual current Iact within the hysteresis band as quickly as possible irrespective of the tracking error value. The switching table for this technique is given in Table 1. 4. Simulation Results The PQ control strategy is employed using the three current control techniques, in a three-phase 415 V system, rated 11 kV at the generation end and transmission lines of length 100 km. A resistor–capacitor circuit (RC) filter is designed to eliminate a band of harmonics at the inverter end to improve power quality. Figure 10 represents the phase voltage at the point of common coupling (230 Vrms). To illustrate the effect of the proposed controller, the load power, power supplied by the grid, and that generated by SPOV are varied and the details of variations are given in Table 2. When the real power demand is 8000 Watts (W), a peak current of 29.55 amps (A) is supplied to the load out of which 14.49 A peak comes from the grid. When the real power demand changes to 16,000W as indicated in Table 2, a peak current of 49 A is supplied to the load out of which 43 A comes from the grid. Both these current wave forms are given in Figs. 11(a) and (b). Similarly, the power demand and the power delivered by the main grid are shown in the Figs. 12(a) and (b). The sum of the power delivered by the main grid and DG source equals the power demand. Figure 9. Space vector representation. The load is rated to demand 8000 W of real power and 12,000 volts-amperes reactive (VAR) of reactive power for the time period of 0–0.1 s. An additional load is connected to the system at 0.1 s, so that the demand rises to 16,000 W and 24,000 VAR for the time period 0.1–0.2 s. Hence the PQ controller (section 2) allows the DG system to inject 8,000 W of real power and 6,000 VAR of Investigations on Hysteresis-Based Current Control Techniques for Grid Connected Photovoltaic Systems 40 reactive power for 0–0.1 s (KP = 1; KQ = 0.5). Once the demand rises, the controller allows the system to deliver its total capacity (10,000 W) as real power (KP = 0.625; KQ = 0). Table 2. Power Allocation. Periods/Parameters 0 – 0.1 0.1- 0.2 PL(W) 8000 16000 QL(VAR) 12000 24000 Pi(W) 6000 10000 Qi(VAR) 6000 0 PG(W) 2000 6000 QG(VAR) 6000 24000 KP 1 0.625 KQ 0.5 0 Figure 10. Voltage at the point of common coup ling. Figure 11(a). Load current required. Figure 11(b). Current drawn from the main grid. Figure 12(a). Real and reactive power demand of the load Figure 12(b). Real and reactive power delivered the main grid. The inverter current waveform and the real and reactive power delivered by the DG source for the three current control techniques are shown in the Figs. 13 and 14, respectively. It can be inferred from the figures that the inverter current and power waveforms of the vector HCC technique are smoother compared to the other two techniques. V. Rajini 41 Figure 13(a). Current delivered by the distributed generation (DG) source—conventional hysteresis current control (HCC). Figure 13(b). Current delivered by the distributed generation (DG) source— Adaptive hysteresis current control (HCC). Figure 13(c). Current delivered by the distributed generation (DG) source—vector hysteresis current control (HCC). Figure 15(a) seems to resemble Fig. 15(b). Figure 16 shows the error bandwidth in the stationary frame. Figure 15 shows that the error band region is more confined and the magnitude of the error vector is restricted in the case of a vector-based technique. Figure 15(b) shows a tuned bandwidth resembling Fig. 7. Figure 14(a). Real and reactive power delivered by the distributed generation (DG) source— conventional hysteresis current control (HCC). Figure 14(b). Real and reactive power delivered by the distributed generation (DG) source—adaptive hysteresis current control (HCC). Figure 14(c). Real and reactive power delivered by the distributed generation (DG) source—vector hysteresis current control (HCC). Investigations on Hysteresis-Based Current Control Techniques for Grid Connected Photovoltaic Systems 40 5. Conclusions An effective real and reactive power control strategy for the distribution side power injection of a PV DG source has been discussed. This PQ control strategy (a) (b) (c) Figure 15(a). Current error in dq plane— conventional hysteresis current control (HCC), (b) adaptive HCC and (c) vector HCC. Figure 16. Adaptive band tuning. Table 3: Comparison of simulation results Parameters Switching loss (Watts) THD of inverter current (%) Conventional 170 5.4 Adaptive HCC 130 3.2 Vector HCC 60 2.1 was executed through a conventional HCC technique. Having observed the disadvantages in implementing a conventional HCC technique, an adaptively tuned HCC and a vector-based HCC were employed in the system, Their simulation results were then compared and the system was analysed for the three current control techniques. Based on the observations carried out on parameters such as total harmonic distortion and the switching loss of the system, the vector-based control technique was found to be more appropriate and effective. References Babu D, KN, Ramaprabha R, Rajini V (2012), Mathematical modeling and simulation of grid connected solar photovoltaic system. International Journal of Electrical and Electronics Engineering 20(5):73–77. Blaabjerg F, Chen Z, Kjaer S (2004), Power electronics as efficient interface in dispersed power generation systems. IEEE Transactions on Power Electronics 19(5):1184–1194. Brod DM, Novotny DW (1985), Current control techniques of VSI-PWM inverter. IEEE Transactions on Industrial Applications 562– 570. Chitti Babu B, Mohapatra M, Jena M, Naik A (2008), Dynamic performance of adaptive hysteresis current controller for mains- connected inverter system. IEEE Transactions on Power Electronics 22(3):197– 306. Esram T, Chapman PL (2007), Comparison of photovoltaic array maximum power point tracking techniques. IEEE Transactions on Energy Conversion 22(2):439–449. Kwon BH, Kim TW, Youn JH (1998), A novel SVM-based hysteresis current controller. IEEE Transactions on Power Electronics 13(2):297– 307. Mohseni M, Islam SM (2010), A new vector-based hysteresis current control scheme for three- phase PWM voltage-source inverters. IEEE Transactions on Power Electronics 25(9). Rahim NA, Mekhilef S (2002), Implementation of three-phase grid connected inverter for photovoltaic solar power generation system. Proceedings IEEE PowerCon, 1:570–573. 42 V. Rajini 43 Salmi T, Bouzguenda M, Gastli A, Masmoudi A (2012), MATLAB/Simulink-based modelling of solar photovoltaic cell. International Journal of Renewable Energy Research 2(2). Vahedi H, Sheikholeslami A, Bina MT (2011), A novel hysteresis bandwidth (NHB) calculation to fix the switching frequency employed in active power filter. IEEE Applied Power Electronics Colloquium (IAPEC), 156. Appendix Transmission line parameters: Resistance: 0.2 ohms/Km Inductance: 0.5mH/Km Transformer rating:11KV/415V Load parameters: From 0-0.1 s: Real power demand: 8000 W Reactive power demand: 12000 VAR From 0.1-0.2 s: Real power demand: 16000 W Reactive power demand: 24000 VAR DC link voltage (Vdc):700 V Adaptive tuning of hysteresis band: Line inductance (L): 5 mH Switching frequency (fs): 12 KHz