Photovoltaic Power Generation
Photovoltaic Power Generation
q Vdc
I pv N p I ph N p I rs exp 1
kTc A N s (2)
Ppv Vdc *I pv f (Vdc , S ,Tc )
IV and PV characteristic
• Open-circuit voltage (Voc):
• Terminal voltage when the terminals of the module
are disconnected.
• Short circuit current (Isc)
• Maximum Power PMPP where the voltage versus current
product is maximum
• VMP≈ 0.8 * Voc
• IMP≈ 0.9 * ISC
• The best conditions, are the "standard test
conditions " happen at Irradiance equal to 1000W/m²,
cells temperature equals to 25°C, and spectral
distribution (Air Mass) AM is equal to 1.5
Maximum Power Point Tracking
• Basic methods include:
• Fractional Open Circuit Voltage
Vmpp ≈ K * Voc
• Fractional Short Circuit Current
Impp ≈ K * Isc
• Look-up Table Technique
• Perturb and Observe
• Incremental Conductance
• Fuzzy Logic
• Neural network based
P&O MPPT algorithm and working principle
• Here dc voltage reference (Vdcref) is perturbed with a fixed amount (∆V) and the next perturb is
decided based on the status of change in power (∆P).
• For correct execution of the P&O MPPT, fMPPT should not be greater than the inverse of 20ms
(settling time).
• The sensing of Ipv and Vpv done at higher frequency than the frequency of updating (the Vdcref),
fMPPT may result in erroneous results since the decision of perturb is taken based on an unsettled
outputs as it may destabilize the system.
702
Start
701
in V
Sensing/Sampling of
dc
Vdc(n) and Ipv (n) 700
V
Vdc Ref
699
Vdc
P ( n) V ( n) * I ( n)
pv dc pv 1.95 1.96 1.97 1.98 1.99 2 2.01 2.02 2.03 2.04 2.05
Time in second
dP ( n ) P ( n ) P ( n 1) (b)
pv pv pv
dV ( n ) V ( n ) V ( n 1)
dc dc dc
VMPP2
NO YES
dPpv ( n )
0
dVdc ( n )
k = -1 k=1
VMPP1
Vdcref ( n 1) Vdc ( n ) k V 1 2 N Samples
Time
(a) (c)
Three step operation in P&O
PB>PC
PB<PC
40 40
B MPP
C B MPP
C
30
VA →VB →VC →VB →VA A
VB →VC →VD →VC →VB D
30
in kW
in kW
20 20
pv
pv
P
P
10 VA VC 10 VB
VD
0 0
0 100 200 300 400 500 600 700 800 900 0 100 200 300 400 500 600 700 800 900
Voltage in V VB VMPP VC
Voltage in V VMPP
Fig. 1 Fig. 2
From fig 1 the dP/dV is greater than zero and results in positive perturb (since PA<PB) when the operation is shifted
either from A to B or from B to A. Similarly, when it is from B to C or from C to B the dP/dV is less than zero and hence
the P&O algorithm initiates a negative perturb (since PB>PC). Hence, it can be concluded that a well executed P&O
algorithm contributes continuously shifting of dc voltage, involving three levels. Similarly for the fig 2 can be
explained.
Incremental conductance
• At the MPPT point of the P~V curve
dp dv di di i
0i v 0
dv dv dv dv v
Comparison of MPPT algorithms
Impact of parallel leakage resistance on IV characteristics
• In ideal PV cell
• 𝐼𝑃𝑉 = 𝐼𝑠𝑐 − 𝐼𝑑𝑖𝑜𝑑𝑒 ; where 𝐼𝑠𝑐 is the short
circuit current and assume 𝑅𝑠 = 0 and 𝑅𝑃 =
∞
• When parallel leakage resistance is
considered
𝑉𝑑𝑐
• 𝐼𝑃𝑉 = 𝐼𝑠𝑐 − 𝐼𝑑𝑖𝑜𝑑𝑒 − ; keeping 𝑅𝑆 = 0.
𝑅𝑃
• This shows that at any given voltage load
𝑉
current will be decreased by ( 𝑑𝑐 ) from the
𝑅𝑃
load current obtained with ideal PV cell
𝑉
• Due to this reduction of current (equal to 𝑑𝑐 ),
𝑅𝑃
1
the IV characteristics shows a slope of ; as
𝑅𝑃
shown in figure.
• For a cell to have less than 1% loss due to
100𝑉𝑜𝑐
the parallel resistance, 𝑅𝑃 >
𝐼𝑠𝑐
Impact of series resistance on IV characteristics
• Series resistance includes
• Contact resistance associated with the bond between the cell and its wire leads
• Resistance of the semiconductor
• In ideal PV cell
• 𝐼𝑃𝑉 = 𝐼𝑠𝑐 − 𝐼𝑑𝑖𝑜𝑑𝑒 ; Assume 𝑅𝑠 = 0 and 𝑅𝑃 = ∞
𝑞𝑉
• Diode current is given by 𝐼𝑑𝑖𝑜𝑑𝑒 = 𝐼𝑟𝑠 (exp − 1)
𝑘𝑇𝑐 𝐴
• Let n be the number of cells connected in series and one of the cells is
partially shaded
• The characteristics will depend on the short circuit current of the partially
shaded cell
• For example if one of the n cells in 50% shaded, the short circuit current of
𝐼𝑆𝐶
that cell reduces to
2
• The IV characteristics is obtained by considering two operating
circumstances
𝐼𝑆𝐶
• First, when current supplied by remaining n-1 cells is less than
2
𝐼𝑆𝐶
• Second, when current supplied by remaining n-1 cells is greater than
2
Impact of Shading
Scenario of one of the series connected cells getting partially shaded (cont.)
𝐼𝑆𝐶
1. When I < (Figure (a))
2
• Diode of the 1st cell is forward biased
𝐼 𝑉
• Current 𝐼𝑑 = ( 𝑆𝐶 − 𝐼 − 𝑑 ) will flow through
2 𝑅𝑃
the diode
• The partially shaded cell contributes positive
voltage
𝐼𝑆𝐶
2. When I > (Figure (b))
2
𝐼
• A current equal to (I − 𝑆𝐶 ) will flow through
2
the parallel resistance
• There is a voltage drop across the resistance.
• The diode becomes reverse biased and carries
no current
• Similar to a complete shading scenario, there
𝑉
is a resultant voltage drop given by ∆𝑉 = +
𝐼𝑆𝐶 𝑛
(𝐼 − )𝑅𝑃
2
Impact of Shading
Change in IV characteristics for various partial shading scenarios
(a) Three bypass diodes, each covering (b) IV characteristics and maximum power (c) Use of bypass diode in a string of modules
one – third of the cells in a module of a module with 3 bypass diodes and one
cell shaded
Shade mitigation using blocking diode
• When strings are connected in parallel
• Instead of supplying current, the shaded string can withdraw current from rest of the parallel
connected strings (Figure (a))
• This problem is solved by using blocking diode; also known as isolation diode
• The blocking diode is placed at top of each string
• The diode blocks the reverse current withdrawn by the shaded string (Figure (b)).
c) IV and PV curves for a PV module under various shading scenarios as in (b) and
each cell having bypass diode connected across it.
Grid connected PV module
• Controlling action
• Active and reactive power control
• DC link voltage control
• MPPT tracking for variable insolation
• Controller configuration
• Single stage inverter
• Inverter handles MPPT, P&Q control
• Double stage inverter
• Boost converter handles MPPT, Inverter handles P&Q control
Decoupled controller philosophy
• The active and reactive power controlled by a Decoupled controller
• For a three phase system total power in dq frame
3 3
𝑆 = 𝑉𝑑 + 𝑖𝑉𝑞 𝐼𝑑 − 𝑖𝐼𝑞 = 𝑉𝑑 𝐼𝑑 + 𝑉𝑞 𝐼𝑞 + 𝑖 𝑉𝑞 𝐼𝑑 − 𝑉𝑑 𝐼𝑞
2 2
• When it is a balanced system, 𝑉𝑞 = 0;
S P iQ
3
2
(Vd I d iVd I q )
• With park’s transformation on the three phase Ipcc and Vpcc its
corresponding Id+iIq and Vd+iVq is found, which needs the information of
the angle and frequency of the supply voltage. This can be achieved
using a PLL(phase locked loop) circuit.
Park’s Transformation
2𝜋 2𝜋
cos 𝜃 cos 𝜃 − cos 𝜃 +
3 3
2 2𝜋 2𝜋
− sin 𝜃 − sin 𝜃 − − sin 𝜃 +
3 3 3
1 1 1
2 2 2
• In this system, the dc-link voltage controller and
reactive-power controller determine reference
components Idref and Iqref respectively.
&
Sliding Mode Control
• Consider the grid connected PV system
• Current equation in dq reference frame is written as
𝑑𝑖 𝑉
𝐿𝑓 𝑑𝑡 = −𝑅𝑓 𝑖𝑑𝑡 + 𝜔𝐿𝑓 𝑖𝑞𝑡 + 𝑣𝑑𝑡 − 𝑣𝑑𝑝𝑐𝑐 ; 𝑣𝑑𝑡 = 𝑚𝑑 𝑑𝑐
𝑑𝑡 2
𝑑𝑖𝑞𝑡 𝑉𝑑𝑐
𝐿𝑓 = −𝑅𝑓 𝑖𝑞𝑡 − 𝜔𝐿𝑓 𝑖𝑑𝑡 + 𝑣𝑞𝑡 − 𝑣𝑞𝑝𝑐𝑐 ; 𝑣𝑞𝑡 = 𝑚𝑞
𝑑𝑡 2
𝑑𝑉𝑑𝑐 1
• By KCL at DC link capacitor = 𝑖𝑃𝑉 − 𝑖𝑑𝑐
𝑑𝑡 𝐶
• Writing these 3 equations in matrix form for reference values
𝑑𝑖𝑑𝑡 𝑅𝑓 𝑖𝑑𝑡 𝑣𝑑𝑝𝑐𝑐 𝑉𝑑𝑐
− + 𝜔𝑖𝑞𝑡 − 0 0
𝑑𝑡 𝐿𝑓 𝐿𝑓 2𝐿𝑓
𝑚𝑑
𝑑𝑖𝑞𝑡 𝑅 𝑖 𝑣 𝑉𝑑𝑐 𝑚𝑞
= − 𝑓 𝑞𝑡 − 𝜔𝑖𝑑𝑡 − 𝑞𝑝𝑐𝑐 + 0 0
𝑑𝑡 𝐿𝑓 𝐿𝑓 2𝐿𝑓 𝑖𝑑𝑐𝑟𝑒𝑓
𝑑𝑉𝑑𝑐 𝑖𝑃𝑉 1
𝑑𝑡 0 0 −
𝑑𝑥
𝐶 𝐶
• or = 𝑓 𝑥 + 𝑔 𝑥 𝑢
𝑑𝑡
Sliding Mode Control (cont.)
𝑑𝑥
• The control input matrix 𝑢 is obtained by 𝑢 = 𝑔 𝑥 −1 − 𝑔 𝑥 −1 𝑓 𝑥
𝑑𝑡
𝑖𝑑𝑡
• The output matrix considered is 𝑦 = 𝑖𝑞𝑡
𝑉𝑑𝑐
• In case of SMC we design a sliding surface so that the control input forces the
instantaneous states to the sliding surface which is achieved by adding a switching
control to the equivalent control i.e. 𝑢 𝑡 = 𝑢𝑒𝑞 𝑡 + 𝑢𝑠𝑤 𝑡 = 𝑢𝑒𝑞 𝑡 −
𝜌 𝑠𝑖𝑔𝑛𝑢𝑚(𝜎); where 𝑢𝑒𝑞 𝑡 = −𝑔 𝑥 −1 𝑓 𝑥 and 𝜌 is a positive constant.
𝜎1 𝑖𝑑𝑡 − 𝑖𝑑𝑡𝑟𝑒𝑓
• The sliding surface followed by the output is 𝜎 = 𝜎2 = 𝑖𝑞𝑡 − 𝑖𝑞𝑡𝑟𝑒𝑓
𝜎3 𝑉𝑑𝑐 − 𝑉𝑑𝑐𝑟𝑒𝑓
• The stability of the controller can be analysed by Lyapunov’s method.
• Another function that is used to avoid hard switching of signum function is hyperbolic
tangent function i.e. 𝑢 𝑡 = 𝑢𝑒𝑞 𝑡 + 𝑢𝑠𝑤 𝑡 = 𝑢𝑒𝑞 𝑡 − 𝜌 𝑡𝑎𝑛ℎ(𝜎)
Results with signum and hyperbolic tangent function
0.8
d-axis current
0.6
0.4
1 2 3 4 5 6 7 8
Time in second
d-axis current
1
0.5
0
1 2 3 4 5 6 7 8
Time in second
0.5
0
1 2 3 4 5 6 7 8
Time in second
(c) System response with hyperbolic tangent switching function
Sliding Mode Control (cont.)
• The reference signals are obtained in the following manner
• 𝑉𝑑𝑐𝑟𝑒𝑓 is obtained from MPPT algorithm
• 𝑖𝑑𝑐𝑟𝑒𝑓 is obtained from the corresponding differential equation
• The 𝑖𝑑𝑡𝑟𝑒𝑓 and 𝑖𝑞𝑡𝑟𝑒𝑓 current references are obtained as per the phasor diagram
• From the phasor diagram
𝑉𝑑𝑐 𝑖𝑑𝑐𝑟𝑒𝑓
𝑖𝑑𝑡𝑟𝑒𝑓 cos 𝜃 sin 𝜃 𝑣𝑝𝑐𝑐
=
𝑖𝑞𝑡𝑟𝑒𝑓 sin 𝜃 − cos 𝜃 𝑄𝑟𝑒𝑓
𝑣𝑝𝑐𝑐
3 phase,
DC
480 V, 60
Hz DC
Circuit
Breaker Multiple PV
loads Array Battery
GATE and its
Grid Side Converter controller
3
𝑃= 𝑣 𝑖 + 𝑣𝑞 𝑖𝑞
2 𝑑𝑑
𝑑𝑖𝑑
𝑅𝑔 𝑖𝑑 + 𝐿𝑔 − 𝜔𝐿𝑔 𝑖𝑞 + 𝑣𝑑 = 𝑣𝑑𝑔
𝑑𝑡
𝑑𝑖𝑞
𝑅𝑔 𝑖𝑞 + 𝐿𝑔 + 𝜔𝐿𝑔 𝑖𝑑 + 𝑣𝑞 = 𝑣𝑞𝑔
𝑑𝑡
𝑑𝑉𝑑𝑐
𝐼𝑑𝑐 =𝐶 + 𝐼𝐿 + 𝐼𝑖𝑛𝑣
𝑑𝑡
I dref Ki 1 1 Id
+ Kp R g L g .s
-
s 1 T a .s
PI Controller PWM Block System transfer
Function
(b)
(a)
(c)
Nominal behaviour of battery in autonomous mode
(a) Solar Irradiation (b) DC-Link Voltage
(c) Battery SOC
Results
Grid connected mode of operation: SOC > SOCref
(b) (c)
(a)
(d)
Power Balance at(e) higher SOC (f)
(a) (b)
• The customer may withdraw power from the grid if the local loads are more
than the PV power generation.
f(p.u.)
participating in frequency control PV participation