![ ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 9, No. 3, June 2019 : 1479 - 1495
1480
unbalance, reactive currents, and neutral current. B. Han et. al and Vinod Khadkikar have investigated UPQC
as one part of active power filter consisting of shunt and series active filters connected in parallel and serves
as superior controller to overcome a number of power quality problems simultaneously. Series component has
been responsible for reducing a number of interference on source side i.e. voltage sag/swell, flicker, unbalanced
voltage, and harmonics. Shunt component has been responsible for addressing a current quality problems i.e.
low power factor, load current harmonics, and unbalanced load. [1, 2]. UPQC based on RE has been
investigated by some researchers. There are two methods used to overcome this problem i.e. using conventional
and artificial intelligence. Shafiuzzaman K.K., et.al have proposed system includes a series inverter, shunt
inverter, and a DG connected to a DC link through a rectifier using PI. The system was capable of increase
source voltage quality i.e. sag and interruption and load current quality, as well as transfer of active power
on/off grid mode [3]. The influence of DG on UPQC performance in reducing sag under conditions of some
phase to ground faults to using DSTATCOM has been implemented by Norshafinash S., et.al. The DG was
effective enough to help UPQC to improve sag [4]. It was connected in series with load resulting better sag
mitigation compared to system without DG. Implementation of UPQC using UVTG method with PI to reduce
sag, swell, voltage/current harmonics has been doneby S. N. Gohil, et.al. Simulation of voltage distortion was
made by adding 5th
and 7th
harmonics at fundamental source voltage, resulting in a reduction of THD source
current and THD load voltage [5].
UPQC supplied by PV panels using boost converter, PI, MPPT P and O, and p-q theory has been
proposed by Yahia Bouzelata at.al [6]. The system was capable of compensate reactive power and reduce
source current/load voltage harmonics, but did not discuss migitation of sag and interuption caused by PV
penetration. Power quality enhancement of sag and source voltage harmonics on grid using UPQC supplied by
PV array connected to DC link using PI compared with FLC has been done by Ramalengswara Rao, et.al.
Combination of UPQC and PV using FLC can improve source voltage THD better than PI [7]. Amirullah et.al
have researched a method for balancing current and line voltage, as a result of DGs of a single phase PV
generator unit in randomly installed at homes through on a three phase four wire 220 kV and 50 Hz distribution
line using BES and three of single phase bidirectional inverter. Both devices was capable of reduce unbalanced
line current/voltage, but both of them were also capable of increase current/voltage harmonics on PCC bus [8].
Power quality migitation of UPQC on microgrid supplied by PV and wind turbine has been implemented
by K S Srikanth et. al. It resulted that PI and FLC was able to improve power quality and reduce distortion in
output power [9]. The UPQC-wind turbine to provide active power to overcome low sag and interruption
voltage to grid has been investigated by H.Toodeji, et.al. The model was used VSC as a rectifier on generator
output and controlled so that maximum power desired can be generated by different speed wind turbines using
PI [10]. The UPQC-wind turbine connected to UPQC DC link was implemented by M. Hosseinpour, et al.
The proposed combination using PI was capable of compensate swell, interruption voltage, and reactive power
both on on/off grid [11]. R.Bhavani, et, al have researched on UPQC controlled by FLC to improve power
quality in a DFIG wind turbine connected grid. FLC can improve power quality i.e. sag voltage and load current
harmonics better than PI [12]. Power quality enhancement on wind turbine and BES with PI connected grid on
PCC bus using UPQC has been implemented by S.RajeshRajan, et al. BES was installed to maintain and
stabilize active power supply under different wind speed [13].
This research will analyze UPQC-BES performance supplied by PV-wind hybrid connected to 3P3W
of 380 volt (L-L) and 50 hertz distribution system.The performance of supply system is compared with two
RE sources i.e. PV and Wind, respectively. BES serves to store excess energy produced by three RE sources
and distribute it to load if necessary, to prevent interruption voltage, and to adjust charging and discharging of
energy in battery. BES is also expected to store excess power produced by three RE combinations and use it as
backup power. FLC is proposed and compared with PI to control variable of DC voltage and DC reference
voltage input to generate reference current source in current hysteresis controller on shunt active filter.
DC voltage controller in shunt active filter and series active filter is used to migitate power quality of load
voltage and source current.Performance of two controllers are used to determine load voltage, source current,
load voltage THD, and source current THD based on IEEE 519. This paper is presented as follow. Section 2
describes proposed method, model of UPQC-BES system supplied by three RE sources i.e. PV, wind,
and PV-wind hybrid, simulation parameters, PV and PMSG wind turbine model, series and shunt active filter,
as well as application of PI and FLC method for proposed model. Section 3 shows results and analysis about
performance of THD analysis on the proposed model of three RE sources connected to DC link of UPQC-BES
system using PI and FLC. In this section, six disturbance scenarios are presented and the results are verified
with Matlab/Simulink. Finally, this paper in concluded in Section 4.](https://image.slidesharecdn.com/0313075-200724025151/85/Matlab-simulink-simulation-of-unified-power-quality-conditioner-battery-energy-storage-system-supplied-by-PV-wind-hybrid-using-fuzzy-logic-controller-2-320.jpg)
![ ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 9, No. 3, June 2019 : 1479 - 1495
1482
2.2. Photovoltaic model
Figure 2 shows the equivalent circuit and V-I characteristic of a solar panel. A solar panel is composed
of several PV cells that have series, parallel, or series-parallel external connections [14].
(a) (b)
Figure 2. Equivalent circuit and V-I characteristic of solar panel
The V-I characteristic of a solar panel is showed in (1):
P
S
t
S
oPV
R
IRV
aV
IRV
III
1exp (1)
where IPV is the photovoltaic current, Io is saturated reverse current, ‘a’ is the ideal diode constant,
Vt=NSKTq-1
is the thermal voltage, NS is the number of series cells, q is the electron charge, K is the Boltzmann
constant, T is the temperature of p–n junction, RS and RP are series and parallel equivalent resistance of the
solar panels. IPV has a linear relation with light intensity and also varies with temperature variations.
Io is dependent on temperature variations. The values of Ipv and Io are calculated as following (2) and (3):
I
G
G
TKII
n
nPVPV )( 1,
(2)
1/)exp( ,
,
tVnOC
InSC
o
aVTKV
TKI
I (3)
In which IPV,n, ISC,n and VOC,n are photovoltaic current, short circuit current and open circuit voltage in
standard conditions (Tn=25 C and Gn=1000 Wm-2
) respectively. KI is the coefficient of short circuit current to
temperature, ∆T=T-Tn is the temperature deviation from standard temperature, G is the light intensity and KVis
the ratio coefficient of open circuit voltage to temperature. Open circuit voltage, short circuit current and
voltage-current corresponding to the maximum power are three important points of I-V characteristic of solar
panel. These points are changed by variations of atmospheric conditions. By using (4) and (5) which are derived
from PV model equations, short circuit current and open circuit voltage can be calculated in different
atmospheric conditions.
n
SCSC
G
G
TKII )( 1 (4)
TKVV VOCOC (5)
2.3. PMSG wind turbine
Wind turbine is one of part of an integrated system, which can be divided into two types i.e. fixed and
variable speed wind turbines. In fixed speed type, rotating speed of turbine is fixed and hence, frequency of
generated voltage remains constant, so it can be directly connected to the network. In this case, maximum
power can not always be extracted by wind. On the other hand, on variable speed, turbine can rotate at different
speeds, so maximum power can be generated in each wind speed by MPPT method [10]. The advantage of
0 5 10 15 20 25 30 35 40 45 50
0
50
100
150
200
Voltage (Volt)
Current(Ampere)
1200 W/m2
1000 W/m2
800 W/m2
600 W/m2
400 W/m2
200 W/m2](https://image.slidesharecdn.com/0313075-200724025151/85/Matlab-simulink-simulation-of-unified-power-quality-conditioner-battery-energy-storage-system-supplied-by-PV-wind-hybrid-using-fuzzy-logic-controller-4-320.jpg)
![Int J Elec & Comp Eng ISSN: 2088-8708
Matlab/simulink simulation of unified power quality conditioner-battery energy storage… (Amirullah)
1483
using a PMSG over a synchronous generator and DFIG machine because it has high efficiency and reliability.
Due to the elimination of rotor external excitation, machine size, and cost also decreases, making PMSG to be
more controlled easily with feedback control system. PMSG has become an attractive solution on wind
generation systems with variable speed wind turbine applications [15]. Figure 3 and 4 shows model of PMSG
wind turbine and wind turbin power characteristic curve.
Figure 3. Model of PMSG wind turbine
Figure 4. Wind turbin power characteristic curve
The output power of wind turbine can be expressed using (6), (7), and (8) [11].
wind
r
V
R
(6)
32
2
1
windpM VCRP
(7)
3
3
5
2
1
M
p
r
M
M CR
P
T (8)
Where λ is the speed-tip ratio, Vwindis wind speed, R is blade radius, ωr is rotor speed (rad/sec), ρ is the air
density, CP is the power coefficient, PM is the output mechanical power, and TM is output torque of wind turbine.
The CP coefficient is dependent on the pitch angle value, at which rotor blade can rotate along axis and tip-
speed ratio λ expressed in (9).
0 0.2 0.4 0.6 0.8 1 1.2 1.4
0.2
0.4
0.6
0.8
1
1.2
1.2 pu
Max. power at base wind speed (5 m/s) and beta = 0 deg
3 m/s
3.4 m/s
3.8 m/s
4.2 m/s
4.6 m/s
5 m/s
5.4 m/s
5.8 m/s
Turbine speed (pu of nominal generator speed)
Turbineoutputpower(pu)
Turbine Power Characteristics (Pitch angle beta = 2 deg)](https://image.slidesharecdn.com/0313075-200724025151/85/Matlab-simulink-simulation-of-unified-power-quality-conditioner-battery-energy-storage-system-supplied-by-PV-wind-hybrid-using-fuzzy-logic-controller-5-320.jpg)
![ ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 9, No. 3, June 2019 : 1479 - 1495
1484
)3.013(
)2(
sin)167.044.0(
PC (9)
Where β is a pitch angle blade. In a fixed pitch type, the value of β is set to a fixed value.
2.4. Control of series active filter
The main function of series active filter is as a sensitive load protection against a number of voltage
interference at PCC bus. The control strategy algorithm of the source and load voltage in series active filter
circuit is shown in Figure 5. It extracts the unit vector templates from the distorted input supply. Furthermore,
the templates are expected to be ideal sinusoidal signal with unity amplitude. The distorted supply voltages are
measured and divided by peak amplitude of fundamental input voltage Vm give in (10) [6].
222
3
2
scsbsam VVVV (10)
A three phase locked loop (PLL) is used in order to generate a sinusoidal unit vector templates with a
phase lagging by the use of sinus function. The reference load voltage signal is determined by multiplying the
unit vector templates with the peak amplitude of the fundamental input voltage Vm. The load reference voltage
(VLa
*
, VLb
*
, Vc
*
) is then compared against to sensed load voltage (VLa, VLb, VLc) by a pulse width modulation
(PWM) controller used to generate the desired trigger signal on series active filter.
Vsa
Vsb
Vcb
K
K
K
Three
Phase
PLL
Function to
get Sin (wt)
terms only
Sin (wt)
Sin (wt - 2π)
Sin (wt + 2π)
V*La
V*Lb
V*Lc
VLa VLb VLc
Gating
Signals
Sensed Load Voltage
Sensed
Source
Voltage
PWM
Voltage
Controller
Vm = Peak fundamental input
voltage magnitude
Figure 5. Control strategy of series active filter
2.5. Control of shunt active filter
The main function of shunt active filter is mitigation of power quality problems on the load side. The
control methodology in shunt active filter is that the absorbed current from the PCC bus is a balanced positive
sequence current including unbalanced sag voltage conditions in the PCC bus or unbalanced conditions or non-
linear loads. In order to obtain satisfactory compensation caesed by disturbance due to non-linear load, many
algorithms have been used in the literature. This research used instantaneous reactive power theory method "p-
q theory". The voltages and currents in Cartesian abc coordinates can be transformed to Cartesian αβ
coordinates as expressed in (11) and (12) [16].
c
b
a
V
V
V
v
v
2/32/31
2/12/11
(11)
Lc
Lb
La
i
i
i
i
i
2/32/31
2/12/11
(12)
The computation of the real power (p) and imaginary power (q) is showed in (13). The real power and
imaginary are measured instantaneously power and in matrix it is form is given as.The presence of oscillating
and average components in instantaneous power is presented in (14) [17].](https://image.slidesharecdn.com/0313075-200724025151/85/Matlab-simulink-simulation-of-unified-power-quality-conditioner-battery-energy-storage-system-supplied-by-PV-wind-hybrid-using-fuzzy-logic-controller-6-320.jpg)
![Int J Elec & Comp Eng ISSN: 2088-8708
Matlab/simulink simulation of unified power quality conditioner-battery energy storage… (Amirullah)
1485
i
i
vv
vv
q
p
(13)
;~ppp qqq ~ (14)
Where p = direct component of real power, p~ = fluctuating component of real power, q = direct
component of imaginary power, q~ = fluctuating component of imaginary power. The total imaginary power
(q) and the fluctuating component of real power are selected as power references and current references and
are utilized through the use of (15) for compensating harmonic and reactive power [18].
q
pp
vv
vv
vvi
i loss
c
c
~
1
22*
*
(15)
The signal lossp , is obtained from voltage regulator and is utilized as average real power. It can also
be specified as the instantaneous active power which corresponds to the resistive loss and switching loss of the
UPQC. The error obtained on comparing the actual DC-link capacitor voltage with the reference value is
processed in FLC, engaged by voltage control loop as it minimizes the steady state error of the voltage across
the DC link to zero. The compensating currents ( *
ci , *
ci ) as required to meet the power demand of load are
shown in (15). These currents are represented in α-β coordinates. The phase current is required to acquire
using (16) for compensation. These source phase currents ( *
sai , *
sbi , *
sci ) are represented in a-b-c axis obtained
from the compensating current in the α-β coordinates presented in (16) [18]. Figure 6 shows a control of shunt
active filter.
*
*
*
*
*
2/32/1
2/32/1
01
3
2
c
c
sc
sb
sa
i
i
i
i
i
(16)
Eq. 11
Eq. 12
Vsa
Vsb
Vsc
iα
iβ
vα
vβ
iLa
iLb
iLc
q
p
LPF
-1
vαβ
Icα
*
Hysterisis
Current
Controller
isa
*
isb
*
isc
*
Gating
Signals
isa isb isc
Sensed Source Current
Fuzzy Logic
Controller
VDC
VDC
*
Sensed
Source
Voltage
Sensed
Load
Current
-
+
lossp
vαβ
-q
Eq. 13 Eq. 15
Icβ
*
Eq. 16
-p
Figure 6. Control strategy of shunt active filter
The proposed model of UPQC-BES system supplied by three RE sources is shown in Figure 1.
From the figure, we can see that PV is connected to the DC link through a DC-DC boost converter circuit.
The PV partially distributes power to the load and the remains is transfered to three phase grid. The load
consists of non linear and unbalanced load. The non-linear load is a diode rectifier circuit with the RL load
type, while the unbalanced load is a three phase RC load with different R value on each phase. In order to
economically efficient, PV must always work in MPP condition. In this research, MPPT method used is P and
O algorithm. The model is also applied for UPQC-BES system which is supplied by wind and PV-wind hybrid
respectively. In order to operate properly, UPQC-BES system device must have a minimum DC link voltage
(Vdc). The value of common DC link voltage depends on the nstantaneous energy avialable to UPQC is defined
by in (17) [16]:](https://image.slidesharecdn.com/0313075-200724025151/85/Matlab-simulink-simulation-of-unified-power-quality-conditioner-battery-energy-storage-system-supplied-by-PV-wind-hybrid-using-fuzzy-logic-controller-7-320.jpg)
![ ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 9, No. 3, June 2019 : 1479 - 1495
1486
m
V
V LL
dc
3
22
(17)
where m is modulation index and VLL is the AC grid line voltage of UPQC. Considering that modulation index
as 1 and for line to line grid voltage (VLL=380 volt), the Vdc is obtained 620, 54 volt and selected as
650 volt.
The input of shunt active filter showed in Figure 5 is DC voltage (Vdc) and reference DC voltage
(Vdc
*
), while the output is lossp by using PI controller. Then, the lossp is as one of input variable to generate the
reference source current (Isa
*
, Isb
*
, and Isc
*
). The reference source current output is then compared to source
current (Isa, Isb, and Isc) by the current hysteresis control to generate trigger signal in IGBT circuit of shunt
active filter. In this research, FLC as DC voltage control algorithm on shunt active filter is proposed and
compared with PI controller. The FLC is capable of reduce oscilation and generate quick convergence
calculation during disturbances. This method is also used to overcome the weakness of PI controller in
determining proportional gain (Kp) and integral gain constant (Ki) which still use trial and error method.
2.6. Fuzzy logic controller
The research is started by determine lossp as the input variable to result the reference source current on
current hysteresis controller to generate trigger signal on the IGBT shunt active filter of UPQC using PI
controller (Kp=0.2 and Ki=1.5). By using the same procedure, lossp is also determined by using FLC.
The FLC has been widely used in recent industrial process because it has heuristic, simpler, more effective and
has multi rule based variables in both linear and non-linear system variations. The main components of FLC
are fuzzification, decision making (rulebase, database, reason mechanism) and defuzzification in
Figure 7. The output membership function is generated using inference blocks and the basic rules of FLC as
shown in Table 1.
Database
Reason
Mechanism
Rulebase
Fuzzy Logic Controller
Fuzzification DefuzzificationVdc
Vdc
*
Vdc-error
errordcV _
lossp
Figure 7. Diagram block of FLC
Table 1. Fuzzy Rule Base
Vdc-
error
NB NM NS Z PS PM PB
∆Vdc-
error
PB Z PS PS PM PM PB PB
PM NS Z PS PS PM PM PB
PS NS NS Z PS PS PM PM
Z NM NS NS Z PS PS PM
NS NM NM NS NS Z PS PS
NM NB NM NM NS NS Z PS
NB NB NB NM NM NS NS Z
The fuzzy rule algorithm collects a number of fuzzy control rules in a particular order. This rule is used
to control the system to meet the desired performance requirements and they are designed from a number of
intelligent system control knowledge. The fuzzy inference of FLC using Mamdani method related to max-min
composition. The fuzzy inference system in FLC consists of three parts: rule base, database, and reasoning
mechanism [19]. The FLC method is performed by determining input variables Vdc (Vdc-error) and delta Vdc
(ΔVdc-error), seven linguistic fuzzy sets, operation fuzzy block system (fuzzyfication, fuzzy rule base and
defuzzification), Vdc-error and ΔVdc-error during fuzzification process, fuzzy rule base table, crisp value to
determine lossp in defuzzification phase. The lossp is one of input variable to obtain compensating currents
( *
ci , *
ci ) in (16). During fuzzification process, a number of input variables are calculated and converted into
linguistic variables based on a subset called membership function. The error Vdc (Vdc-error) and delta error Vdc
(ΔVdc-error) are proposed input variable system and output variable is lossp . To translate these variables, each
input and output variable is designed using seven membership functions: Negative Big (NB), Negative Medium
(NM), Negative Small (NS), Zero (Z), Positive Small (PS), Positive Medium (PM) and Positive Big (PB). The
membership functions of crisp input and output are presented with triangular and trapezoidal membership
functions. The value of Vdc-error range from -650 to 650, ΔVdc-error from -650 to 650, and lossp from -100 to 100.
The input and output MFs are shown in Figure 8.](https://image.slidesharecdn.com/0313075-200724025151/85/Matlab-simulink-simulation-of-unified-power-quality-conditioner-battery-energy-storage-system-supplied-by-PV-wind-hybrid-using-fuzzy-logic-controller-8-320.jpg)
![ ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 9, No. 3, June 2019 : 1479 - 1495
1494
REFERENCES
[1] B. Han, et.al, “Combined Operation of Unified Power Quality Conditioner With Distributed Generation,” IEEE
Transactions on Power Delivery, Vol. 21, No. 1, pp. 330-338, Januari 2006.
[2] Vinod Khadkikar, “Enhanching Electric Power Quality UPQC: A. Comprehensive Overview,” IEEE Transactions
on Power Electronics, Vol. 27, No. 5, pp. 2284-2297, May 2012.
[3] Shafiuzzaman K. Khadem, et.al.,“Integration of UPQC for Power Quality Improvement in Distribution Genertation
Network,” ISGT Europe 2011, Manchester, United Kingdom, December 2011.
[4] NorshafinashSaudinet.al,“Study on The Effect of Distributed Generation towards Unified Power Quality Conditioner
Performance in Mitigating Voltage Sags,”IEEE International Conference on Power and Energy (PECon), Kota
Kinibalu, Sabah, Malaysia, pp. 695-7002-5, December 2012.
[5] S. N. Gohil, et.al.,“Three Phase Unified Power Quality Conditioner (UPQC) for Power Quality Improvement by
using UVTG technique,” 2013 International Conference on Renewable Energy and Sustainable Energy (ICRESE),
Coimbatore,pp 151 – 156,5-6 Dec. 2013.
[6] Yahia Bouzelata, et.al.,“Design and Simulation of Unified Power Quality Conditioner Fed by Solar Energy,”
International Journal of Hydrogen Energy Vol. 40, Elsevier Ltd, pp. 15267-15277, 2015.
[7] K.Ramalingeswara Rao et.al.,“Improvement of Power Quality using Fuzzy Logic Controller In Grid Connected
Photovoltaic Cell Using UPQC,”International Journal of Power Electronics and Drive System (IJPEDS) Vol. 5, No.
1, pp. 101~111 ISSN: 2088-8694, July 2014.
[8] Amirullah, et.al.,“Multi Units of Single Phase Distributed Generation Combined With Battery Energy Storage for
Phase Balancing in Distribution Network,” Vol. 78: 10–4, eISSN 2180–3722, Universiti Teknologi Malaysia (UTM)
Publisher, pp. 27–33, 2016.
[9] K.S. Srikanth, et.al.,,“Improvement of Power Quality for Microgrid Using Fuzzy Based UPQC
Controller,”International Conference on Electrical, Electronics, Signals, Communication and Optimization
(EESCO), Visakhapatnam, 24-25, pp. 1-6, Jan 2015.
[10] H. Toodeji, et.al., ,“Power Management and Performance Improvement in Integrated System of Variable Speed Wind
Turbin and UPQC,”International Conference onClean Electrical Power, 9-11, pp. 609-614, June 2009.
[11] M. Hosseipour, et.al., “Design and Simulation of UPQC to Improve Power Quality and Transfer Wind Energy to
Grid”, Journal of Applied Sciences 8 (21), pp. 3770-3782,2008.
[12] R. Bhavani, et.al, ,“Fuzzy Connected UPQC for Power Quality Enhancement in DFIG based Grid Connected Wind
Power System”, International Conference on Circuit, Power and Computing Technologies (ICCPCT), 19-20,
Nagercoil, pp. 1-7, March 2015.
[13] S. Rajesh Rajan, “Power Quality Improvement in Grid Connected Wind Energi System Using UPQC”, International
Journal of Research in Engineering and Technology (IJRET), Vol. 1, Issue 1, pp. 13-19, June 2013.
[14] Ali Reza Reisi, et.al, “Combined Photovoltaic and Unified Power Quality Controller to Improve Power Quality”,
Solar Energy 88, pp.154-162, 2013.
[15] Mallikarjuna and K. Amaresh, “Implementation of UPQC in Micro Grid Applications with Renewable Energy
Resources”, International Journal of Scientific Engineering and Technology Research, Vol. 4, Issue.36, pp. 7734-
7739, September-2015.
[16] Yash Pal, et.al.,“A Comparative Analysis of Different Magnetic Support Three Phase Four Wire Unified Power
Quality Conditioners – A Simulation Study”, Electrical Power and Energy System 47,pp. 437-447, 2013.
[17] Swapnil Y. Kamble and Madhukar M. Waware, “Unified Power Quality Conditioner for Power Quality
Improvement”, International Multi Conference on Automation, Computing, Communication, Control and
Compressed Sensing (iMac4s), Kottayam, India, 22-23 March 2013, pp. 432-437, 2013.
[18] MihirHembram and Ayan Kumar Tudu, “Mitigation of Power Quality Problems Using Unified PowerQuality
Conditioner (UPQC)”, Proceedings of the 2015 Third International Conference on Computer, Communication,
Control and Information Technology (C3IT), (2015), Hooghly, India, 7-8, pp.1-5, February 2015.
[19] Amirullah and AgusKiswantono, “Power Quality Enhancement of Integration Photovoltaic Generator to Grid under
Variable Solar Irradiance Level using MPPT-Fuzzy,” International Journal of Electrical and Computer Engineering
(IJECE),IAES Publisher,Vol. 6, No. 6, pp.2629-2642, December 2016.
APPENDIX
Three phase grid: RMS voltage 380 volt (L-L), 50 Hz, line impedance: Rs=0.1 Ohm Ls=15 mH; series
and shunt active filter: series inductance Lse=0.015 mH; shunt inductance Lsh=15 mH; injection transformers:
rating 10 kVA, 50 Hz, turn ratio (N1/N2)=1:1; non linear load: resistance RL=60 ohm, inductance LL=0.15 mH,
load impedance Rc=0.4 ohm and Lc=15 mH; unbalance load: resistance R1=24 ohm, R2=12 ohm, and R3=6
ohm, capacitance C1,C2, C3=2200 μF; DC-link: voltage VDC=650 volt and capacitance CDC=3000 μF; battery
energe storage: type=nickel metal hybrid, DC voltage=650 volt, rated capacity=200 Ah, initial SOC=100%,
inductance L1=6 mH, capacitance C1=200 μF; photovoltaic: active power=0.6 kW temperature=250
C,
irradiance=1000 W/m2
; PMSG wind turbine active power=0.6 kW, voltage=380 volt, 50 Hz, wind
speed=5 m/s, picth angle=2; PI controller: Kp=0.2, Ki=1.5; fuzzy model: method=mamdani,
composition=max-min; input membership function: error (Vdc)=trapmf, trimf delta error (∆Vdc)= trapmf, trimf;
output membership function: lossp =trapmf,trimf.](https://image.slidesharecdn.com/0313075-200724025151/85/Matlab-simulink-simulation-of-unified-power-quality-conditioner-battery-energy-storage-system-supplied-by-PV-wind-hybrid-using-fuzzy-logic-controller-16-320.jpg)
Matlab/simulink simulation of unified power quality conditioner-battery energy storage system supplied by PV-wind hybrid using fuzzy logic controller
- 1. International Journal of Electrical and Computer Engineering (IJECE) Vol. 9, No. 3, June 2019, pp. 1479~1495 ISSN: 2088-8708, DOI: 10.11591/ijece.v9i3.pp1479-1495 1479 Journal homepage: http://iaescore.com/journals/index.php/IJECE Matlab/simulink simulation of unified power quality conditioner-battery energy storage system supplied by PV-wind hybrid using fuzzy logic controller Amirullah1 , Ontoseno Penangsang2 , Adi Soeprijanto3 1,2,3 Department of Electrical Engineering, Faculty of Electrical Technology, ITS Surabaya, Indonesia 1 Study Program of Electrical Engineering, Faculty of Engineering, University of Bhayangkara Surabaya, Indonesia Article Info ABSTRACT Article history: Received May 2, 2018 Revised Nov 11, 2018 Accepted Dec 2, 2018 This paper presents performance analysis of Unified Power Quality Conditioner-Battery Energy Storage (UPQC-BES) system supplied by Photovoltaic (PV)-Wind Hybrid connected to three phase three wire (3P3W) of 380 volt (L-L) and 50 hertz distribution system. The performance of supply system is compared with two renewable energy (RE) sources i.e. PV and Wind, respectively. Fuzzy Logic Controller (FLC) is implemented to maintain DC voltage across the capacitor under disturbance scenarios of source and load as well as to compare the results with Proportional Intergral (PI) controller. There are six scenarios of disturbance i.e. (1) non-linear load (NL), (2) unbalance and nonlinear load (Unba-NL), (3) distortion supply and non-linear load (Dis-NL), (4) sag and non-linear load (Sag-NL), (5) swell and non-linear load (Swell-NL), and (6) interruption and non-linear load (Inter- NL). In disturbance scenario 1 to 5, implementation of FLC on UPQC-BES system supplied by three RE sources is able to obtain average THD of load voltage/source current slightly better than PI. Furthermore under scenario 6, FLC applied on UPQC-BES system supplied by three RE sources gives significantly better result of average THD of load voltage/source current than PI. This research is simulated using Matlab/Simulink. Keywords: Battery energy storage Photovoltaic Power quality PV-Wind Hybrid Total harmonic distortion UPQC Wind turbine Copyright © 2019Institute of Advanced Engineering and Science. All rights reserved. Corresponding Author: Amirullah, Department of Electrical Engineering, Institut Teknologi Sepuluh Nopember (ITS), Kampus ITS Keputih, Sukolilo, Surabaya 60111, Indonesia. Email: amirullah14@mhs.ee.its.ac.id, amirullah@ubhara.ac.id, am9520012003@yahoo.com 1. INTRODUCTION PV and wind are the most RE distributed generations (DGs) because they are able to convert sunlight and wind into power. PV and solar are the potential DGs sources since it only need sunlight to generate electricity, where the resources are available in abundance, free and relatively clean. Indonesia has enormous energy potential from the sun because it lies on the equator. Almost all areas of Indonesia get sunlight about 10 to 12 hours per day, with an average intensity of irradiation of 4.5 kWh/m2 or equivalent to 112.000 GW. The potential of wind energy in Indonesia generally has a speed between 4 to 5 m/s and classified as medium scale with potential capacity of 10 to 100 kW. The weakness of PV and wind turbine besides able to generate power, they also produces a number voltage and current harmonics resulted by presence of several types of PV and wind turbine devices and power converters as well as to increase a number of non-linear loads connected to the grid,so finally resulting in the decrease in power quality. In order to overcome and improve power quality due to presence of non-linear loads and integration of PV and wind turbine to grid, UPQC is a proposed. UPQC serves to compensate for source voltage quality problemsi.e. sag, swell unbalance, flicker, harmonics, and load current quality problems i.e. harmonics,
- 2. ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 3, June 2019 : 1479 - 1495 1480 unbalance, reactive currents, and neutral current. B. Han et. al and Vinod Khadkikar have investigated UPQC as one part of active power filter consisting of shunt and series active filters connected in parallel and serves as superior controller to overcome a number of power quality problems simultaneously. Series component has been responsible for reducing a number of interference on source side i.e. voltage sag/swell, flicker, unbalanced voltage, and harmonics. Shunt component has been responsible for addressing a current quality problems i.e. low power factor, load current harmonics, and unbalanced load. [1, 2]. UPQC based on RE has been investigated by some researchers. There are two methods used to overcome this problem i.e. using conventional and artificial intelligence. Shafiuzzaman K.K., et.al have proposed system includes a series inverter, shunt inverter, and a DG connected to a DC link through a rectifier using PI. The system was capable of increase source voltage quality i.e. sag and interruption and load current quality, as well as transfer of active power on/off grid mode [3]. The influence of DG on UPQC performance in reducing sag under conditions of some phase to ground faults to using DSTATCOM has been implemented by Norshafinash S., et.al. The DG was effective enough to help UPQC to improve sag [4]. It was connected in series with load resulting better sag mitigation compared to system without DG. Implementation of UPQC using UVTG method with PI to reduce sag, swell, voltage/current harmonics has been doneby S. N. Gohil, et.al. Simulation of voltage distortion was made by adding 5th and 7th harmonics at fundamental source voltage, resulting in a reduction of THD source current and THD load voltage [5]. UPQC supplied by PV panels using boost converter, PI, MPPT P and O, and p-q theory has been proposed by Yahia Bouzelata at.al [6]. The system was capable of compensate reactive power and reduce source current/load voltage harmonics, but did not discuss migitation of sag and interuption caused by PV penetration. Power quality enhancement of sag and source voltage harmonics on grid using UPQC supplied by PV array connected to DC link using PI compared with FLC has been done by Ramalengswara Rao, et.al. Combination of UPQC and PV using FLC can improve source voltage THD better than PI [7]. Amirullah et.al have researched a method for balancing current and line voltage, as a result of DGs of a single phase PV generator unit in randomly installed at homes through on a three phase four wire 220 kV and 50 Hz distribution line using BES and three of single phase bidirectional inverter. Both devices was capable of reduce unbalanced line current/voltage, but both of them were also capable of increase current/voltage harmonics on PCC bus [8]. Power quality migitation of UPQC on microgrid supplied by PV and wind turbine has been implemented by K S Srikanth et. al. It resulted that PI and FLC was able to improve power quality and reduce distortion in output power [9]. The UPQC-wind turbine to provide active power to overcome low sag and interruption voltage to grid has been investigated by H.Toodeji, et.al. The model was used VSC as a rectifier on generator output and controlled so that maximum power desired can be generated by different speed wind turbines using PI [10]. The UPQC-wind turbine connected to UPQC DC link was implemented by M. Hosseinpour, et al. The proposed combination using PI was capable of compensate swell, interruption voltage, and reactive power both on on/off grid [11]. R.Bhavani, et, al have researched on UPQC controlled by FLC to improve power quality in a DFIG wind turbine connected grid. FLC can improve power quality i.e. sag voltage and load current harmonics better than PI [12]. Power quality enhancement on wind turbine and BES with PI connected grid on PCC bus using UPQC has been implemented by S.RajeshRajan, et al. BES was installed to maintain and stabilize active power supply under different wind speed [13]. This research will analyze UPQC-BES performance supplied by PV-wind hybrid connected to 3P3W of 380 volt (L-L) and 50 hertz distribution system.The performance of supply system is compared with two RE sources i.e. PV and Wind, respectively. BES serves to store excess energy produced by three RE sources and distribute it to load if necessary, to prevent interruption voltage, and to adjust charging and discharging of energy in battery. BES is also expected to store excess power produced by three RE combinations and use it as backup power. FLC is proposed and compared with PI to control variable of DC voltage and DC reference voltage input to generate reference current source in current hysteresis controller on shunt active filter. DC voltage controller in shunt active filter and series active filter is used to migitate power quality of load voltage and source current.Performance of two controllers are used to determine load voltage, source current, load voltage THD, and source current THD based on IEEE 519. This paper is presented as follow. Section 2 describes proposed method, model of UPQC-BES system supplied by three RE sources i.e. PV, wind, and PV-wind hybrid, simulation parameters, PV and PMSG wind turbine model, series and shunt active filter, as well as application of PI and FLC method for proposed model. Section 3 shows results and analysis about performance of THD analysis on the proposed model of three RE sources connected to DC link of UPQC-BES system using PI and FLC. In this section, six disturbance scenarios are presented and the results are verified with Matlab/Simulink. Finally, this paper in concluded in Section 4.
- 3. Int J Elec & Comp Eng ISSN: 2088-8708 Matlab/simulink simulation of unified power quality conditioner-battery energy storage… (Amirullah) 1481 2. RESEARCH METHOD 2.1. Proposed method Figure 1 shows proposed model in this research. The RE sources based DGs used i.e. PV, Wind Turbine, and PV-Wind Turbine Hybrid connected to 3P3W distribution system with 380 volt (L-L) and 50 Hz frequency, through UPQC-BES system. PV array produce power under fixed temperature and radiation as well as connect to UPQC-DC link through a DC/DC boost converter. The maximum power point tracking (MPPT) method with Pertub and Observer (P and O) algorithms helps PV produce maximum power and generate output voltage, as the input voltage for DC/DC boost converter. The converter serves to adjust duty-cycle value and output voltage of PV as its input voltage to produce an output voltage corresponding to UPQC DC link voltage. Wind turbine type used is a permanent magnetic synchronous generator (PMSG) with variable speed and fixed voltage which generating power and is connected to UPQC DC link circuit through AC/DC bridge rectifier. The rectifier helps to change AC PSMG stator output voltage to DC voltage through LC circuit that serves to filter and smooth it before connected to UPQC-DC link.BES connected to the UPQC-DC link circuit serves as energy storage and is expected to overcome interruption voltage and overall help UPQC performance to improve voltage and current quality on source and load bus. Simulation parameters proposed in this study is shown in Appendix Section. Power quality analysis is performed on PV, Wind, PV-Wind Hybrid respectively, connected to 3P3W system through UPQC DC-link (on-grid) using BES circuit. Single phase circuit breakers (CBs) are used to connect and disconnect PV, Wind, and Hybrid PV-Wind respectively with UPQC DC-link. There are six disturbance scenarios i.e. (1) NL, (2) Unba-NL, (3) Dis-NL, (4) Sag-NL, (5) Swell-NL, and (6) Inter-NL. In scenario 1, the model is connected a non-linear load with RL and LL of 60 Ohm and 0.15 mH respectively. In scenario 2, the model is connected to non-linear load and during 0.3 s since t=0.2 s to t=0.5 s connected to unbalance three phase load with R1, R2, R3 as 6 Ohm, 12 Ohm, 24 Ohm respectively, and value of C1, C2, C3 as 2200 μF. In scenario 3, the model is connected to non-linear load and source voltage generating 5th and 7th harmonic components with individual harmonic distortion values of 5% and 2% respectively. In scenario 4, the model is connected to non-linear load and source experiences a sag voltage disturbance of 50% for 0.3 s between t=0.2 s to t=0.5 s. In scenario 5, the model is connected to a non-linear load and source experiences a swell voltage disturbance of 50% for 0.3 s between t=0.2 s to t=0.5 s. In scenario 6, the model is connected to non-linear load and source experiences an interruption voltage interference of 100% for 0.3 s between t=0.2 s to t=0.5 s. FLC is used as a DC voltage control in a shunt active filter to improve the power quality of the load voltage and current source and compare it with PI controller. Each disturbance scenario uses a PI controller and FLC so that the total of 12 disturbances. The result analysis of research was carried out i.e. (1) voltage and current on source or poin common coupling (PCC) bus, (2) voltage and current on load bus, (3) harmonic voltage and harmonic current on source bus and (4) harmonic voltage and harmonic current on load bus. The final phase is to compare performance of UPQC-BES system on-grid supplied by PV, Wind, PV-Wind Hybrid respectively using two controllers to improve power quality of load voltage and source current under six disturbance conditions. Shunt Active Filter Cdc Vdc Battery Energy Storage Boost Converter ishb ishc isha iLa iLb iLc CB1 3 Phase Grid LsRs Lse - + + + - - vca vcb vcc Cr Rr LsRs LsRs Lse Lse Vsa isa isb isc Vsb Vsc Series Active Filter LcRc LcRc LcRc VLa VLc VLb C1 R1 C2 R2 C3 R3 LL RL Rectifier Unbalance Load Non Linear Load PV Array Lsh Lsh Lsh Cpv L CL L1 C1 ipv + - vpv PCC Bus Load Bus PMSG Wind Rectifier isw-a isw-b isw-c CB2 L CL UPQC-BES System Figure 1. Proposed model of UPQC-BES system supplied by PV, Wind, and PV-Wind Hybrid
- 4. ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 3, June 2019 : 1479 - 1495 1482 2.2. Photovoltaic model Figure 2 shows the equivalent circuit and V-I characteristic of a solar panel. A solar panel is composed of several PV cells that have series, parallel, or series-parallel external connections [14]. (a) (b) Figure 2. Equivalent circuit and V-I characteristic of solar panel The V-I characteristic of a solar panel is showed in (1): P S t S oPV R IRV aV IRV III 1exp (1) where IPV is the photovoltaic current, Io is saturated reverse current, ‘a’ is the ideal diode constant, Vt=NSKTq-1 is the thermal voltage, NS is the number of series cells, q is the electron charge, K is the Boltzmann constant, T is the temperature of p–n junction, RS and RP are series and parallel equivalent resistance of the solar panels. IPV has a linear relation with light intensity and also varies with temperature variations. Io is dependent on temperature variations. The values of Ipv and Io are calculated as following (2) and (3): I G G TKII n nPVPV )( 1, (2) 1/)exp( , , tVnOC InSC o aVTKV TKI I (3) In which IPV,n, ISC,n and VOC,n are photovoltaic current, short circuit current and open circuit voltage in standard conditions (Tn=25 C and Gn=1000 Wm-2 ) respectively. KI is the coefficient of short circuit current to temperature, ∆T=T-Tn is the temperature deviation from standard temperature, G is the light intensity and KVis the ratio coefficient of open circuit voltage to temperature. Open circuit voltage, short circuit current and voltage-current corresponding to the maximum power are three important points of I-V characteristic of solar panel. These points are changed by variations of atmospheric conditions. By using (4) and (5) which are derived from PV model equations, short circuit current and open circuit voltage can be calculated in different atmospheric conditions. n SCSC G G TKII )( 1 (4) TKVV VOCOC (5) 2.3. PMSG wind turbine Wind turbine is one of part of an integrated system, which can be divided into two types i.e. fixed and variable speed wind turbines. In fixed speed type, rotating speed of turbine is fixed and hence, frequency of generated voltage remains constant, so it can be directly connected to the network. In this case, maximum power can not always be extracted by wind. On the other hand, on variable speed, turbine can rotate at different speeds, so maximum power can be generated in each wind speed by MPPT method [10]. The advantage of 0 5 10 15 20 25 30 35 40 45 50 0 50 100 150 200 Voltage (Volt) Current(Ampere) 1200 W/m2 1000 W/m2 800 W/m2 600 W/m2 400 W/m2 200 W/m2
- 5. Int J Elec & Comp Eng ISSN: 2088-8708 Matlab/simulink simulation of unified power quality conditioner-battery energy storage… (Amirullah) 1483 using a PMSG over a synchronous generator and DFIG machine because it has high efficiency and reliability. Due to the elimination of rotor external excitation, machine size, and cost also decreases, making PMSG to be more controlled easily with feedback control system. PMSG has become an attractive solution on wind generation systems with variable speed wind turbine applications [15]. Figure 3 and 4 shows model of PMSG wind turbine and wind turbin power characteristic curve. Figure 3. Model of PMSG wind turbine Figure 4. Wind turbin power characteristic curve The output power of wind turbine can be expressed using (6), (7), and (8) [11]. wind r V R (6) 32 2 1 windpM VCRP (7) 3 3 5 2 1 M p r M M CR P T (8) Where λ is the speed-tip ratio, Vwindis wind speed, R is blade radius, ωr is rotor speed (rad/sec), ρ is the air density, CP is the power coefficient, PM is the output mechanical power, and TM is output torque of wind turbine. The CP coefficient is dependent on the pitch angle value, at which rotor blade can rotate along axis and tip- speed ratio λ expressed in (9). 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0.2 0.4 0.6 0.8 1 1.2 1.2 pu Max. power at base wind speed (5 m/s) and beta = 0 deg 3 m/s 3.4 m/s 3.8 m/s 4.2 m/s 4.6 m/s 5 m/s 5.4 m/s 5.8 m/s Turbine speed (pu of nominal generator speed) Turbineoutputpower(pu) Turbine Power Characteristics (Pitch angle beta = 2 deg)
- 6. ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 3, June 2019 : 1479 - 1495 1484 )3.013( )2( sin)167.044.0( PC (9) Where β is a pitch angle blade. In a fixed pitch type, the value of β is set to a fixed value. 2.4. Control of series active filter The main function of series active filter is as a sensitive load protection against a number of voltage interference at PCC bus. The control strategy algorithm of the source and load voltage in series active filter circuit is shown in Figure 5. It extracts the unit vector templates from the distorted input supply. Furthermore, the templates are expected to be ideal sinusoidal signal with unity amplitude. The distorted supply voltages are measured and divided by peak amplitude of fundamental input voltage Vm give in (10) [6]. 222 3 2 scsbsam VVVV (10) A three phase locked loop (PLL) is used in order to generate a sinusoidal unit vector templates with a phase lagging by the use of sinus function. The reference load voltage signal is determined by multiplying the unit vector templates with the peak amplitude of the fundamental input voltage Vm. The load reference voltage (VLa * , VLb * , Vc * ) is then compared against to sensed load voltage (VLa, VLb, VLc) by a pulse width modulation (PWM) controller used to generate the desired trigger signal on series active filter. Vsa Vsb Vcb K K K Three Phase PLL Function to get Sin (wt) terms only Sin (wt) Sin (wt - 2π) Sin (wt + 2π) V*La V*Lb V*Lc VLa VLb VLc Gating Signals Sensed Load Voltage Sensed Source Voltage PWM Voltage Controller Vm = Peak fundamental input voltage magnitude Figure 5. Control strategy of series active filter 2.5. Control of shunt active filter The main function of shunt active filter is mitigation of power quality problems on the load side. The control methodology in shunt active filter is that the absorbed current from the PCC bus is a balanced positive sequence current including unbalanced sag voltage conditions in the PCC bus or unbalanced conditions or non- linear loads. In order to obtain satisfactory compensation caesed by disturbance due to non-linear load, many algorithms have been used in the literature. This research used instantaneous reactive power theory method "p- q theory". The voltages and currents in Cartesian abc coordinates can be transformed to Cartesian αβ coordinates as expressed in (11) and (12) [16]. c b a V V V v v 2/32/31 2/12/11 (11) Lc Lb La i i i i i 2/32/31 2/12/11 (12) The computation of the real power (p) and imaginary power (q) is showed in (13). The real power and imaginary are measured instantaneously power and in matrix it is form is given as.The presence of oscillating and average components in instantaneous power is presented in (14) [17].
- 7. Int J Elec & Comp Eng ISSN: 2088-8708 Matlab/simulink simulation of unified power quality conditioner-battery energy storage… (Amirullah) 1485 i i vv vv q p (13) ;~ppp qqq ~ (14) Where p = direct component of real power, p~ = fluctuating component of real power, q = direct component of imaginary power, q~ = fluctuating component of imaginary power. The total imaginary power (q) and the fluctuating component of real power are selected as power references and current references and are utilized through the use of (15) for compensating harmonic and reactive power [18]. q pp vv vv vvi i loss c c ~ 1 22* * (15) The signal lossp , is obtained from voltage regulator and is utilized as average real power. It can also be specified as the instantaneous active power which corresponds to the resistive loss and switching loss of the UPQC. The error obtained on comparing the actual DC-link capacitor voltage with the reference value is processed in FLC, engaged by voltage control loop as it minimizes the steady state error of the voltage across the DC link to zero. The compensating currents ( * ci , * ci ) as required to meet the power demand of load are shown in (15). These currents are represented in α-β coordinates. The phase current is required to acquire using (16) for compensation. These source phase currents ( * sai , * sbi , * sci ) are represented in a-b-c axis obtained from the compensating current in the α-β coordinates presented in (16) [18]. Figure 6 shows a control of shunt active filter. * * * * * 2/32/1 2/32/1 01 3 2 c c sc sb sa i i i i i (16) Eq. 11 Eq. 12 Vsa Vsb Vsc iα iβ vα vβ iLa iLb iLc q p LPF -1 vαβ Icα * Hysterisis Current Controller isa * isb * isc * Gating Signals isa isb isc Sensed Source Current Fuzzy Logic Controller VDC VDC * Sensed Source Voltage Sensed Load Current - + lossp vαβ -q Eq. 13 Eq. 15 Icβ * Eq. 16 -p Figure 6. Control strategy of shunt active filter The proposed model of UPQC-BES system supplied by three RE sources is shown in Figure 1. From the figure, we can see that PV is connected to the DC link through a DC-DC boost converter circuit. The PV partially distributes power to the load and the remains is transfered to three phase grid. The load consists of non linear and unbalanced load. The non-linear load is a diode rectifier circuit with the RL load type, while the unbalanced load is a three phase RC load with different R value on each phase. In order to economically efficient, PV must always work in MPP condition. In this research, MPPT method used is P and O algorithm. The model is also applied for UPQC-BES system which is supplied by wind and PV-wind hybrid respectively. In order to operate properly, UPQC-BES system device must have a minimum DC link voltage (Vdc). The value of common DC link voltage depends on the nstantaneous energy avialable to UPQC is defined by in (17) [16]:
- 8. ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 3, June 2019 : 1479 - 1495 1486 m V V LL dc 3 22 (17) where m is modulation index and VLL is the AC grid line voltage of UPQC. Considering that modulation index as 1 and for line to line grid voltage (VLL=380 volt), the Vdc is obtained 620, 54 volt and selected as 650 volt. The input of shunt active filter showed in Figure 5 is DC voltage (Vdc) and reference DC voltage (Vdc * ), while the output is lossp by using PI controller. Then, the lossp is as one of input variable to generate the reference source current (Isa * , Isb * , and Isc * ). The reference source current output is then compared to source current (Isa, Isb, and Isc) by the current hysteresis control to generate trigger signal in IGBT circuit of shunt active filter. In this research, FLC as DC voltage control algorithm on shunt active filter is proposed and compared with PI controller. The FLC is capable of reduce oscilation and generate quick convergence calculation during disturbances. This method is also used to overcome the weakness of PI controller in determining proportional gain (Kp) and integral gain constant (Ki) which still use trial and error method. 2.6. Fuzzy logic controller The research is started by determine lossp as the input variable to result the reference source current on current hysteresis controller to generate trigger signal on the IGBT shunt active filter of UPQC using PI controller (Kp=0.2 and Ki=1.5). By using the same procedure, lossp is also determined by using FLC. The FLC has been widely used in recent industrial process because it has heuristic, simpler, more effective and has multi rule based variables in both linear and non-linear system variations. The main components of FLC are fuzzification, decision making (rulebase, database, reason mechanism) and defuzzification in Figure 7. The output membership function is generated using inference blocks and the basic rules of FLC as shown in Table 1. Database Reason Mechanism Rulebase Fuzzy Logic Controller Fuzzification DefuzzificationVdc Vdc * Vdc-error errordcV _ lossp Figure 7. Diagram block of FLC Table 1. Fuzzy Rule Base Vdc- error NB NM NS Z PS PM PB ∆Vdc- error PB Z PS PS PM PM PB PB PM NS Z PS PS PM PM PB PS NS NS Z PS PS PM PM Z NM NS NS Z PS PS PM NS NM NM NS NS Z PS PS NM NB NM NM NS NS Z PS NB NB NB NM NM NS NS Z The fuzzy rule algorithm collects a number of fuzzy control rules in a particular order. This rule is used to control the system to meet the desired performance requirements and they are designed from a number of intelligent system control knowledge. The fuzzy inference of FLC using Mamdani method related to max-min composition. The fuzzy inference system in FLC consists of three parts: rule base, database, and reasoning mechanism [19]. The FLC method is performed by determining input variables Vdc (Vdc-error) and delta Vdc (ΔVdc-error), seven linguistic fuzzy sets, operation fuzzy block system (fuzzyfication, fuzzy rule base and defuzzification), Vdc-error and ΔVdc-error during fuzzification process, fuzzy rule base table, crisp value to determine lossp in defuzzification phase. The lossp is one of input variable to obtain compensating currents ( * ci , * ci ) in (16). During fuzzification process, a number of input variables are calculated and converted into linguistic variables based on a subset called membership function. The error Vdc (Vdc-error) and delta error Vdc (ΔVdc-error) are proposed input variable system and output variable is lossp . To translate these variables, each input and output variable is designed using seven membership functions: Negative Big (NB), Negative Medium (NM), Negative Small (NS), Zero (Z), Positive Small (PS), Positive Medium (PM) and Positive Big (PB). The membership functions of crisp input and output are presented with triangular and trapezoidal membership functions. The value of Vdc-error range from -650 to 650, ΔVdc-error from -650 to 650, and lossp from -100 to 100. The input and output MFs are shown in Figure 8.
- 9. Int J Elec & Comp Eng ISSN: 2088-8708 Matlab/simulink simulation of unified power quality conditioner-battery energy storage… (Amirullah) 1487 (a) (b) (c) Figure 8. Membership functions (a) Vdc-error, (b) ΔVdc-error, and (c) lossp After the Vdc-error and ΔVdc-error are obtained, then two input membership functions are converted to linguistic variables and uses them as input functions for FLC. The output membership function is generated using inference blocks and the basic rules of FLC as shown in Table 1. Finally the defuzzification block operates to convert generated lossp output from linguistic to numerical variable again. Then it becomes input variable for current hysteresis controller to produce trigger signal on the IGBT circuit of UPQC shunt active filter to reduce source current and load voltage harmonics. While simultaneously, it also improve power quality of 3P3W system under six scenarios due to the integration of three RE sources to UPQC-BES system. 3. RESULTS AND ANALYSIS The analysis of proposed model is investigated through determination of six disturbance scenarios i.e. (1) NL, (2) Unba-NL, (3) Dis-NL, (4) Sag-NL, (5) Swell-NL, and (6) Inter-NL. Each scenario of UPQC uses PI controller and FLC so total number of disturbance are 12 scenarios. By using Matlab/Simulink, the model is then executed according to the desired scenario to obtain curve of source voltage (Vs), load voltage (VL), compensation voltage (Vc), source current (Vs), load current (IL), and DC voltage DC link (Vdc). Then, THD value of source voltage, source current, load voltage, and load current in each phase as well as average THD value (Avg THD) are obtained base on the curves. THD in each phase is determined in one cycle started at t = 0.35 s. The results of average of source voltage, source current, load voltage, and load current of 3P3W system using UPQC-BES system supplied by three RE sources i.e. PV, wind, and PV-wind hybrid are presented in Table 2, 3, and 4. Next THD in each phase and average THD are showed in Table 5, 6 and 7. Table 2 shows UPQC-BES system supplied PV connected 3P3W system using PI and FLC, disturbance scenarios 1 to 5 produce an average load voltage above 307 V. While in scenario 6, FLC produces a higher average load voltage of 304.1 V than if using a PIof 286.7 V. If reviewed from average source current with PI, the highest and lowest average source currents are generated by disturbance scenarios 2 and 4 of 28.15
- 10. ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 3, June 2019 : 1479 - 1495 1488 A and 7,246 A. Whereas if using FLC the highest and lowest average source currents are achieved in scenario 2 and 6 of 28.84 A and 3.804 A. Table 3 shows UPQC-BES supplied by wind connected to 3P3W with PI and FLC, scenarios 1 to 5 is able to obtain a stable load voltage above 308 V. The difference is that in scenario 6, PI generates a load voltage of 274.8 V, otherwise if using FLC, load voltage increases to 306.4 V. If reviewed from source current using PI, the highest and lowest average source currents are resulted by scenarios 2 and 4 of 28.28 A and 7.417 A. Rather, if using FLC, the highest and lowest average source currents are achieved in scenario 2 and 6 of 28.82 A and 3,420 A. Table 4 shows UPQC-BES supplied PV-wind hybrid connected 3P3W using PI and FLC, scenarios 1 to 5 produce an average load voltage above 307 V. While in scenario 6, FLC generates an average load voltage 305.9 V higher than if using PI control of 283.9 V. If reviewed from average source current with PI, the highest and lowest average source currents are resulted by scenarios 2 and 4 of 28.21 A and 6,773 A. Otherwise if using FLC the highest and lowest average source currents are achieved in scenario 2 and 6 of 28.82 A and 3.640 A respectively. Table 2. Voltage and Current of 3P3W System Using UPQC-BES System Supplied by PV Table 3. Voltage and Current of 3P3W System Using UPQC-BES System Supplied by Wind Table 4. Voltage and Current of 3P3W System Using UPQC-BES System Supplied PV-Wind Hybrid
- 11. Int J Elec & Comp Eng ISSN: 2088-8708 Matlab/simulink simulation of unified power quality conditioner-battery energy storage… (Amirullah) 1489 Table 5. Harmonics of 3P3W System Using UPQC-BES System Supplied by PV Table 6. Harmonics of 3P3W System Using UPQC-BES System Supplied by Wind Table 7. Harmonics of 3P3W System Using UPQC-BES System Supplied by PV-Wind Hybrid Table 5 shows that average THD of load voltage (VL) of UPQC-BES system supplied by PV in 3P3W system for scenarios 1 to 5 using PI control is within limits in IEEE 519. The highest and lowest average THD load voltages are achieved under scenario 6 and 2 as 25.25% and 2.34% respectively. PI controller is also able to mitigate average THD source voltage in scenario 6 from not accessible (NA) to 25.25% on the load side. The highest and lowest average THD of source current are achieved in scenario 6 and 2 as 230.82% and 2.41%. Table 5 also indicates that average THD of load voltage of UPQC system supplied by PV with BES using FLC in scenarios 1 to 5, has fulfilled limits in IEEE 519. The highest and lowest average THD of load voltage are achieved under scenario 6 and 2 of 10.27% and 2.38%. The use of FLC method is also able to reduce average THD on source voltage in scenario 6 from NA to 10.27% on load side. The highest and lowest average THD of source current are achieved in scenario 6 and 2 of 32.30% and 2.413%. Table 6 shows that average THD of load voltage of UPQC-BES system supplied by wind in 3P3W system for interference scenarios 1 to 5 using PI is within limits prescribed in IEEE 519. The highest and lowest
- 12. ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 3, June 2019 : 1479 - 1495 1490 average THD load voltages is achieved under scenario 6 and scenario 2 as 28.14% and 1.90% respectively. PI controller is also able to reduce average THD source voltage in scenario 6 from NA to 28.14% on the load side. The highest and lowest average THD of source current are achieved in scenario 6 and 2 as 186.95% and 2.24% respectively. Table 6 also shows that average THD of load voltage of UPQC-BES supplied by wind using FLC in scenarios 1 to 5, has fulfilled limits prescribed in IEEE 519. The highest and lowest average THD of load voltage are achieved under scenario 6 and 4of 8.33% and 1.90%. FLC method is also able to reduce average THD on source voltage in scenario 6 from NA to 7.33% on load side. The highest and lowest average THD source current are achieved in scenario 6 and 2 of 27.33% and 2.22%. Table 7 shows that average THD of load voltage (VL) of UPQC-BES system supplied by PV-wind hybrid in 3P3W system for scenarios 1 to 5 using PI is within limits in IEEE 519. The highest and lowest average THD load voltages are achieved under scenario 6 and 2 as 27.8% and 2.22% respectively. PI is also able to mitigate average THD source voltage in scenario 6 from NA to 27.8% on load side. The highest and lowest average THD of source current are achieved in scenario 6and 2 as 302.45% and 2.407% respectively. Table 7 also indicates that average THD of load voltage of UPQC-BES system supplied by PV-wind hybrid using FLC for scenarios 1 to 5, has fulfilled limits in IEEE 519. The highest and lowest average THD of load voltage are achieved under scenario 6 and 2 of 8.48% and 2.26%. FLC method is also able to reduce average THD on source voltage in scenario 6 from NA to 8.48% on load side. The highest and lowest average THD of source current are achieved in scenario 6and 2 of 28.53% and 2.446%. Overall for UPQC-BES systemsystem supplied by three RE sources in scenarios 1 to 5 using PI and FLC is able to improve average THD of source current better on average THD of load current. Figure 9 presents UPQC-BES system performance supplied by three RE sources using FLC in scenario 6. Figure 9 a(i) shows that in scenario 6, UPQC-BES system supplied by PV at t=0.2 s to t=0.5 s, average source voltage (VS) falls as 100% to 0.4062 V. During the disturbance, PV is able to generate power to UPQC DC link and injecting full average compensation voltage (VC) in Figure 9 a(iii) through injection transformer on series active filter so that average load voltage (VL) in Figure 9 a(ii) remains stable at 304.1 V. As long as fault period, although nominal of average source current (IS) in Figure 9 a(iv) drops to 3.804 A, combination of PV and BES is able to generate power, store excess energy of PV, and inject current into load through shunt active filter so that average load current (IL) in Figure 9 a(v) remains as 8.421 A. Figure 9 b(i) presents on UPQC-BES supplied by wind at t=0.2 s to t=0.5 s average source voltage (VS) drops 100% to 0.3828 V. During the disturbance, wind is able to generate power to UPQC DC link and injecting full average compensation voltage (VC) in Figure 9 b(iii) through injection transformer on series active filter so that average load voltage (VL) in Figure 9 b(ii) remains stable at 306.4 V. During fault period, although nominal of average source current (IS) falls to 3.420 A, combination of wind and BES is able to generate power, store excess energy of wind, and inject current into load through shunt active filter so that average load current (IL) in Figure 9 b(v) remains as 8.569 A. Figure 9 c(i) indicates on UPQC-BES supplied by PV-Wind Hybrid at t=0.2 s to t=0.5 s average source voltage (VS) drops 100% to 0.4175 V. During the disturbance, PV-wind hybrid is able to generate power to UPQC DC link and injecting full average compensation voltage (VC) in Figure 9 c(iii) through injection transformer on series active filter so that average load voltage (VL) in Figure 9 b(ii) remains stable at 305.9 V. As long as disturbance period, although nominal of average source current (IS) falls to 3.640 A, combination of PV-wind hybrid and BES is able to generate power, store excess energy of wind, and inject current into load through shunt active filter so that average load current (IL) in Figure 9 c(v) remains as 8.488 A. Figure 10 shows spectra of load voltage harmonics on phase A of UPQC-BES system supplied by three RE sources using FLC in scenario 6. Figure 11 shows performance of average THD of load voltage and source current on UPQC-BES system supplied by three RE sources. Figure 11(a) shows that in scenario 1 to 5, the implementation of FLC on UPQC-BES system supplied by three RE sources is able to obtain average THD of load voltage slightly better than PI controller and both method have already met the limit in IEEE 519. Further under scenario 6, FLC applied on UPQC-BES system supplied by three RE sources gives significantly better result of average THD of load voltage than PI controller. In six disturbance scenarios, both PI controller and FLC applied on UPQC-BES system supplied by three RE sources, PV is able to obtain the highest average THD of load voltages. Figure 11(b) shows that in disturbance scenario 1 to 5, implementation of FLC on UPQC-BES system supplied by three RE sources is able to obtain average THD of source current slightly better than PI controller. Furthermore under scenario 6, FLC applied on UPQC-BES system supplied by three RE sources gives significantly better result of average THD of source current than PI controller. Both PI controller and FLC on UPQC-BES system supplied by three RE sources in six disturbance scenarios, PV is able to obtain the highest average THD of source current.
- 13. Int J Elec & Comp Eng ISSN: 2088-8708 Matlab/simulink simulation of unified power quality conditioner-battery energy storage… (Amirullah) 1491 Figure 9. UPQC-BES system performance using FLC in scenario 6 (Inter-NL): (a) UPQC-BES + PV; (b) UPQC-BES + Wind; (c) UPQC-BES + PV-Wind Hybrid (continue)
- 14. ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 3, June 2019 : 1479 - 1495 1492 Figure 9. UPQC-BES system performance using FLC in scenario 6 (Inter-NL): (a) UPQC-BES + PV; (b) UPQC-BES + Wind; (c) UPQC-BES + PV-Wind Hybrid (a) UPQC-BES + PV (b) UPQC-BES + Wind (c) UPQC-BES + PV-Wind Hybrid Figure 10. Spectra of load voltage harmonics on phase A of UPQC-BES using FLC in scenario 6 (Inter-NL)
- 15. Int J Elec & Comp Eng ISSN: 2088-8708 Matlab/simulink simulation of unified power quality conditioner-battery energy storage… (Amirullah) 1493 (a) Average load voltage harmonics (b) Average source current harmonics Figure 11. Performance of UPQC-BES supplied by three RE sources using PI and FLC 4. CONCLUSION Comparative performance analysis of UPQC-BES system supplied by three RE sources i.e. PV, wind, and PV-wind hybrid respectively using PI controller and FLC have been discussed. In disturbance scenario 1 to 5, implementation of FLC on UPQC-BES system supplied by three RE sources is able to obtain average THD of load voltage slightly better than PI and both methods have already met the limit in IEEE 519. Under scenario 6, FLC applied on UPQC-BES system supplied by three RE sources gives significantly better result average THD of load voltage than PI. In six disturbance scenarios, both PI and FLC applied on UPQC- BES system supplied by three RE sources, PV is able to obtain the highest average THD of load voltage. In interference scenario 1 to 5, FLC method on UPQC-BES system supplied by three RE sources is able to obtain average THD of source current slightly better than PI. Furthermore under scenario 6, FLC applied on UPQC-BES system supplied by three RE sources gives significantly better result of average THD of source current than PI. Both PI and FLC on UPQC-BES system supplied by three RE sources in six scenarios, PV is able to obtain the highest average THD of source current. ACKNOWLEDGEMENTS The authors would like to acknowledge to Ministry of Research, Technology, and Higher Education, Republic of Indonesia, for financial support by BPP-DN Scholarships to pursue Doctoral Program in Department of Electrical Engineering, Faculty of Electrical Technology, ITS Surabaya. NL Unba-NL Dis-NL Sag-NL Swell-NL Inter-NL 0 50 100 150 200 250 300 Disturbance Scenarios AverageTHDSourceCurrent(%) PV-PI PV-FLC Wind-PI Wind-FLC PV+Wind Hybrid-PI PV+Wind Hybrid-FLC
- 16. ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 3, June 2019 : 1479 - 1495 1494 REFERENCES [1] B. Han, et.al, “Combined Operation of Unified Power Quality Conditioner With Distributed Generation,” IEEE Transactions on Power Delivery, Vol. 21, No. 1, pp. 330-338, Januari 2006. [2] Vinod Khadkikar, “Enhanching Electric Power Quality UPQC: A. Comprehensive Overview,” IEEE Transactions on Power Electronics, Vol. 27, No. 5, pp. 2284-2297, May 2012. [3] Shafiuzzaman K. Khadem, et.al.,“Integration of UPQC for Power Quality Improvement in Distribution Genertation Network,” ISGT Europe 2011, Manchester, United Kingdom, December 2011. [4] NorshafinashSaudinet.al,“Study on The Effect of Distributed Generation towards Unified Power Quality Conditioner Performance in Mitigating Voltage Sags,”IEEE International Conference on Power and Energy (PECon), Kota Kinibalu, Sabah, Malaysia, pp. 695-7002-5, December 2012. [5] S. N. Gohil, et.al.,“Three Phase Unified Power Quality Conditioner (UPQC) for Power Quality Improvement by using UVTG technique,” 2013 International Conference on Renewable Energy and Sustainable Energy (ICRESE), Coimbatore,pp 151 – 156,5-6 Dec. 2013. [6] Yahia Bouzelata, et.al.,“Design and Simulation of Unified Power Quality Conditioner Fed by Solar Energy,” International Journal of Hydrogen Energy Vol. 40, Elsevier Ltd, pp. 15267-15277, 2015. [7] K.Ramalingeswara Rao et.al.,“Improvement of Power Quality using Fuzzy Logic Controller In Grid Connected Photovoltaic Cell Using UPQC,”International Journal of Power Electronics and Drive System (IJPEDS) Vol. 5, No. 1, pp. 101~111 ISSN: 2088-8694, July 2014. [8] Amirullah, et.al.,“Multi Units of Single Phase Distributed Generation Combined With Battery Energy Storage for Phase Balancing in Distribution Network,” Vol. 78: 10–4, eISSN 2180–3722, Universiti Teknologi Malaysia (UTM) Publisher, pp. 27–33, 2016. [9] K.S. Srikanth, et.al.,,“Improvement of Power Quality for Microgrid Using Fuzzy Based UPQC Controller,”International Conference on Electrical, Electronics, Signals, Communication and Optimization (EESCO), Visakhapatnam, 24-25, pp. 1-6, Jan 2015. [10] H. Toodeji, et.al., ,“Power Management and Performance Improvement in Integrated System of Variable Speed Wind Turbin and UPQC,”International Conference onClean Electrical Power, 9-11, pp. 609-614, June 2009. [11] M. Hosseipour, et.al., “Design and Simulation of UPQC to Improve Power Quality and Transfer Wind Energy to Grid”, Journal of Applied Sciences 8 (21), pp. 3770-3782,2008. [12] R. Bhavani, et.al, ,“Fuzzy Connected UPQC for Power Quality Enhancement in DFIG based Grid Connected Wind Power System”, International Conference on Circuit, Power and Computing Technologies (ICCPCT), 19-20, Nagercoil, pp. 1-7, March 2015. [13] S. Rajesh Rajan, “Power Quality Improvement in Grid Connected Wind Energi System Using UPQC”, International Journal of Research in Engineering and Technology (IJRET), Vol. 1, Issue 1, pp. 13-19, June 2013. [14] Ali Reza Reisi, et.al, “Combined Photovoltaic and Unified Power Quality Controller to Improve Power Quality”, Solar Energy 88, pp.154-162, 2013. [15] Mallikarjuna and K. Amaresh, “Implementation of UPQC in Micro Grid Applications with Renewable Energy Resources”, International Journal of Scientific Engineering and Technology Research, Vol. 4, Issue.36, pp. 7734- 7739, September-2015. [16] Yash Pal, et.al.,“A Comparative Analysis of Different Magnetic Support Three Phase Four Wire Unified Power Quality Conditioners – A Simulation Study”, Electrical Power and Energy System 47,pp. 437-447, 2013. [17] Swapnil Y. Kamble and Madhukar M. Waware, “Unified Power Quality Conditioner for Power Quality Improvement”, International Multi Conference on Automation, Computing, Communication, Control and Compressed Sensing (iMac4s), Kottayam, India, 22-23 March 2013, pp. 432-437, 2013. [18] MihirHembram and Ayan Kumar Tudu, “Mitigation of Power Quality Problems Using Unified PowerQuality Conditioner (UPQC)”, Proceedings of the 2015 Third International Conference on Computer, Communication, Control and Information Technology (C3IT), (2015), Hooghly, India, 7-8, pp.1-5, February 2015. [19] Amirullah and AgusKiswantono, “Power Quality Enhancement of Integration Photovoltaic Generator to Grid under Variable Solar Irradiance Level using MPPT-Fuzzy,” International Journal of Electrical and Computer Engineering (IJECE),IAES Publisher,Vol. 6, No. 6, pp.2629-2642, December 2016. APPENDIX Three phase grid: RMS voltage 380 volt (L-L), 50 Hz, line impedance: Rs=0.1 Ohm Ls=15 mH; series and shunt active filter: series inductance Lse=0.015 mH; shunt inductance Lsh=15 mH; injection transformers: rating 10 kVA, 50 Hz, turn ratio (N1/N2)=1:1; non linear load: resistance RL=60 ohm, inductance LL=0.15 mH, load impedance Rc=0.4 ohm and Lc=15 mH; unbalance load: resistance R1=24 ohm, R2=12 ohm, and R3=6 ohm, capacitance C1,C2, C3=2200 μF; DC-link: voltage VDC=650 volt and capacitance CDC=3000 μF; battery energe storage: type=nickel metal hybrid, DC voltage=650 volt, rated capacity=200 Ah, initial SOC=100%, inductance L1=6 mH, capacitance C1=200 μF; photovoltaic: active power=0.6 kW temperature=250 C, irradiance=1000 W/m2 ; PMSG wind turbine active power=0.6 kW, voltage=380 volt, 50 Hz, wind speed=5 m/s, picth angle=2; PI controller: Kp=0.2, Ki=1.5; fuzzy model: method=mamdani, composition=max-min; input membership function: error (Vdc)=trapmf, trimf delta error (∆Vdc)= trapmf, trimf; output membership function: lossp =trapmf,trimf.
- 17. Int J Elec & Comp Eng ISSN: 2088-8708 Matlab/simulink simulation of unified power quality conditioner-battery energy storage… (Amirullah) 1495 BIOGRAPHIES OF AUTHORS Amirullah was born in Sampang East Java Indonesia, in 1977. He received bachelor and master degree in electrical engineering from University of Brawijaya Malang and ITS Surabaya, in 2000 and 2008, respectively. He also worked as a lecturer in University of Bhayangkara Surabaya. He is currently working toward the doctoral degree, in electrical engineering in Power System and Simulation Laboratory (PSSL) ITS Surabaya. His research interest includes power distribution modeling and simulation, power quality, harmonics migitation, design of filter/PFC, and RE. Ontoseno Penangsang was born in Madiun East Java Indonesia, in 1949. He received bachelor degree in electrical engineering from ITS Surabaya, in 1974. He received M.Sc and Ph.D degree in Power System Analysis from University of Wisconsin, Madison, USA, in 1979 and 1983, respectively. He is currently a professsor at Department of Electrical Engineering and the head of PSSL ITS Surabaya. He has a long experience and main interest in power system analysis (with renewable energy sources), design of power distribution, power quality, and harmonic mitigation in industry. Adi Soeprijanto was born in Lumajang East Java Indonesia, in 1964. He received bachelor in electrical engineering from ITB Bandung, in 1988. He received master of electrical engineering in control automatic from ITB Bandung. He continued his study to Doctoral Program in Power System Control in Hiroshima University Japan and was finished it’s in 2001. He is currently a professsor at Department of Electrical Engineering and member of PSSL in ITS Surabaya. His main interest includes power system analysis, power system stability control, and power system dynamic stability. He had already achieved a patent in optimum operation of power system.