Syncronous DQ PV Grid Control
Syncronous DQ PV Grid Control
Advisor:
Mengesha Mamo (Associate Professor)
October, 2019
By
Wondimagegn Minda
Graduate Committee
Name Signature
Place: Addis Ababa Institute of Technology, Addis Ababa University, Addis Ababa
This thesis has been submitted for examination with my approval as a university advisor.
i|Page
Abstract
Grid connected solar PV systems are increasing dramatically from year to year due to the
increase of global interest in renewable energy sources and growth in energy demand. In thesis
work, the designing and modeling of grid connected PV system using synchronously rotating
reference frame (SRRF) has been carried out. In this thesis, the important points and literature
regarding the elements in a PV system, an inverter and grid connection system have been studied
and reviewed.
The photovoltaic model is used to simulate actual PV array characteristics, and an MPPT method
using incremental conductance method is proposed in order to control the DC-DC boost
converter to extract the maximum power from the PV array. The inverter used in this thesis work
is three-phase two-level voltage source inverter (VSI) which is controlled using synchronously d-
q reference frame with space vector pulse width modulation (SVPWM). The control structure for
the inverter is designed by synchronous reference frame with phase locked loop (PLL). The
phase locked loop is used to acquire the available necessary information of grid voltage. An
LCL-filter is used to interconnect the output of inverter to the grid and used to attenuate the total
harmonics distortion (THD) present in the output current and voltage of the inverter. The system
model is developed using MATLAB/SIMULINK environment.
The efficiency of the overall system is around 97% which is improved more than the previous
related work whose efficiency is around 96%. The total harmonic distortion of the grid current is
0.88% which is less than the previous work that is 1.42%. For resistive load change, the THD of
the grid current and inverter current are 0.88%, 0.91% and 0.93%, and 0.86%, 0.87% and 0.87%
for a load of 1.5kW, 2.5kW and 3.5kW respectively. For a nonlinear load change, the THD of
the grid current and inverter current are 0.61%, 0.51% and 0.50%, and 0.12%, 0.12% and 0.12%
for a load of 100Ω, 125Ω and 320Ω and C = 3000µF respectively.
Keywords: DC-DC Converter, MPPT controller, Inverter control scheme, Grid synchronization,
Phase locked loop (PLL), Grid connected solar PV system.
ii | P a g e
Table of Contents
Contents Page
Acknowledgement ........................................................................................................................... i
Abstract ........................................................................................................................................... ii
Table of Contents ........................................................................................................................... iii
List of Tables ................................................................................................................................. vi
List of Figures ............................................................................................................................... vii
List of Abbreviations ...................................................................................................................... x
Chapter One:Introduction ........................................................................................................... 1
1.1. Background of the Thesis ..................................................................................................... 1
1.2. Statement of problem ........................................................................................................... 2
1.3. Objective of the thesis .......................................................................................................... 2
1.3.1. General Objective .......................................................................................................... 2
1.4. Scope of the thesis ................................................................................................................ 3
1.5. Methodology ........................................................................................................................ 3
1.6. Thesis Organization.............................................................................................................. 5
Chapter Two:Literature Review ................................................................................................. 6
2.1. Overview of Photovoltaic Cell ............................................................................................. 7
2.1.1. Types of PV Cells .......................................................................................................... 7
2.2. Types of Solar PV System ................................................................................................... 8
2.2.1. Grid connected PV system ............................................................................................ 8
2.2.2. Stand-Alone Photovoltaic Systems ............................................................................... 8
2.3. Maximum Power Point Tracking ......................................................................................... 9
2.3.1. Perturb and Observe Algorithm ..................................................................................... 9
2.3.2. Incremental Conductance Algorithm........................................................................... 10
2.4. DC-DC Boost converter ..................................................................................................... 12
2.5. DC-AC Inverter in PV System ........................................................................................... 12
2.5.1. Current Source Inverter ............................................................................................... 12
2.5.2. Voltage Source Inverter ............................................................................................... 13
2.6. Active and Reactive Powers and DC-Link Voltage Controller ......................................... 13
iii | P a g e
2.6.1. Clarke‟s and Park‟s Transformation ............................................................................ 13
2.7. Grid Synchronization ......................................................................................................... 14
2.7.1. Phase Locked Loop (PLL) ........................................................................................... 14
2.8. Power Decoupling .............................................................................................................. 15
2.9. Modulation Techniques for Inverter .................................................................................. 16
2.10. Filter as Grid Interface ..................................................................................................... 16
Chapter Three:System Modeling of Photovoltaic Inverter .................................................... 18
3.1. Photovoltaic Cell ................................................................................................................ 18
3.2. Modelling of PV Modules .................................................................................................. 18
3.3. Modeling of Boost converter.............................................................................................. 25
3.3.1. Design of DC-DC converter ........................................................................................ 26
3.3.2. Incremental Conductance MPPT ................................................................................. 28
3.4. DC-AC Inverter .................................................................................................................. 29
3.4.1. Design of LCL-Filter ................................................................................................... 30
3.5. Space Vector PWM Strategy Control ................................................................................ 33
3.5.1. Principle of Space Vector PWM ................................................................................. 33
Chapter Four:Grid Connected PV System............................................................................... 38
4.1. Control Structure of Grid Connected Inverter.................................................................... 38
4.2. Control Structure of the Inverter ........................................................................................ 40
4.2.1. Representation of Three Phase Variables in Stationary RF ........................................ 41
4.2.2. Representation of Three Phase Variables in SRRF ..................................................... 42
4.2.3. Grid Synchronization ................................................................................................... 43
4.2.4. Control Loops for Inverter Control ............................................................................. 44
4.3. Grid with Load ................................................................................................................... 47
4.4. Block Diagram of Overall System ..................................................................................... 48
Chapter Five:Simulation Results and Discussion .................................................................... 50
5.1. Introduction ........................................................................................................................ 50
5.2. PV Array ............................................................................................................................ 50
5.3. DC-DC Converter .............................................................................................................. 54
5.4. DC-Link Voltage ................................................................................................................ 55
5.5. DC-AC Inverter .................................................................................................................. 57
iv | P a g e
5.6. Three Phase Voltage and Current Waveforms ................................................................... 60
5.7. Active and Reactive Reference Currents ........................................................................... 62
5.8. Output Power of the System .............................................................................................. 63
5.9. Effect of Resistive Load Change ........................................................................................ 64
5.10. Effect of Non-Linear Load Change .................................................................................. 67
Chapter Six:Conclusion and Future Work .............................................................................. 71
6.1. Conclusion.......................................................................................................................... 71
6.2. Future Work ....................................................................................................................... 72
References ..................................................................................................................................... 73
Appendices .................................................................................................................................... 78
Appendix A: MATLAB CODES to obtain I-V and P-V Curve ............................................... 78
Appendix B: MATLAB code to design inverter filter .............................................................. 80
v|Page
List of Tables
Table 3.1. Design parameters for PV array[30] ............................................................................ 21
Table 3.2. Design parameters for boost converter ........................................................................ 27
Table 3.3. Boost converter calculated parameters ........................................................................ 27
Table 3.4. Inverter and grid parameters ........................................................................................ 33
Table 3.5. Sector Identification ..................................................................................................... 35
Table 3.6. Seven segment switching sequence [43] ..................................................................... 37
Table 4.1. Grid and controller parameters used for simulation .................................................... 49
Table 5.1. A summary for different values of resistive load ....................................................... 67
Table 5.2. A summary for different values of nonlinear load ...................................................... 70
vi | P a g e
List of Figures
Figure 1.1 General Block diagram of Grid Connected PV system ................................................. 4
Figure 2.1 Block diagram of Grid connected solar PV system....................................................... 8
Figure 2.2 Block diagram of Standalone PV system with Battery storage ..................................... 9
Figure 2.3 Flowchart of Perturb and Observe Algorithm [13] ..................................................... 10
Figure 2.4 Flowchart of Incremental Conductance Algorithm [13] ............................................. 11
Figure 2.5 Incremental Conductance Method for MPPT.............................................................. 12
Figure 2.6 Clarke‟s Transformation .............................................................................................. 13
Figure 2.7 Park‟s Transformation ................................................................................................. 14
Figure 2.8 PLL loop [18] .............................................................................................................. 15
Figure 2.9 Capacitor location in single stage ................................................................................ 16
Figure 2.10 Capacitor location in multistage system .................................................................... 16
Figure 2.11 L, LC and LCL filters [22] ........................................................................................ 17
Figure 3.1 PV Cell, Module and Array [23] ................................................................................ 18
Figure 3.2 The Equivalent Circuit of Photovoltaic cell [25] ........................................................ 19
Figure 3.3 Block diagram of temperature conversion .................................................................. 21
Figure 3.4 Block diagram of photo current ................................................................................... 22
Figure 3.5 Block diagram of diode saturation current .................................................................. 22
Figure 3.6 Block diagram of reverse saturation current ............................................................... 23
Figure 3.7 Block diagram of shunt current ................................................................................... 23
Figure 3.8 Block diagram of output current.................................................................................. 24
Figure 3.9 Block diagram of Solar PV System for a single PV module ...................................... 24
Figure 3.10 Block diagram of Solar PV System for a PV array ................................................... 25
Figure 3.11 Circuit diagram of Boost converter ........................................................................... 26
Figure 3.12 Block Diagram of Incremental Conductance ............................................................ 29
Figure 3. 13 Three phase two-level inverter ................................................................................. 30
Figure 3. 14 Circuit diagram of LCL-filter [38] ........................................................................... 30
Figure 3. 15 Bode Plot of LCL filter: (a) without Damping Resistor (b) with Damping Resistor 32
Figure 3. 16 Representation of space vector in complex plane [44] ............................................. 34
Figure 3.17 Principle of space vector time calculation ................................................................. 36
Figure 4.1 Control Structure in Natural reference Frame [18] ..................................................... 39
vii | P a g e
Figure 4.2 General Structure for Implementation of Control in Stationary RF [18] .................... 39
Figure 4.3 General Control Structure in Synchronous Reference Frame [18] ............................. 40
Figure 4.4 Current Control Loop; Cross coupling [47] ................................................................ 40
Figure 4.5 Internal Structure of Phase Locked Loop .................................................................... 44
Figure 4.6 DC-Link Voltage Controller........................................................................................ 45
Figure 4.7 Internal Current Controller .......................................................................................... 46
Figure 4.8 Overall Voltage and Current controller for Active and Reactive Power Control ....... 47
Figure 4.9 Model of Grid Connected Load ................................................................................... 47
Figure 4.10 The Complete Model of the Grid Connected PV System ......................................... 49
Figure 5.1 Effect of temperature changes on I-V and P-V Curve ................................................ 50
Figure 5.2 Effect of solar irradiance changes on I-V and P-V Curve ........................................... 51
Figure 5.3 I-V and P-V Curve at different temperature found from datasheet ............................. 51
Figure 5.4 I-V and P-V Curve at different irradiance found from datasheet ................................ 52
Figure 5.5 PV array characteristic curve for I-V and P-V under STC .......................................... 53
Figure 5.6 I-V and P-V Curve for PV Array Measured from MATLAB/Simulink ..................... 53
Figure 5.7 Temperature and Irradiance at different level ............................................................. 54
Figure 5.8 Output of PV array ...................................................................................................... 55
Figure 5.9 DC-link Voltage and Reference Voltage Signal ......................................................... 56
Figure 5.10 Error Signal Between DC-link Voltage and Reference Voltage ............................... 56
Figure 5.11 Inverter Output Voltage Vab ..................................................................................... 57
Figure 5.12 Output of Inverter before and after filter for Current waveform ............................... 58
Figure 5.13 Output of Inverter before and after filter for Voltage waveform .............................. 58
Figure 5.14 FFT Analysis of Inverter Output Current .................................................................. 59
Figure 5.15 FFT Analysis of Inverter Output Voltage ................................................................. 59
Figure 5.16 Output Waveform of Three Phase Voltage for Inverter and Grid ............................. 60
Figure 5.17 Output Waveform of Three Phase Current for Inverter and Grid ............................. 61
Figure 5.18 FFT Analysis of Grid Current ................................................................................... 61
Figure 5.19 FFT Analysis of Grid Voltage ................................................................................... 62
Figure 5.20 Direct and quadrature current component ................................................................. 63
Figure 5.21 Active and Reactive Power injected into Grid .......................................................... 64
Figure 5.22 DC-Link voltage at resistive load of 2.5kW and 3.5kW ........................................... 65
viii | P a g e
Figure 5.23 FFT Analysis of Inverter Current at load of 2.5kW and 3.5kW................................ 65
Figure 5.24 FFT Analysis of Grid Current at load of 2.5kW ....................................................... 66
Figure 5.25 FFT Analysis of Grid Current at load of 3.5kW ....................................................... 66
Figure 5.26 Non-linear Load......................................................................................................... 67
Figure 5.27 DC-Link voltage at non-linear load of 125Ω and 320Ω and C = 3000µF ................ 68
Figure 5.28 FFT Analysis of Inverter Current at non-linear load of 100Ω, 125Ω and 320Ω and C
= 3000µF ....................................................................................................................................... 68
Figure 5.29 FFT Analysis of Grid Current at non-linear load of 100Ω and C = 3000µF ............ 69
Figure 5.30 FFT Analysis of Grid Current at non-linear load of 125Ω and C = 3000µF ............ 69
Figure 5.31 FFT Analysis of Grid Current at non-linear load of 320Ω and C = 3000µF ............ 70
ix | P a g e
List of Abbreviations
PV Photovoltaic
DC Direct current
AC Alternating current
d-q Direct-quadrature
RF Reference frame
x|Page
Grid Connection and Synchronization Inverter Control for Solar Power using SRRF
2019
Chapter One
Introduction
MPPT
Controller
(IC)
PWM
Generation
PWM
Generation
Inverter
Phase Locked
Control Using
Loop(PLL)
SRRF
Chapter Two
Literature Review
The recent work mainly focuses on the grid interfacing behaviour, control performance and
operation after connecting to renewable energy. The integration of renewable energy sources to
the utility grid depends on the scale of power which is being generated. Different researchers
have been worked on the grid connected PV system to improve the quality of power in different
ways as well as to enhance the stability and efficiency of the system by employing various types
of techniques and methods.
A. Angela Sheril and M. Ramesh Babu (Middle-East Journal of Scientific Research 25 (4): 864-
870, 2017) have been demonstrated that the voltage, current and frequency of the grid connected
photovoltaic system could be synchronized and controlled by using phase locked loops and
stationary reference frame technique. But the active and reactive power flow from the output of
inverter to the grid has not been controlled. Since the two powers are not being controlled
independently, the stability and efficiency cannot be improved [6].
Manel Hammami, Rabeh Abbassi and Souad Chebbi (Conference Paper, January, 2014)
proposed that the efficiency of the photovoltaic system is the serious issue in the world, so the
rotating reference frame technique has been applied to control the active and reactive power
depending on the currents in the d-q rotating reference system in order to improve the efficiency
of the system [7]. However, the new proposed method which is going to apply depends on both
voltage and current in the rotating system for controlling purpose.
Bhavesh M. Jesadia and Prof. Indrajit N.Trivedi have been applied instantaneous reactive power
theory as a control strategy. According to this theory during the sunlight system send active
power to grid and same time compensate the reactive power of load and compensate the
harmonics. During no sunlight available system only compensate reactive power of load and
harmonics. In this method, the synchronization control would not be applied to control the phase
angle of the voltage of the grid [8]. So this leads the system to the faulty conditions. However,
the aim of the proposed system is to control the required power of the system in parallel way
using the technique depending on voltage and current.
multiple layers of amorphous silicon cells which leads to higher rates of efficiency (up to 8%) for
these technologies; even though it is quite costly.
generally used to supply certain DC and/or AC electrical loads. In order to operate the load
during sunlight hours as well as night hours fed by PV power, electrical energy storing devices
(batteries) are required, which makes the system costlier. Figure 2.2 depicts about the block
diagram of standalone PV system where PV array is connected to DC and AC load through
charge controller.
Start
Measure Vnew
and Inew
Calculate power,
Pnew
Yes No
Pnew>Pold
Vnew>Vold Vnew>Vold
Yes No Yes No
Start
Sense PV
array V and I
Calculate V
and I
V=0
No Yes
Yes
I/ V=-I/V Yes
I=0
No No
I/ V>-I/V
I>0
Yes
Yes No
No
Update
voltage and
current
Return
components are also contained in the transformation to the arbitrary reference frame. This
transformation is a well-known three-phase to two-phase transformation in synchronous system
analysis [19]. It is presented in the figure below.
The Phase Detector: the phase detector outputs a signal that is proportional to the phase
difference between the reference signal and the signal generated by the transformation
block (abc→ dq).
The Loop Filter: This block presents a low-pass filter, usually PI controller, to attenuate
the high-frequency AC components from the phase detector output. The output of this
block is the grid frequency.
VCO Controller: This block integrates the grid frequency to obtain the phase shift of the
grid voltage (Ө).
Chapter Three
System Modeling of Photovoltaic Inverter
A solar cell is basically a p-n junction which is fabricated in a thin wafer of semiconductor. The
electromagnetic radiation of solar energy can be directly converted to electricity through the
photo voltaic effect. When it is being exposed to the sunlight, the photons with energy can be
greater than the band-gap energy of the semiconductor creates some electron-hole pairs
proportional to the incident irradiation [28]. The equivalent circuit of the photovoltaic cell is as
shown below in figure 3.2.
Where G is solar radiation, Isc is short circuit current of the PV cell, Tref is the nominal
temperature and Ki is the temperature coefficient of short circuit current of the cell at 25oC and
1000w/m2.
From equivalent-circuit, the diode current of the cell can be expressed as:
[ ] (3.2)
(3.3)
Where V is the external voltage and I is the external current of the photovoltaic cell.
VT is thermal voltage and Io is the dark saturation current of the diode.
The dark saturation current of the diode is mathematically expressed as:
( )
( )
( * (3.4)
[ ]
By applying all the above equations, the external current can be found as:
( )
[ ] (3.7)
The Simulink modeling is done using the above equations in the following manner.
The conversion of temperature from degree Celsius to kelvin is shown in figure 3.3.
(3.9)
The DC-link capacitor can be obtained using the following formula [21],
(3.11)
The calculated values of each component used for this design are provided below.
Table 3.3. Boost converter calculated parameters
The conventional MPPT algorithms are using to obtain the maximum power point.
Several algorithms can be used in order to implement the MPPT as follows: perturb & observe,
incremental conductance, parasitic capacitance and constant voltage, but only the first two are
the most frequently used.
The incremental conductance algorithm has been chosen as a MPPT control strategy in this
thesis work. This algorithm has advantages compared to perturb and observe as it can determine
when the MPPT has reached the MPP, where perturb and observe oscillates around the MPP.
The incremental conductance can track rapidly the increase and decrease of irradiance conditions
with higher accuracy than perturb and observe. One disadvantage of this algorithm is the
increased complexity when compared to perturb & observe [35].
Using the P-V curve described in figure 2.5 under chapter 2, the following equation can be
analyzed:
at MPP,
In addition, the following equations can be obtained using eq. (3.12) by applying chain rule for
the derivative of product yields,
(3.13)
The DC-DC converter is controlled by incremental conductance technique in order to achieve the
maximum power point (MPP). The incremental conductance method is based on the observation,
the slope of the PV array curve is zero at MPP In addition to this, the curve is increasing from
left of MPP and decreasing from the right of MPP.
at MPP,
Here is the block diagram of the incremental conductance method used to extract the maximum
power from the PV array.
Where, m is the modulation index, VDC is the DC-link voltage and VLL(rms) is the line to line
voltage which is chosen for this design.
The design procedure for the LCL-filter can be described by the following steps [38-41]:
Step 1: the filter capacitance is selected based on the limit of the reactive power absorbed by the
capacitor. Generally, it should be less than 5% of the rated power. The value of capacitance is
given by [39]:
(3.16)
Where, Cb is the base capacitance and x is the percentage of reactive power absorbed at rated
condition.
The base capacitance of the system can be calculated using the following formula:
(3.17)
(3.18)
Where, VLLrms is the rms value of line to line voltage, Prated is the rated power of the system, Zb is
the base impedance and fg is the grid frequency.
Step 2: the inverter side inductor is used to reduce the ripple of the inverter current and it is
calculated using the relationship between the ripple current and inverter inductor [39]:
(3.19)
Where, ∆I is the ripple current, usually limited from 10% to 20% of rated current, fs is switching
frequency.
The ripple current of the inverter (∆I) can be selected as 20% of base current. The base current
(Ib) of the system is given by:
(3.20)
√
Step 3: the grid side inductor can be found based on the required ripple attenuation factor (RAF)
at switching frequency. RAF is the ratio of the grid current to the inverter output current as
follows:
( )
| | (3.22)
( )
RAF at switching frequency is given by eq. (3.22) which can be rewritten as eq. (3.23) to
calculate the grid inductor.
| | (3.23)
| | (3.24)
Step 4: the resonance frequency should be far away from the grid frequency and switching
frequency. A suitable range of fres is given by the following equations [38]:
√ (3.25)
(3.26)
Step 5: to avoid the resonance phenomenon, a damping resistor (Rf) is usually connected in
series with filter capacitor. The value of the damping resistor can be calculated using the
following formula [42]:
( )
(3.27)
Figure 3.15 depicts the effect of the damping resistor on the transfer function of the LCL-filter
which illustrates the bode plot of the filter transfer function without using damping resistor and
the performance of damping resistor how it prevents the resonance frequency in the filter.
(a) (b)
Figure 3.15 Bode Plot of LCL filter: (a) without Damping Resistor (b) with
Damping Resistor
After calculating the required parameter values of the inverter filter design, it is provided in the
following table.
angular velocity of where the direction of rotation depends on the phase sequence of the
voltages. If it has a positive phase sequence, then it rotates in the counter clockwise direction.
Otherwise, it rotates in the clockwise direction with a negative phase sequence. The three phase
voltages could be described with only two components, α and β, in a two dimensional plane. The
magnitude of each active vector is 2 . The active vectors are apart and describe a hexagon
boundary. The locus of the circle projected by the space reference vector depends on V0, V1, V2,
V3, V4, V5, V6 and V7. According to the principle of space-vector PWM, the alpha-beta (αβ)–
coordinate system can be realized by Clarke transformation using eq. (3.28).
[ ] [ ]
√ √
[ ]
The 3-phase quantities are reduced to two phase quantities, the first one is direct axis (real axis)
and other one is quadrature axis (imaginary axis). These two quantities are represented as the
magnitude (Vref) and angle (α). Vref rotates with the speed of angular velocity (ω= f). The
SVPWM is the combination of the reference voltage and 6 switching states of the inverter. The
basic switching vector and sector diagram is shown in figure 3.16, forms a hexagon which can be
seen as consisting of six sectors spanning 60°each. The reference vector which represents
three-phase sinusoidal voltage is generated using SVPWM by switching between two
nearest active vectors and zero vectors [44].
The magnitude of the reference voltage and angle ( of the above equation used to implement
the SVPWM are calculated as follows,
| | √ (3.29)
( *
Determination of Sector
It is very necessary to understand in which sector the reference output voltage lies in order to
determine the switching time and sequence. The phase voltages corresponding to eight switching
states, six nonzero vectors and two zero vectors at the origin. Depending on the reference voltage
and , the angle of the reference vector can be used to determine the sector as shown in
Table 3.5.
Table 3.5. Sector Identification
Sector Degrees
1
2
3
4
5
6
adjacent basic vectors differ in only one bit, so that only one of the upper transistors switches is
closed when the switching pattern moves from one vector to an adjacent one. The two vectors
are time weighted in a sample period to produce the desired output voltage.
To calculate the time of application of different vectors, consider figure 3.17 which is depicting
the position of different available space vectors and the reference vector in the first sector.
√ ( )
√ ( ) (3.31)
Space vector Switching state On state switch Off state switch Vector definition
V0 000 S4, S6, S2 S1, S3, S5
V1 100 S1, S6, S2 S4, S3, S5
Chapter Four
Grid Connected PV System
The Distributed power generation systems (DPGS) can be broadly classified into two categories
which are standalone and grid connected PV systems. In standalone PV systems, storage
batteries are used which stores the energy obtained from these intermittent sources. These
batteries serve as input to the inverter which converts DC power into AC power.
In the grid connected PV systems, energy storage devices may not be required as the energy
obtained from the sun light directly. In the case of a solar PV system, the sun light can be
directly converted to AC and transferred to the grid. To achieve this, the inverter must be
controlled in advanced way using control loops so that the expected power can be transferred to
the grid and also the inverter‟s output terminals cannot directly be connected to the grid because
of the harmonics that are available in the current due to the switching nature of the inverter. A
filter is used to reduce the total harmonic distortions (THD) from the inverter and grid current.
synchronization to power flow management and pulse width modulation of the inverter is taken
care by the control structure of the inverter.
In this model, the control going to be used is synchronously rotating reference frame which is
used to transform voltage and current quantities from natural reference frame (abc-frame) to
synchronously rotating reference frame (dq-frame).
√ √
[ ] [ ] (4.1)
[ ]
√ √
[ ] [ ] (4.2)
[ ]
The scaling factor 2/3 is being considered for amplitude invariant transformations. In a balanced
three phase systems, Ia+Ib+Ic=0 and I0 can be neglected, which is resulting in simplified
transformation form as follows,
[ ] [ √ √
] [ ] (4.3)
From the three phase currents, the α-β current components can be calculated as can be seen in eq.
(4.3). The α-β components of the voltages can be calculated from line to line voltages using eq.
(4.4).
[ ] [ √ √
] [ ] (4.4)
The simplified inverse Clarke transformation for converting the quantities back to natural
reference frame (abc-frame) is given as
[ ] √
[ ] (4.5)
[ √ ]
[ ] √
[ ] (4.6)
[ √ ]
( ) ( )
[ ] ( ) ( ) [ ] (4.7)
[ ]
( ) ( )
[ ] ( ) ( ) [ ] (4.8)
[ ]
Where θ=ωt is the reference phase angle. This transformation requires the information of the
reference phase angle. The phase angle information can be obtained by using the phase locked
loop (PLL).
Since the system is balanced, the above eq. (4.7) and eq. (4.8) are simplified in the following
transformation matrix form as follows,
( ) ( )
[ ] [ ] [ ] (4.9)
( ) ( )
( ) ( )
* + [ ] [ ] (4.10)
( ) ( )
[ ] [ ( ) ( )] [ ] (4.11)
( ) ( )
[ ] [ ( ) ( )] * + (4.12)
( ) ( )
In this thesis work, the dq0 to abc transformation helps the space vector pulse width modulation
(SVPWM) to generate the pulse gates for the switching nature of the inverter in order to control
the voltage source inverter (VSI).
(4.13)
(4.14)
Where Kp and Ki are proportional gain, C is the boost converter capacitor and Ts is the switching
period of the inverter switch.
Current Control Loop
This controller is used to regulate the injected current in synchronously dq-reference frame. The
PI controller gain is selected using the formula,
(4.15)
(4.16)
Where L is the filter inductor, R is the filter resistance which should be very minimum that is
1ohm‟s and is the time constant given by L/R.
The external voltage loop controller is used to keep the DC-link voltage constant. The output of
the active current reference is generated by the DC-link voltage controller.
The reactive power controller sets the reactive current reference which is set to zero. The main
function of the voltage loop controller is to change the active reference current so that the power
available from the PV array can be matched with the power injected to the grid. Figure 4.6 shows
the external voltage loop controller.
The internal current controller loop is used to regulate the grid currents that are converted into
direct current (Id) and quadrature current (Iq). The internal current loop controller is worked
based on the active current reference and reactive current reference. But, the reactive power
current reference must be taken as zero in order to supply only active power to the grid by
providing the desired voltage signal to achieve the unity power factor. Figure 4.7 shows the
internal current loop controller.
Figure 4.8 Overall Voltage and Current controller for Active and Reactive Power Control
Chapter Five
Simulation Results and Discussion
5.1. Introduction
In this chapter, the results obtained from the model which has been designed in chapter 3 and
chapter 4 are presented. This shows all the simulations for the PV array with the characteristics
curve, DC-DC boost converter with MPPT controller, reduction of THD of the inverter and the
inverter control scheme and the results of the overall system under different atmospheric
conditions. The complete grid connected PV model and simulation results were being done using
MATLAB/SIMULINK.
5.2. PV Array
The results of the designed model have similarity with that of the results of the PV array which
are provided in the datasheet in table 3.1. The graphs of I-V and P-V characteristics curve of the
PV array are obtained which are quite similar to the actual behavior of a real PV array under the
effect of various levels of irradiance and temperature which are depicted in figure 5.1 and figure
5.2. The MATLAB codes used to find the PV array characteristics are available in appendix A.
Figure 5.2 Effect of solar irradiance changes on I-V and P-V Curve
It is clear that the simulation results of the modeled PV module match with the PV module in the
datasheet by comparing the designed model results with the datasheet results.
Figure 5.3 I-V and P-V Curve at different temperature found from datasheet
Figure 5.4 I-V and P-V Curve at different irradiance found from datasheet
The sun power SPR-305 module has a maximum power of 305W under standard test conditions.
The PV array is constructed using a PV string which consists of five module connected in series
in order to generate a power of around 1.5KW as shown in figure 5.5. The measured
characteristics curves of I-V and P-V for PV array using MATLAB/Simulink under standard test
conditions are shown in figure 5.6.
Figure 5.5 PV array characteristic curve for I-V and P-V under STC
Figure 5.6 I-V and P-V Curve for PV Array Measured from MATLAB/Simulink
The simulation result of the system model in chapter 3 is implemented in Simulink environment.
The simulation running time is 2.5 second. The input signals of the PV array can be either
constant or time-varying irradiance and temperature which are illustrated in figure 5.7. However,
the system is operated at different solar irradiance and temperature level.
The magenta color indicates the operating temperature that has a value of 25oC and 50oC
respectively.
The red color shows the operating solar irradiation that has a value of 1000W/m2, 250W/m2 and
again 1000W/m2 respectively.
Figure 5.10 Error Signal Between DC-link Voltage and Reference Voltage
Figure 5.12 Output of Inverter before and after filter for Current waveform
Figure 5.13 Output of Inverter before and after filter for Voltage waveform
Figure 5.16 Output Waveform of Three Phase Voltage for Inverter and Grid
Figure 5.17 illustrates the results of the three phase current which are injected into the inverter
and grid. As can be seen in the figure 5.17, the three phase grid voltage is in phase with that of
the three phase inverter current. The magnitude of current of each phase is around 2.142A rms.
The magnitudes of harmonics available in the injected current with respect to the fundamental
frequency/harmonic order for phase A and the total harmonics distortion (THD) of the injected
grid current and voltage, and Fast Fourier Transform (FFT) analysis of the grid current and
volatge as frequency spectrum is shown in figure 5.18 and figure 5.19.
Figure 5.17 Output Waveform of Three Phase Current for Inverter and Grid
Figure 5.23 FFT Analysis of Inverter Current at load of 2.5kW and 3.5kW
Resistive Load THD of Inverter Current (%) THD of Grid Current (%)
1.5kW 0.86 0.88
2.5kW 0.87 0.91
3.5kW 0.87 0.93
Figure 5.27 DC-Link voltage at non-linear load of 125Ω and 320Ω and C = 3000µF
Figure 5.28 FFT Analysis of Inverter Current at non-linear load of 100Ω, 125Ω and 320Ω and C
= 3000µF
Figure 5.29 FFT Analysis of Grid Current at non-linear load of 100Ω and C = 3000µF
Figure 5.30 FFT Analysis of Grid Current at non-linear load of 125Ω and C = 3000µF
Figure 5.31 FFT Analysis of Grid Current at non-linear load of 320Ω and C = 3000µF
Table 5.2. A summary for different values of nonlinear load
Non-linear load and THD of Inverter Current (%) THD of Grid Current (%)
C = 3000µF
100-Ω 0.12 0.61
125-Ω 0.12 0.51
320-Ω 0.12 0.50
Chapter Six
Conclusion and Future Work
6.1. Conclusion
In this thesis work, the task of designing and modelling of the three phase grid connected solar
PV system using synchronously rotating reference frame has been successfully carried out. The
proposed system comprises of the PV array, DC-DC boost converter, Incremental conductance
based MPPT, 3-ɸ VSI, LCL-filter.
The PV array was modelled accurately that was constructed using MATLAB/SIMULINK
environment. The accuracy of the model for PV array was verified by comparing the results of
the model with the datasheet information. The DC-DC boost converter was designed and used to
boost the output voltage of the PV array to the desired level. The control algorithms have been
proposed in order to maximize the output power of the PV array. The most important parts of the
system are the MPPT, DC-link regulator, active and reactive power control, total harmonics
eliminations and synchronization control. Simulation results show that the incremental
conductance based MPPT method tracks the maximum power point of the PV array to extract
maximum power for a given environmental condition. A three phase voltage source inverter was
controlled using SRRF with SVPWM to convert the DC voltage to an AC voltage in which its
output current is synchronized with the grid voltage. An LCL-filter was designed which used to
reduce and attenuate the high frequency harmonics at the output side of the inverter. This filter is
used to interconnect the inverter output with the grid which provides the filtering of harmonics
available in the current due to the switching nature of the inverter.
The control structure of the inverter is done based on two loops such as external loop and
internal loop which are used to perform as DC-link voltage control and current control loop in
order to maintain the power flow in the system and power quality and harmonics present in the
current respectively. The simulation results show that the proposed system operates efficiently
and effectively under various weather conditions. The efficiency of the overall system is around
97% which is improved more than the previous related work whose efficiency is around 96%.
The total harmonic distortion of the grid current is 0.88% which is less than the previous work
that is 1.42%. For resistive load change, the THD of the grid current and inverter current are
0.88%, 0.91% and 0.93%, and 0.86%, 0.87% and 0.87% for a load of 1.5kW, 2.5kW and 3.5kW
respectively. For a nonlinear load change, the THD of the grid current and inverter current are
0.61%, 0.51% and 0.50%, and 0.12%, 0.12% and 0.12% for a load of 100Ω, 125Ω and 320Ω and
C = 3000µF respectively.
References
[1] V. Katinas, M. Marčiukaitis and M. Tamašauskienė, “Analysis of the wind turbine noise
emissions and impact on the environment”, Renewable and Sustainable Energy Reviews, vol. 58,
pp. 825-831, 2016.
[2] I. Onakpoya, J. O‟Sullivan, M. Thompson and Heneghan, “The effect of wind turbine noise
on sleep and quality of life: A systematic review and meta-analysis of observational studies”,
Environment International, vol. 82, pp. 1-9, 2015.
[3] R. Fu et al., “U.S. Solar Photovoltaic System Cost Benchmark: Q1 2016” ,National
Renewable Energy Laboratory, 2016.
[4] Marcelo Gustavo Molina,”Modeling and Control of Grid-connected Solar Photovoltaic
Systems”, Institute of Electrical Energy, National University of San Juan−CONICET, Argentina.
[5] V. Jagadeesh Babu and J. Arun Venkatesh, ”Control of Grid Connected Hybrid Energy
System using Synchronous Reference Frame Technique”, International Journal of Advanced
Research in Electrical, Electronics and Instrumentation Engineering, Vol. 3, Special Issue 4,
May 2014.
[6] A. Angela Sheril and M. Ramesh Babu,”Synchronization Control of Grid Connected
Photovoltaic System”, Middle-East Journal of Scientific Research 25 (4): 864-870, 2017.
[7] Manel Hammami, Rabeh Abbassi and Souad Chebbi, “Design and simulation of a grid-
connected photovoltaic system” Conference Paper, January, 2014.
[8] Bhavesh M. Jesadia and Prof. Indrajit N.Trivedi,”Analysis Of PV System Embeded To
Distribution Grid For Active & Reactive Power Supply To Grid”, International Journal of
Advance Engineering and Research Development, Volume 2, Issue 5, May, 2015.
[9] Tomas Markvart (2001, „‟Solar Electricity‟‟, John wiley & sons Baffin‟s Lane chichester
west Sussex P019IUD, England,p-25.
[10] Green et al., “Solar cell efficiency tables (version 49),” Progress in
Photovoltaics: Research and Applications, vol. 25, no. 1. pp. 3–13, Jan. 2017.
[11] S. Sumathi, L. Ashok Kumar, and P. Surekha, „‟Solar PV and Wind Energy Conversion
Sytems‟‟.
[12] M. A. Elgendy, B. Zahawi and D. J. Atkinson, "Evaluation of perturb and observe MPPT
algorithm implementation techniques," 6th IET International Conference on Power Electronics,
Machines and Drives (PEMD 2012), Bristol, pp. 1-6, 2012.
[13] M. Eltawil and Z. Zhao, "MPPT techniques for photovoltaic applications", Renewable and
Sustainable Energy Reviews, vol. 25, pp. 793-813, 2013.
[14] M. G. G. Villalva et al., “Comprehensive Approach to Modeling and Simulation of
Photovoltaic Arrays,” IEEE Transactions on Power Electronics, vol. 24, no. 5. pp.
1198–1208, May 2009.
[15] Rashid, Power Electronics: Circuits, Devices, and Applications, 3rd ed. Prentice
Hall, 2003.
[16] Mohan et al., Power Electronics: Converters, Applications and Design, 3rd ed.
John Wiley & Sons, Ltd, 2003.
[17] J. Jana, H. Saha and Das Bhattacharya, "A review of inverter topologies for single-phase
grid-connected photovoltaic systems", Renewable and Sustainable Energy Reviews, vol. 72, pp.
1256-1270, 2017.
[18] F. Blaabjerg et al., “Overview of Control and Grid Synchronizations of DPG Systems,”
IEEE Transactions on Industrial Electronics, vol.53, no. 5. pp. 1398–1409, 2006.
[19] Digital Signal Processing Based Electromechanical Motion Control, 2004.
[20] L. N. Arruda, S. M. Silva “Phase Locked Loop structures for utility Grid connected
systems,” in Proc. IEEE-IAS Annu. Meeting, vol. 4, pp. 2655–2660, 2002.
[21] S. Kjaer, J. Pedersen and Blaabjerg, "A Review of Single-Phase Grid-Connected Inverters
for Solar Photovoltaic Systems", IEEE Transactions on Industry Applications, vol. 41, no. 5, pp.
1292-1306, 2005.
[22] IEC 61727. Solar Photovoltaic systems – characteristics of the utility interface 2004.
[23] IEEE 1547. Standards for Interconnecting DRs with Electric Power Systems, 2003.
[24] IEEE Standard 929. IEEE Recommended Practice for Utility Interface of solar PV Systems,
2000.
[25] M. Prodanovic and T. Green, "Control and filter design of three-phase inverters for high
power quality grid connection," in IEEE Transactions on Power Electronics, vol. 18, no. 1, pp.
373 380, 2003.
[26] M. Rashid, Power electronics handbook, 3rd ed. Elsevier, 2011.
[27] G. M. Masters, Renewable and Efficient Electric Power Systems. Hoboken, NJ,
USA: John Wiley & Sons, Inc., 2004.
[28] Samlexsolar, “Solar PV Cell Module, array”.
[29] Xuan H. Nguyen & Minh Ph. Nguyen, „‟Mathematical Modeling of PV Cell, PV Module
and PV Array‟‟, 2015.
[30] Omar Mohammed Benaissa, Samir Hadjeri and Sid Ahmed Zidi,”Modeling and Simulation
of Grid Connected PV Generation System using Matlab/Simulink‟‟, 2017.
[31] Liu X, Wang P & Loh PC (2010), „‟Optimal Coordination Control for Standalone PV
System with Nonlinear Load‟‟, Proceedings of 2010 9th international power and energy
conference, Singapore, 27-29October 2010, 5pp.
[32] Grandi G, Rossi C, Ostojic D et al (2009) A new multilevel conversion structure for grid-
connected PV applications. IEEE Trans Ind Electron 56(11):4416–4426.
[33] Mohan N, Under land TM (2007) Power electronics: converters, applications, and design.
Wiley, New York.
[34] F. Gao, D. Li, P. C. Loh, Yi Tang and Peng Wang, "Indirect dc-link voltage control of two
stage single-phase PV inverter," 2009 IEEE Energy Conversion Congress and Exposition, San
Jose, CA, pp. 1166-1172, 2009.
[35] A.A.Nouredine and A.A. Yacine,” simulation and implementation of single-stage three
phase grid connected PV system using SVM-controlled three-level inverter‟‟, 2017.
[36] K. Zhou and Wang, “Relationship between space-vector modulation and three-phase
carrier-based PWM: a comprehensive analysis [three-phase inverters],” IEEE Transactions on
Industrial Electronics, vol. 49, no. 1. pp. 186–196, 2002.
[37] W.F.Zhang et al., “Comparison of Three SVPWM Strategies,” Journal of
Electronic Science and Technology, vol. 5, no. 3. Editorial Board of Journal of
University of Electronic Science and Technology, pp. 283–287, 2007.
[38] M.-Y.Park et al., “an LCL filter design for grid-connected PCS using THD and RAF,” The
2010 International Power Electronics Conference - ECCE ASIA -. IEEE, pp. 1688–1694, Jun.
2010.
[39] X. Wei et al., “Design of LCL filter for wind power inverter,” in 2010 World NonGrid-
Connected Wind Power and Energy Conference, 2010, pp. 1-6.
[44] Suresh Lakhimsetty, “Simulation of space vector PWM for voltage source inverter using
matlab/Simulink‟‟, 2014.
[45] D. Candusso, L. Valero, A. Walter, S. Bacha, E. Rulliere and B. Raison, "Modeling, control
and simulation of a fuel cell based on power supply system with energy management," IEEE
2002 28th Annual Conference of the Industrial Electronics Society. IECON 02, pp. 1294-1299
vol.2, 2002.
[46] C. Ramos, A. Martins, and Carvalho, “Current control in the grid connection of the double
output induction generator linked to a variable speed wind turbine,” in Proc. IEEE IECON,
vol.2, pp. 979–984, 2002.
[47] T. Huang, X. Shi, Y. Sun and D. Wang, "Three-phase photovoltaic grid-connected inverter
based on feed-forward decoupling control," 2013 International Conference on Materials for
Renewable Energy and Environment, Chengdu. pp. 476-480, 2013.
[48] Matlab/Simulink from mathwork for 2018b.
[49] A. Vijayakumari, A. Devarajan and N. Devarajan, "Decoupled control of grid connected
inverter with dynamic online grid impedance measurements for micro-grid applications",
International Journal of Electrical Power & Energy Systems, vol. 68, pp. 1-14, 2015.
[50] L. Ma, W. Ran and T. Q. Zheng, "Modeling and control of three phase grid-connected
photovoltaic inverter," Control and Automation (ICCA), 2010 8th IEEE International
Conference on, Xiamen, 2010, pp. 2240-2245.
Appendices
Appendix A: MATLAB CODES to obtain I-V and P-V Curve
MATLAB code to obtain the PV characteristics curves for I-V and P-V curve
%constant and variables
q=1.6*10^-19; %electron charge
Ego=1.1; %band-gap energy of semiconductor
K=1.380504*10^-23; %Boltzmann constant
Ki=0.061745; %Temperature coefficient of short-circuit current
A=0.94504; %ideality factor
Isc=5.96; %short-circuit current
Voc=64.2; %open-circuit voltage
Rs=0.37154; %series resistance
Rsh=269.5934; %shunt resistance
Ns=96; %number of cells connected in series
Np=1; %number of cells connected in parallel
Tref=298; %nominal temperature
Tk= [100 75 50 25]; % different temperatures
Top=273+Tk; %operating temperature
G=1000; %solar irradiance
Vpv= (0:70); %output voltage of PV system
for i=1:4
Iph=(Isc+Ki*(Top(i)-Tref))*G/1000;
Irs=Isc/(exp((q*Voc)/(A*Ns*K*Top(i))-1));
Io=Irs*(Top(i)/Tref )^3*(exp((q*Ego*((1/Tref)-(1/Top(i))))/(A*K)));
Ish=Vpv/Rsh;
Ipv=Iph-Io*(exp((q*Vpv)/(A*Ns*K*Top(i))-1)-Ish);
Ppv=Vpv.*Ipv;
figure(1)
plot(Vpv, Ppv,'linewidth',2)
title('P-V Curve for PV Module')
82 | P a g e