Power Quality Analysis of Hybrid Renewable Energy System
Power Quality Analysis of Hybrid Renewable Energy System
Power Quality Analysis of Hybrid Renewable Energy System
discussions, stats, and author profiles for this publication at: http://www.researchgate.net/publication/276310414
READS
41
3 AUTHORS:
Rinchin Mosobi
Toko Chichi
2 PUBLICATIONS 0 CITATIONS
2 PUBLICATIONS 0 CITATIONS
SEE PROFILE
SEE PROFILE
Sarsing Gao
North Eastern Regional Institute of Science a
10 PUBLICATIONS 21 CITATIONS
SEE PROFILE
Sarsing Gao
2015 The Author(s). This open access article is distributed under a Creative Commons Attribution
(CC-BY) 4.0 license.
Page 1 of 15
Keywords: Renewable energy system; maximum power point tracking; boost converter;
solar photovoltaic; wind energy; self-excited induction generator; STATCOM; total
harmonic distortion
1. Introduction
Unlike the conventional energy sources, the non-conventional energy sources are clean, reliable,
and abundant in nature. The environmental degradation such as pollution, global warming, and
greenhouse gas emissions which are caused by conventional sources of energy and accelerated
by ever-growing industrial activities throughout the world is a concern for all. The current
researches, therefore, lay emphasis on harnessing renewable energy sources (RES) for generating
electricity to supply power especially, to rural consumers where grid connection is not available.
For such locations, decentralized power generation using available dispersed RES is a better and
workable solution. It is shown that power generation from combined RES such as wind and
hydro could make a power system more cost-effective and environmental friendly (Nejad, Radzi,
Kadir, & Hizam, 2012). A technological innovation, however, is needed to utilize these energy
sources to an optimum level and to obtain greater efficiency. The combined system is likely to
improve the generating capacity as well the reliability of the power supply. A feasibility study of
solarwinddiesel hybrid system is carried out in (Lipu, Uddin, & Miah, 2013) with promising
results but for the increased cost of diesel oil and associated environmental pollution caused
due to operation of diesel engine. Thus, the system considered may not be of much help. As
reported, a windmicrohydro hybrid system when operated in grid connected mode using
back-to-back connected power converters may improve the power system reliability (Goel,
Singh, Murthy, & Kishore, 2011). The reliability and cost-effectiveness of integrated renewable
energy system at remote and distant places is evaluated using optimization technique (Bansal,
Khatod, & Saini, 2014) with fairly good results. However, it is found that not many works have
been done on the power quality analysis of integrated renewable energy system using solar,
wind, and microhydro energy sources. This system could be very successful especially in the
sub-tropical region where there is sufficient rainfall.
A solar PV arrays with maximum power point tracking (MPPT) can generate electricity on large
scale (Gera, Rai, Parvej, & Soni, 2013). With the development of large wind turbine generators and
advancement in solid-state devices, wind energy system (WES) has become one of the most viable
options of generating energy. However, due to the intermittent nature of wind speed, it calls for
tracking of maximum power point of operation (Villalva, Gazoli, & Ruppert, 2009; Villalva, Gazoli, &
Filho, 2009). Power generation from mini- or microhydro is a well-proven technology (International
Renewable Energy Agency (IRENA), 2012).
In this paper, a HRES comprising photovoltaic (PV) system, wind-turbine generator system, and a
microhydro system is proposed. The proposed system is a combination of 15 2 PV arrays with
variable temperature and irradiance, a variable-speed wind-turbine coupled to self -excited induction
generator (SEIG) and a constant power microhydro system coupled to SEIG. The performances of the
hybrid system under various loads namely, induction motor (IM) load, resistive (R), and R-L load are
presented. The reactive power compensation using static synchronous compensator (STATCOM) is
presented. The STATCOM is observed to be working well thereby regulating the terminal voltage and
reducing the harmonic contents of the system voltage and current.
2. Solar PV system
2.1. System description
A PV cell is basically a semiconductor diode whose p-n junction is exposed to light. When the irradiance falls on the surface of the solar cell, it absorbs photons and excites electron-hole pair separation to produce electromotive force.
Page 2 of 15
Figure 1 shows the equivalent circuit of the ideal PV cell including series (Rs) and parallel (Rp) resistances which actually represents the contact and junction resistances. From the basic equation of
semiconductors, the IV characteristics of the ideal PV cell may be written as
)
]
qV
I = Ipv,cell Io,cell exp
1
aKT
(1)
Id
where Ipv,cell is the current generated by the incident light, Id is the diode current, Io,cell is the reverse
saturation current or leakage current of the diode, q is the electron charge, k is Boltzmann constant,
T (in Kelvin) is the temperature of the p-n junction, and a is the diode ideality factor. Figure 2 shows
the ideal IV characteristics of a solar PV module.
[
(
)
]
V + RS I
V + Rs I
I = Ipv Io exp
1
Vt a
Rp
(2)
where, Ipv and Io are the PV and saturation currents of the array and Vt =
Ns KT
q
of the array with Ns connected in series. If the array is composed of Npp parallel connection of
modules, the PV and saturation currents may be expressed as Ipv = Ipv,cell Npp, Io = Io,cell Npp. In
Equation 2, Rs is the equivalent series resistance of the array and Rp is the equivalent parallel
resistance. Rs and Rp are expressed as
Nss
Npp
Rs and
Nss
Npp
IV curve shown in Figure 3. This curve depends on internal characteristics of the device such as Rs
and Rp and on external influences such as irradiance level and temperature. The solar PV array in this
paper is modeled for Npp=15 and Nss=2. The light-generated current of the PV cell depends linearly
on solar irradiation and is influenced by temperature as,
(
)
G
Ipv = Ipv,n + KI T
Gn
(3)
Figure 2. IV characteristics of
the PV cell.
Page 3 of 15
Figure 3. IV characteristics of
the PV array.
where, Ipv,n is the light-generated current at the nominal condition, T=T Tn, where T and Tn are the
actual and nominal temperatures, G (in W/m) is the irradiance on the device surface and Gn is the
nominal irradiance.
The diode saturation current Io and its dependence on the temperature may be expressed as
Io =
exp
Isc,n + KI T
( V + K T )
oc,n
aVt
(4)
where, KI is current coefficient, KV is voltage coefficient, Isc,n is short circuit current at nominal
condition, Voc,n is open circuit voltage at nominal condition, a is ideality factor, and Vt is thermal
voltage.
Figures 4 and 5 show the simulated values of IV and PV characteristics of the PV array considered
in this work. From Figure 4, it is observed that with decreasing irradiance, there is a negligible change
in the maximum output voltage and there is a marginal decrement in output current and a noticeable
drop in the maximum power point (MPP).
Figure 4. (a) IV and (b) PV
curve of a PV array with
variable irradiance.
Page 4 of 15
Figure 5 shows that for a fixed irradiance value as the cell temperature decreases, a small
decrement in the output current occurs with an increment in output voltage and corresponding
increase in MPP.
Table 1 shows the parameters of the PV array obtained from datasheet of Kyocera KC200GT at
standard test condition (STC) of 25C and 1,000W/m. The KC Series solar modules of Kyocera have
multicrystalline silicon solar cell with an efficiency of over 16%. The same PV modules have been
considered in this study (Kyocera KC200GT, xxxx).
KC200GT
32.9V
Isc
8.21A
KV
0.1230V/K
KI
0.0032A/K
Pmax
200.143W
Rp
415.405
Rs
0.221
1.3
Ns
54
Page 5 of 15
where, Voc = open circuit voltage, Isc = short circuit current, KV = temperature coefficient of Voc,
KI = temperature coefficient of Isc, Pmax = maximum power, Rs = series resistance, Rp = parallel
resistance, a=ideality factor and Ns=number of cells in series. Combining Equations 3 and 4, a new
expression is obtained as
Nss
Nss
V
+
R
V
+
R
s
s
Npp
Npp
I 1
I = Ipv Npp Io Npp exp
I
Vt aNss
N
Rp N ss
pp
(5)
where, Ipv is PV current of the array, Io is diode saturation current, a is ideality factor, Vt is thermal
voltage, I is output current and V is output voltage.
The power-conditioning system includes a boost converter with MPPT, a three-phase CC-VSI, and an
LC filter. The input to boost converter is an unregulated DC voltage obtained directly from the PV
array and therefore, it is likely to fluctuate due to variations in the incident solar radiation, also called
irradiance and temperature. Thus, the average DC output voltage must be regulated to get the
desired constant value throughout its operation. The MPPT used is based on perturbation and
observation (P & O) technique (Kondawar & Vaidya, 2012). For a PV array with operating voltage as
dP
= 0 with its sign defined by
V and current as I, the power is P=VI. At the maximum power point,
dV
1 dP
dI
I
=
+
V dV
dV V
(6)
Figure 6 shows the flowchart of P & O algorithm. The MPPT feeds the desired PV array voltage to the
boost converter through change in duty cycle as a gating signal (Villalva et al., 2009). A change in
applied voltage and current results in changed value of power. If the power increases, voltage and
current increases in the same direction and vice versa.
The boost converter shown in Figure 7 works in two modes:
Mode 1 (0 < t ton):
Page 6 of 15
Id
Vin
Switch
Vi
Pulses
Mode 1 begins when IGBTs is switched on at t=0 and terminates at t=ton. The inductor current IL(t)
greater than zero and ramp up linearly. The inductor voltage is Vi.
(
)
Vi ton + Vi Vo toff = 0
(7)
where, Vi is the input voltage, Vo is the average output voltage, ton is the switching-on time of the
IGBT, and toff is the switching-off time of the IGBT.
Dividing both sides by Ts and rearranging the terms yield
T
Vo
1
= s =
toff
Vi
1D
(8)
Rs
+ jXls + ZL + Zr
F
where, ZL =
jXC RL
(F 2 RL jFXC )
Is = 0
and Zr =
(9)
jXm Rr Xm Xlr (Fv)
j (Xm +Xlr )(Fv)+Rr
Figure 8. Configuration of
three-phase VSI with LC filter.
+
Q1
Q3
Q5
AC Bus
D1
D3
D5
Vdc
Q4
Q6
D4
FILTER
Q2
D6
D2
Rs
+ jXls + ZL + Zr = 0
F
(10)
where, Vg is the air-gap voltage, Vph is the phase voltage, v is the speed, and F is the frequency.
For any speed v, Equation 10 may be solved for Xm and F using NewtonRaphson iteration method
which is ideal for finding solution of non-linear equations or any other known methods of solving
non-linear equations with the starting values being taken as the unsaturated value of Xm, and F
equal to v. Using these values of Xm and F, the total performance of the machine can be evaluated in
conjunction with the measured variation of Xm with Vg/F for the generator being considered.
In wind energy conversion system, a variable speed variable frequency scheme is proposed due to its
ability of efficient energy conversion. The system is operated over a wide range of speed without
employing pitch control mechanism for the wind turbine. Voltage build-up is achieved by connecting
excitation capacitors across the stator terminals. According to wind speed cube law, the power,
P developed by the wind turbine of blade diameter, d at a wind speed, s is given by (Raina & Malik, 1983).
P=
1
d2 Cp s3
8
(11)
where, is the density of air and Cp is the power coefficient which is the ratio of shaft power to wind
power.
The microhydro system is modeled with a hydro turbine driven SEIG. The input to the turbine is
from flowing stream with proper civil works in order to maintain constant flow so that the input
power to the hydro turbine remains constant. Finally, the AC output from wind and microhydro
systems is integrated with the PV system and the combination is brought out to a common AC bus.
the DC link capacitor voltage is monitored. This is an index of the imbalance in real power or indirectly
it is an indication of direct axis component of current, Id. Next, the peak value of the line-to-line
voltage from the integrated system is computed. This is compared with the reference or expected
value. The difference between these two quantities is an indication of the reactive power required by
the system or indirectly this will be the amount of quadrature current, Iq to be supplied by the system.
These two axes reference currents namely Id and Iq are converted into three-phase form by Inverse
Parks transformation. The cos(t) and sin(t) terms needed for Parks transformation are derived
with the help of a phase locked loop which is fed with unit templates of line voltages from the source.
The three-phase reference currents thus obtained are compared with the actual currents from the
integrated system in a hysteresis current controller to yield the firing signals for the six devices in the
STATCOM. Figure 10 shows the schematic diagram of the CC-VSI unit (Gao, Bhuvaneswari, Murthy, &
Kalla, 2014; Goel et al., 2011).
Inverse Parks transformation (dq0-abc transformation) is represented as
Idq0 = TIabc
cos (t)
2
=
sin (t)
3
cos t 2
3
sin t 2
3
2
2
cos t + 2
3 Ia
I
sin t + 2
3
b
2
Ic
(12)
The proposed scheme of integration of solar, wind, and microhydro is modeled and simulated in
Matlab/Simulink environment. The three different RES are connected to a common AC bus on to
which varying loads are connected (Mohammed, 2006; Sharaf & El-Sayed, 2009). The power quality
issues (Anees, 2012) like the system voltage and current regulation and harmonic reductions have
been analyzed with the help of a STATCOM.
Figure 11 shows the hybrid system proposed in this study. The simulink model of the system is
shown in Figure 12. Simulation is carried out initially with the three different sources independently
under no load conditions and eventually with the hybrid system through three-phase AC bus connected
to three-phase IM, R, and RL loads. The results are presented through Figures 1316.
Figure 10. Schematic diagram
of current-controlled voltage
source inverter (CC-VSI).
Vt (ref)
Vab
Vbc
Vca
Vt
Unit Voltage
Template
Generator
PI Controller
Iabc
Pulses to VSI
Current
Regulator
Iq
dqo
abc
Id
PI Controller
+
-
Vdc(ref)
Vdc
Voltage Source Inverter
(VSI)
Page 9 of 15
It is clearly observed that the three different RES generates little over rated voltage of 459V (rms)
or 650 V (peak) under no load conditions. The frequency is maintained at 50 Hz. As seen from
Figure16, on connecting loads namely, IM, R, and R-L loads at 1.0, 1.2 and 1.4s, respectively, there
is a voltage drop in the system. The loads are switched out at 1.6s and the system rated voltage is
regained after sometime. A STATCOM is, therefore, connected to the system at 0.6s and the loads
are connected at the same instants. It is observed that the system voltage and frequency
is maintained at its rated value; in this case the voltage is 415V (rms) or 586V (peak). Figure 17
shows the voltage and current waveforms of the integrated system with STATCOM. It is observed
that the AC bus voltage is regulated with the use of STATCOM though a small amount of transient
current occurs at the time of switching in the STATCOM. The overall performance of the system is
observed to have been improved.
AC BUS
PHOTOVOLTAIC
SYSTEM
PV Array
DC-DC
Boost
Converter
DC-AC
Inverter
Filter
LC
Filter
MICRO HYDRO
SYSTEM
River
IG
Blade
Hydro
Turbine
AC-DC-AC
Conversion
Filter
WIND ENERGY
CONVERSION SYSTEM
STATCOM
LOAD 1
LOAD 2
LOAD 3
Page 10 of 15
(a)
(b)
1000
Voltage (V)
500
0
-500
-1000
0.05
0.1
0.15
Time (sec)
0.2
0.25
0.3
0.05
0.1
0.15
Time (sec)
0.2
0.25
0.3
0.05
0.1
0.15
0.2
0.25
0.3
0.05
0.1
0.15
TIME (SEC)
0.2
0.25
0.3
0.05
0.1
0.2
0.25
0.3
0.05
0.1
0.2
0.25
0.3
Current (A)
100
50
0
-50
-100
500
0
-500
50
CURRENT (A)
VOLTAGE (V)
(a)
0
-50
Voltage (V)
(b)
500
0
-500
Current (A)
0.15
50
0
-50
0
0.15
Time (sec)
Page 11 of 15
500
Voltage (V)
0
-500
0.2
0.21
0.22
0.23
0.24
0.21
0.22
0.23
0.24
0.25
0.26
0.27
0.28
0.29
0.3
0.26
0.27
0.28
0.29
0.3
Current (A)
-5
0.2
500
0
-500
0.9
1.1
1.2
1.3
Time (sec)
1.4
1.5
1.6
1.7
1.1
1.2
1.3
Time (sec)
1.4
1.5
1.6
1.7
Current (A)
20
-20
0.9
1000
Voltage (V)
500
0
-500
-1000
0.6
0.8
1.2
Time (sec)
1.4
1.6
1.8
0.6
0.8
1.2
Time (sec)
1.4
1.6
1.8
50
Current (A)
0.25
Time (sec)
1000
Voltage (V)
-50
The total harmonic distortion (THD) is an important figure of merit used to quantify the level of
harmonics in voltage and current waveform. There exist two different definitions of THD, in the first
one the harmonic content is compared to the fundamental and in the second one to the rms value
of the waveform (Similovitz, 2005).
The THD (compared to fundamental) of the system voltages and currents under load conditions are
analyzed here and they are observed to be well within IEEE Standard 5191992. Figures 18 and 19
show the THD of system voltages and currents under IM, R, and R-L loads.
Page 12 of 15
Page 13 of 15
6. Conclusions
In this work, an attempt has been made to investigate in detail the possibility of HRES as one of the
most effective ways of decentralized power generation. Besides increasing the generating capacity of
the energy sources, the scheme is envisaged to greatly liberate the rural populace of the perennial
energy demands. With the energy sources connected together, the reliability of the power system
increases. The power quality aspects are studied and with the use of STATCOM the voltage and current
regulations improves. The system voltage is well maintained at 415 V (rms) under different load
conditions. The harmonic content of the system voltages and currents is greatly reduced when loads
are connected. The THD of the voltages and currents are of the order of 0.11.3% and 0.070.7% for
IM, R, and RL loads, respectively. The voltage and frequency has been maintained constant under
different dynamic conditions. The integrated system along with the addition of STATCOM can be seen
as a viable option for supplying electricity to far flung areas.
Funding
The authors received no direct funding for this research.
Author details
Rinchin W. Mosobi1
E-mail: [email protected]
Toko Chichi1
E-mail: [email protected]
Sarsing Gao1
E-mail: [email protected]
1
Department of Electrical Engineering, North Eastern Regional
Institute of Science and Technology, Nirjuli, Arunachal
Pradesh 791 109, India.
Citation information
Cite this article as: Power quality analysis of hybrid
renewable energy system, Rinchin W. Mosobi, Toko Chichi &
Sarsing Gao, Cogent Engineering(2015), 2: 1005000.
References
Anees, A. S. (2012). Grid integration of renewable energy
sources: Challenges, issues and possible solutions. In
5th India International Conference on Power Electronics
(IICPE) (pp. 16). Delhi: IEEE.
Bansal, M., Khatod, D. K., & Saini, R. P. (2014, February).
Modeling and optimization of integrated renewable
energy system for a rural site. In International Conference
on Reliability, Optimization and Information Technology
(ICROIT) (pp. 2528). Faridabad.
http://dx.doi.org/10.1109/ICROIT.2014.6798289
Gao, S., Bhuvaneswari, G., Murthy, S. S., & Kalla, U. K.
(2014, March). Efficient voltage regulation scheme for
three-phase self-excited induction generator feeding
single-phase load in remote locations. IET Renewable
Power Generation, 8, 100108.
http://dx.doi.org/10.1049/iet-rpg.2012.0204
Gera, R. K., Rai, H. M., Parvej, Y., & Soni, H. (2013, JanuaryJune).
Renewable energy scenario in India: Opportunities and
challenges. Indian Journal of Electrical And Biomedical
Engineering, 1, 1016.
Goel, P. K., Singh, B., Murthy, S. S., & Kishore, N. (2011, April).
Isolated wind-hydro hybrid system using cage generator
and battery storage. IEEE Transactions on Industrial
Electronics, 58, 11411153.
Hasaneen, B. M., & Mohammed, A. A. E. (2008, March). Design
and simulation of DC/DC boost converter. In IEEE Power
System Conference, MEPCON (pp. 335340). Aswan.
International Renewable Energy Agency (IRENA). (2012,
June). Renewable energy technologies: Cost analysis of
hydropower, 1(3/5). Retrieved from www.irena.org/
Kondawar, S. S., & Vaidya, U. B. (2012, August). A comparison
of two MPPT techniques for pv system in Matlab simulink.
Page 14 of 15
Appendix
PV arrayTemperature(K)Solar
irradiance(W/m2)
Boost converter
L=1mH, C=50F
L=10H, C=350F
LC filter
Wind speed(m/s )
Wind turbine
11, 10, 12, 13, 9, 11, 14Pn=1.5e 6, =0, =8.1, Cp=0.48, base power=1.67e 6VA,
base wind speed=12m/s
SEIG
2015 The Author(s). This open access article is distributed under a Creative Commons Attribution (CC-BY) 4.0 license.
You are free to:
Share copy and redistribute the material in any medium or format
Adapt remix, transform, and build upon the material for any purpose, even commercially.
The licensor cannot revoke these freedoms as long as you follow the license terms.
Under the following terms:
Attribution You must give appropriate credit, provide a link to the license, and indicate if changes were made.
You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.
No additional restrictions
You may not apply legal terms or technological measures that legally restrict others from doing anything the license permits.
Cogent Engineering (ISSN: 2331-1916) is published by Cogent OA, part of Taylor & Francis Group.
Publishing with Cogent OA ensures:
High visibility and discoverability via the Cogent OA website as well as Taylor & Francis Online
Input from, and dialog with, expert editors and editorial boards
Page 15 of 15