Benseddik Djaloud
Benseddik Djaloud
Benseddik Djaloud
Mmoire de Master
En vue de lobtention du diplme de master
Prsent par :
Othmane BENSEDDIK Fathi DJALOUD
THEME
Soutenu le : 27 /06/2012
Devant le jury compos de :
2011/2012
Remerciement
Avant tout nous tenons nos remerciements notre dieu de nos avoir
donn la force et le courage.
Nos plus chaleureux remerciements pour tous ceux qui de prs et de loin
ont contribu la ralisation de cette mmoire.
Introduction gnrale .............................. 1
I
2.5 Simulation du l hacheur .. 27
2.5.1 hacheur srie . 27
2.5.2 Hacheur parallle .................................. 31
2.6 Conclusion ... 33
Conclusion gnrale. 61
Bibliographique ... 63
Annexe . 65
II
Liste des tableaux
IV
Liste des figures
V
Figure 2.18 La valeur moyenne de la tension aux bornes de la charge .. 32
Figure 2.19 La valeur moyenne de courant 33
VI
Introduction gnrale
1.1 Introduction
L'nergie solaire photovoltaque provient de la transformation directe d'une partie du
rayonnement solaire en nergie lectrique. Cette conversion d'nergie s'effectue par le biais
d'une cellule dite photovoltaque (PV) base sur un phnomne physique appel effet
photovoltaque. Ce chapitre prsente les concepts dont la connaissance est ncessaire la
comprhension du fonctionnement des cellules photovoltaques constitues de semi-
conducteur en silicium.
On commencera par brve rappelle sur le principe de la conversion de l'nergie solaire en
nergie lectrique. Nous prsenterons ensuite la modlisation de la chane de conversion
photovoltaque puis dcrirons les modles mathmatiques et nous montrerons ensuite
l'influence de la temprature et l'clairement sur le rendement. A la fin on termine par
conclusion.
h : la constante de planck,
C : la vitesse de la lumire.
Ainsi, plus la longueur d'onde est courte, plus l'nergie du photon est grande [4].
UKMO 2012 3
Chapitre 1 Etude du systme photovoltaque
Le soleil met un rayonnement lectromagntique figure (1.1) compris dans une bande de
longueur donde variant de 0,22 10 microns ( m) [3]. Lnergie associe ce rayonnement
solaire se dcompose approximativement ainsi :
Irradiance
(w.m , .m)
UKMO 2012 4
Chapitre 1 Etude du systme photovoltaque
Rs
+
Figure 1.2: Schma lectrique d'une cellule photovoltaque
Les rsistances Rs et Rsh permettent de tenir en compte des pertes lies aux dfauts de
fabrication .Rs reprsente les diverses rsistances de contact et de connexion tandis que Rsh
caractrise les courants de fuite dus diode et aux effets de bord de la jonction [7].
UKMO 2012 5
Chapitre 1 Etude du systme photovoltaque
Photone
Zone dope N
dplacement
d'lectron V
I
Zone dope P
Jonction PN
Silicium monocristallin
Le silicium cristallin est actuellement loption la plus populaire pour les cellules
commerciales, bien que beaucoup dautres matriaux soient disponibles. Le terme cristallin
implique que tous les atomes dans le matriau PV actif font partie dune structure cristalline
simple o il nya aucune perturbation dans les arrangements ordonns des atomes.
UKMO 2012 6
Chapitre 1 Etude du systme photovoltaque
Nouvelle technologie
On utilise de plus en plus de matriaux organiques dans le domaine de loptolectronique,
avec des perspectives dlectronique organique voire molculaire, pour lclairage laide de
diodes lectroluminescentes organiques (OLED : Organic Light- Emitting Diode). Bien que
les optimisations des matriaux mettre en uvre ne soient pas les mmes, le domaine du
photovoltaque bnficie depuis quelques annes des avances technologiques de
loptolectronique. Ainsi, bien que cette filire soit vraiment rcente, les progrs annuels sont
spectaculaires. Les matriaux organiques, molculaires ou polymriques, base de carbone,
dhydrogne et dazote, sont particulirement intressants en termes dabondance, de cot, de
poids et de mise en uvre [11].
Le tableau (1.1) prsente les avantages et les inconvnients pour les technologies les plus
utiliser d'une cellule photovoltaque.
Type Silicium mono Silicium poly Amorphe
cristallin Cristallin
Dure de vie 35 ans 35 ans < 10 ans
Bon rendement en Souplesse Prix moins
Avantage Bon rendement en soleil soleil direct (mois que lev que les
direct le monocristallin mais cristallins
plus que l'amorphe) Bon rendement en diffus
Mauvais rendement en Mauvais rendement en Mauvais rendement en
Inconvnient soleil diffus (temps soleil diffus (temps plein soleil.
nuageux...),prix lev nuageux...),prix lev
Icc
1 cellule
NP cellule cellule
Icc Vco
Ns
0 Vco Vsco
UKMO 2012 8
Chapitre 1 Etude du systme photovoltaque
Avec:
I Np cellule Icc.Np
Np cellule en parallle
Icc Np
1cellule cellule Vco
Icc
0 Vco V
Id
Iph V
Pour un gnrateur PV idel, la tension aux bornes de la rsistance est gale celle aux
bornes de la diode :
"
La diode tant un lment non linaire, sa caractristique I-V est donne par la relation :
*$
#$ #% &'() & , , -
*+
Avec:
#% : Le courant de saturation inverse de la diode.
: la tension au borne de diode
.: =KT/q potentielle thermique
Donc la relation (1.7) sera :
*$
#)/ #% &'() & , , 0
*+
UKMO 2012 10
Chapitre 1 Etude du systme photovoltaque
Iph Rp V
Donc on a :
* 1 23 #
& ,
2)
* 1 # 23
4 &'() & , ,
*+
Avec:
I : Le courant fourni par la cellule
6
5 777
: Le photo-courant dpendant de lclairement (G).
UKMO 2012 11
Chapitre 1 Etude du systme photovoltaque
8
4 >?@ AB (1.15)
9:;<= F
CDE
Pour calculer I on considre Rp =inf donc Ip=0 dans ce cas, l'quation (1.11) devient:
IJK LM
4 <'() <G H3 NOP
F F
UKMO 2012 12
Chapitre 1 Etude du systme photovoltaque
Donc :
UK N IV
23 [?\] ^_ (1.18)
UI KW XYZ
` [a
Telle que:
UK NIV
23 [?^] \ -
UIbcb KW XYZ< F
` [a
Toutes les constantes dans les quations ci-dessus peuvent tre dtermines en utilisant les
donnes de fabricants de panneaux photovoltaques.
La mthode choisie pour la simulation de ce modle est la mthode de Newton Raphson qui
est dcrit comme suit :
e (N
(dJ Sd 0
e (N
Avec : * *f%
* 1 #l 23
#3f #l #% < F
l *+
dJ d 23 * 1 #l 23
#% <l *+F m'() l *+ n
5 o 5 opRg m 1 qr o opRg n
UKMO 2012 13
Chapitre 1 Etude du systme photovoltaque
6
#)/ j 1 sr o opRg < F
777
5 1 sr o opRg
4
54 1 sv o opRg
RS & j . ,
i o Q
w v 0 '() h & 1 , k !
js
Telle que :
Nss : Nombre de module connect en srie (dans notre cas =1)
Npp : Nombre de module connect en parallle (dans notre cas =1)
NS : Nombre de cellule connect en srie.
UKMO 2012 14
Chapitre 1 Etude du systme photovoltaque
Le courant Iph:
[G]
Ki 1000
[Ipv]
[T]
[dT] Ipvn
Tref
Le courant I0:
Ki
Vocn
[dT]
Iscn
[Io]
Kv u
e
[Vta]
[dT]
1
le courant Im:
[V]
Rs
[Npp]
u [Ipv]
[Nss] e
[Npp]
1 [Im]
[I]
[Nss]
q/(n*k*Ns) [Io]
[Vta]
[Npp]
[T]
UKMO 2012 15
Chapitre 1 Etude du systme photovoltaque
Pour la figure (1.8) Nous remarquons que le courant dpend de la temprature puisque le
courant augmente lgrement mesure que la temprature augmente, on constate que la
temprature influe ngativement sur la tension de circuit ouvert. Quand la temprature augmente
la tension de circuit ouvert diminue. Et par contre la puissance maximale du gnrateur subit une
diminution lorsque la temprature augmente figure (1.9).
SP 75 Module photovoltaique
6
G=1000w/m2
4
Courant (A)
T=0C
1 T=25C
T=50C X: 17.35 X: 23.84
T=75C Y: 0.009769 Y: 0.01203
0
0 5 10 15 20X: 19.54 25
Tension (V) Y: 0.0007647
UKMO 2012 16
Chapitre 1 Etude du systme photovoltaque
SP 75 Module photovoltaique
80
G=1000w/m2
70
60
50
Puissance (W) 40
30
20
T=0C
T=25C
10 T=50C
T=75C
0
0 5 10 15 20 25
Tension (V)
UKMO 2012 17
Chapitre 1 Etude du systme photovoltaque
SP 75 module photovoltaique
6
X: 0 T=25C
Y: 4.8
5
X: 0
Y: 3.84
4
Courant (A)
X: 0
Y: 2.88
3
X: 0
Y: 1.92
2 G=200w/m2
X: 0 G=400w/m2
Y: 0.96
G=600w/m2
1
G=800w/m2
G=1000w/m2
0
0 5 10 15 20 25
Tension (V)
SP 75 Module photovoltaique
80
70 T=25C
60
50
Puissance (W)
40
30
G=200w/m2
20 G=400w/m2
G=600w/m2
10 G=800w/m2
G=1000w/m2
0
0 5 10 15 20 25
Tension (V)
UKMO 2012 18
Chapitre 1 Etude du systme photovoltaque
1.10 Conclusion
Dans ce chapitre nous avons prsent les principales caractristiques et les technologiques
des lments constitutifs dun gnrateur PV et nous avons montr comment augmenter le
courant ou la tension d'un gnrateur photovoltaque ainsi nous avons montr bien
l'influence de la temprature et l'clairement sur le rendement de la cellule, et on constate que
la puissance ne dponde pas seulement de la temprature mais dponde aussi de lclairement.
UKMO 2012 19
Chapitre 2 Les convertisseurs statiques
2.1 Introduction
Les convertisseurs sont des appareils servent transformer la tension continue fournie par
les panneaux ou les batteries pour l'adapter des rcepteurs fonctionnant en une tension
continue diffrente ou une tension alternative.
Ltude du convertisseur est intressante dans la mesure o il est utilis dans la plupart des
nouveaux types de sources de production dnergie disperse connecte au rseau (olienne,
photovoltaque, pile combustible).
Onduleur
Redresseur
Dans ce chapitre nous prsentons les dfrents types des convertisseurs statiques quon
peut utiliser dans le systme photovoltaque. Premirement on commence par les
convertisseurs (DC-AC) et sont dfrents type, et en termine par les convertisseurs (DC-DC)
[15].
2.2 Les onduleurs
tension continue, en tension alternative contrle de faon trs prcise. La source de tension
continue soit un arognrateur ou des panneaux solaires.
La commande de l'onduleur est base sur la prdiction de la tension de sortie dun pas en
avant que nous appellerons Dead beat control par cette commande, la tension de sortie de
londuleur est force de suivre une rfrence sinusodale chantillonne pour la production
la sortie du filtre une onde proche dune sinusode avec un taux de distorsion harmonique trs
rduit.
UKMO 2012 20
Chapitre 2 Les convertisseurs statiques
Onduleur Triphas.
Le schma de principe dun tel onduleur mont en demi-pont est montr sur la figure (2.2).
E/2 D1 S1
.
Charge
0 Vo(t)
.
S2
E/2 D2
.
Lors de la fermeture de linterrupteur S1, la tension aux bornes de la charge serait donc de
+ E/2, et prend la valeur E/2 quand le second interrupteur, S2 est ferm
UKMO 2012 21
Chapitre 2 Les convertisseurs statiques
La conduction simultane des deux interrupteurs est vite par llaboration dune
commande adquate qui tient compte des diffrentes caractristiques des imperfections de ces
interrupteurs de puissance (temps de monte tr, temps de descente tf et temps de stockage ts).
Les diodes D1. D2, dites de rcupration, assurent la conduction dun courant ngatif en
cas de dphasage de ce dernier par rapport la tension aux bornes de la charge.
Deux interrupteurs du mme bras ne peuvent pas conduire simultanment, cause dun
court-circuit de la source de tension continue.
.
.
D1 S1 D3 S3
.
Charge .
E
Vo(t)
.
S4 S2
D4 D2
.
UKMO 2012 22
Chapitre 2 Les convertisseurs statiques
.
E S1 S5 S3
.
u v w
.
S4 S2 S6
E
.
Figure 2.4: Schma de Principe dun Onduleur Triphas En Pont
Figure (2.4) illustre la topologie dun onduleur triphas six interrupteurs de puissance. Le
dcalage entre les signaux de commande est de 60 [16].
UKMO 2012 23
Chapitre 2 Les convertisseurs statiques
L K1
K2 C R
Ve Vt Vs
(2.1)
UKMO 2012 24
Chapitre 2 Les convertisseurs statiques
Fonctionnement
L'interrupteur K1 est ferm pendant la fraction T de la priode de dcoupage T. La
source d'entre fournit l'nergie la charge R au travers de l'inductance L.
Lors du blocage du transistor, la diode K2 assure la continuit du courant dans l'inductance.
L'nergie emmagasine dans cette inductance est alors dcharge dans le condensateur et la
rsistance de la charge.
Les formes d'ondes en conduction continue sont reprsentes la figure (2.6). En rgime permanent, la
valeur moyenne de la tension aux bornes de l'inductance est nulle, ce qui impose la relation suivante :
(2.2)
(2.3)
K2 L
Ve C R Vs
K1
Le hacheur srie est un convertisseur direct DCDC, La source d'entre est de type
tension continue et la charge de sortie continue de type source de courant. L'interrupteur K1
peut tre remplac par un transistor puisque le courant est toujours positif et que les
commutations doivent tre commandes (au blocage et l'amorage).
Fonctionnement
UKMO 2012 25
Chapitre 2 Les convertisseurs statiques
(2.4)
Malgr le rendement lev du convertisseur buck dans les systmes avec des sources de
puissance conventionnelles, le convertisseur boost peut tre plus appropri aux systmes
photovoltaques avec le suiveur du point de puissance maximale (MPPT) puisque le
convertisseur fonctionne au mode de courant continu extrayant autant de puissance que
possible partir des cellules solaires. Par consquent le rendement nergtique du
convertisseur boost peut tre plus grand que le convertisseur buck. Le convertisseur boost est
gnralement employ pour obtenir une tension plus leve de sortie, tandis que le
convertisseur buck est employ pour abaisser la tension de sortie [18].
UKMO 2012 26
Chapitre 2 Les convertisseurs statiques
Scope2
Pulse
Generator
g m
i i
+ a k +
- -
Current Measurement1 Gto1 Current Measurement
i
-
Current Measurement2 Series RLC Branch +
v In Mean
-
+
DC Voltage Source
Voltage Measurement Scope
Mean Value
Diode1
In Mean
Scope3
Mean Value1
a) Charge rsistive R
Le tableau ci-dessous reprsente la variation de la valeur moyenne de la tension et du courant
en fonction de rapport cyclique :
UKMO 2012 27
Chapitre 2 Les convertisseurs statiques
Uc
25
20
15
tension(v)
10
0
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
temp(t)
Ucmoy
14
12
10
courant(A)
0
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
temp(t)
UKMO 2012 28
Chapitre 2 Les convertisseurs statiques
Icmoy
0.45
0.4
0.35
0.3
courant(A )
0.25
0.2
0.15
0.1
0.05
0
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
temp(t)
UKMO 2012 29
Chapitre 2 Les convertisseurs statiques
La tension Uc
25
20
15
T ension(v)
10
-5
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(s)
La tension Ucmoy
15
10
Tension(V)
0
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(s)
UKMO 2012 30
Chapitre 2 Les convertisseurs statiques
Le courant Icmoy
0.5
0.45
0.4
0.35
0.3
ourant(A)
0.25
C
0.2
0.15
0.1
0.05
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Temps(s)
La tension de sortie du hacheur nest pas continue mais toujours positive. Lorsque la
priode est assez faible (frquence de 100 1000 Hz) la charge ne voit pas les crneaux
mais la valeur moyenne de la tension.
On remarque quelle que soit la nature de la charge, on aura Le hacheur
srie est bien abaisseur de tension (hacheur dvolteur).
i
+ i
- +
-
ipv ipv1
L
R
+
v +
- v
Rs -
vpv C
C1
DC Voltage Source
In Mean
I
Pulse Mean Value
Generator
In Mean
V
Mean Value1
UKMO 2012 31
Chapitre 2 Les convertisseurs statiques
Uc
25
24.5
Tonsion(V)
24
23.5
23
0 0.002 0.004 0.006 0.008 0.01
Temp(T)
Ucmoy
70
60
50
T o n s io n (V )
40
30
20
10
0
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temp(T)
UKMO 2012 32
Chapitre 2 Les convertisseurs statiques
Icmoy
1.4
1.2
c o u ra n t(A )
0.8
0.6
0.4
0.2
0
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temp(T)
On remarque que la tension de sortie ne pas le mme que l'entre ce qui exprime le rle
de convertisseur lvateur (boost).
2.6 Conclusion
Dans ce chapitre nous avons prsent les caractristiques lectriques des convertisseurs
statiques (convertisseur DC-DC/DC-AC).nous avons montr par simulation que la valeur
moyenne de la tension de sortie peut tre ajuste en jouant sur la valeur du rapport cyclique.
Nous avons bas sur ltude des convertisseurs du courant continu (boost-buck) pour le but
d'utiliser dans le chapitre suivante.
UKMO 2012 33
Chapitre 3 Les mthodes d'optimisation
3.1 Introduction
Le cot lev des gnrateurs photovoltaques peut justifier ladjonction dun
appareillage annexe, mme sophistiqu, permettant de grer au mieux lnergie disponible. En
particulier, sur la caractristique courant-tension du gnrateur (qui dpend des conditions
dclairement, de temprature, de vieillissement) il existe un point de fonctionnement o la
puissance dbite est maximale.
Puissance charge
UKMO 2012 34
Chapitre 3 Les mthodes d'optimisation
Figure 3.2 : Points de fonctionnement d'un GPV en connexion directe, en fonction de la charge.
Comme le montre dans la figure (3.2), un GPV peut tre connect directement trois
types de charges :
- une charge de type source de tension continue,
- une charge de type source de courant continue,
- une charge purement rsistive.
UKMO 2012 35
Chapitre 3 Les mthodes d'optimisation
Adaptateure
I1 I2
Convertisseur V
Charge
Ve statique Vs DC
GPV
GPV
Figure 3.3 : Etage d'adaptation d'un GPV-charge
=1-
UKMO 2012 36
Chapitre 3 Les mthodes d'optimisation
Puisque le rapport est infrieur 1 ( < 1), le convertisseur ne joue le rle dun lvateur
que si la charge RS remplit la condition suivante:
Dans les conditions optimales et pour une charge RS donne, la rsistance interne du panneau
( Rpv = Ropt ) et le rapport cyclique ( = opt ) obissent donc lquation:
50 2
600W/m
Courant (A)
40 600W/m2
2 3
400W/m
30
400W/m2
2
20 200W/m2
200W/m2
10 1
0
0 5 10 15 20 25 0
Tension (V) 0 5 10 15 20 25
Tension (V)
La Figure (3.5) prsente le schma de principe dun module photovoltaque dot dun
tage dadaptation DC-DC entre le GPV et la charge de sortie. Cette structure correspond un
systme plus communment appel systme autonome. Il permet le plus souvent dalimenter
une batterie servant de stocke lnergie pour une charge qui ne supporte pas les fluctuations
de tension. Cet tage dadaptation dispose dune commande MPPT (Maximum Power Point
UKMO 2012 37
Chapitre 3 Les mthodes d'optimisation
Tracking) qui lui permet de rechercher le PPM que peut fournir un panneau solaire
photovoltaque. Lalgorithme de recherche MPPT peut tre plus ou moins complexe en
fonction du type dimplantation choisi et des performances recherches. Cependant au final,
tous les algorithmes performants doivent jouer sur la variation du rapport cyclique du
convertisseur de puissance associ [4].
I1 I2
V1 V2
Rapport
GPV
pour assurer le fonctionnement d'un GPV a son point de puissance maximale (PPM),des
contrleur MPPT sont utiliss ,ces contrleurs sont destins a minimiser l'erreur entre la
Puissance de fonctionnement et puissance maximale de rfrence variable en fonction des
conditions climatiques baser sur des mthode optimale (voir la suite)[21].
La figure (3.6) illustre trois cas de perturbations. Suivant le type de perturbation, le point
de fonctionnement bascule du point de puissance maximal PPM1 vers un nouveau point P1 de
fonctionnement plus ou moins loign de loptimum.
Dans ce cas est pour une variation densoleillement il suffit de rajuster la valeur du
rapport cyclique pour converger vers le nouveau point de puissance maximum PPM2. Dans le
cas b pour une variation de charge on peut galement constater une modification du point de
fonctionnement qui peut retrouver une nouvelle position optimale grce laction dune
commande. Enfin dans le dernier cas c de variation de point de fonctionnement peut se
UKMO 2012 38
Chapitre 3 Les mthodes d'optimisation
la variation la variation
de ! PPM2 de P1
variation
P1 variation de de la
l'clairement charge
E2 E1
VPV(v) VPV(v)
(a) (b)
Charge: constante
T: variable
L'clairement: constante
PPV(w) Caractristique de
la variation PPM1 la charge
de ! PPM2
P1
La variation de
E1 E1 la Temprature
T2 T1
VPV(v)
(c)
Figure 3.6 : Recherche et recouvrement du Point de Puissance Maximale a) suite une variation
d'clairement, b) suite une variation de charge, c) suite une variation de temprature
UKMO 2012 39
Chapitre 3 Les mthodes d'optimisation
"#$$
K1= "%&
'() * (3.8)
UKMO 2012 40
Chapitre 3 Les mthodes d'optimisation
+,--
K2 = / '() * (3.9)
+ .
Cette mthode prsente les mmes avantages et inconvnients que ceux de la mthode de la
tension de circuit ouvert.
UKMO 2012 41
Chapitre 3 Les mthodes d'optimisation
Ppv
PPM
p<0
Puissance(w) p>0
v>0 v<0
Vpv
La figure (3.8) reprsente lalgorithme classique associ une commande MPPT de type
P&O, o lvolution de la puissance est analyse aprs chaque perturbation de tension. Pour
ce type de commande, deux capteurs (courant et tension du GPV) sont ncessaires pour
dterminer la puissance du PV chaque instant [24].
UKMO 2012 42
Chapitre 3 Les mthodes d'optimisation
Oui
!"
0 0 0 0
Telle que :
03 789: 789:56
UKMO 2012 43
Chapitre 3 Les mthodes d'optimisation
donc trouver un compromis entre la prcision et la rapidit. Ce qui rend cette commande
difficile optimiser [24].
789 89 89
;><=
0@ ;"<=
UKMO 2012 44
Chapitre 3 Les mthodes d'optimisation
;$<=
Ppv =0
"<=
;$<= ;$<=
>0 PPM <0
"<= "<=
Puissance(w)
Vpv
Figure 3.9 : Trajectoire par Incrmentation de Conductance
0
0
Non Oui
0 "
Oui
0><= 5><= 0A3 "
0"<= "<=
Non
Non Oui
Oui Oui
0><= 5><=
! 0A3 !"
0"<= "<=
0 0 0 0
# $
;$<=
B : Le point de fonctionnement est gauche du PPM.
;"<=
;$<=
;"<=
0 : Le point de fonctionnement sur le PPM.
Charge
Vpv
+
% &
- +
Vre
Figure 3.11 : Mthode contre raction de la tension
Les inconvnients de cette configuration sont les mmes que pour la mthode de
connexion directe (gnrateur PV + charge de profil). C'est--dire, le systme ne peut pas
s'adapter aux conditions variables d'environnement, telles que l'irradiante et la temprature.
Les avantages de cette technique c'est une technique simple, conomique et emploie
seulement une commande de boucle de feedback [28].
UKMO 2012 46
Chapitre 3 Les mthodes d'optimisation
Dans toutes les mthodes (MPPT) prsentes jusqu'ici, une mesure de Ipv et Vpv a t
employe pour obtenir des information sur la puissance de sortie du panneau. Cette nouvelle
mthode (MPPT) utiliser seulement une mesure de courant pour obtenir l'information sur le
point oprationnel et actuel du systme figure (3.13).
L'avantage de cette mthode est qu'elle est simple, contrairement aux mthodes mentionnes
prcdemment [12].
I1 I2
Capteur du
courant
Contrleur
MPP
3.7 Conclusion
A travers de ce chapitre, nous avons montr les limites de la connexion directe et la
caractristique I-V d'un panneau photovoltaque n'est pas linaire.
Nous avons vu galement les algorithmes MPPT ont t couverts en dtail. Comme ont t
faits connatre, il y a beaucoup de manires distinguer et de grouper les techniques pour
suivre le MPP au gnrateur PV. Cependant, en ce chapitre, les techniques directes et
indirectes taient choisies et dveloppes en dtail.
UKMO 2012 47
Chapitre 3 Les mthodes d'optimisation
D'autre part, la technique directe offre les avantages suivants: ni une grande base de donnes,
ni une grande mmoire ncessaire pour calculer le MPP; elle est totalement souple en ce qui
concerne le profil de charge, le suiveur de MPP est indpendant de la variation des paramtres
du gnrateur PV.
UKMO 2012 48
Chapitre 4 Simulation du systme photovoltaque
4.1 Introduction
Figure 4.1: Schma synoptique du systme PV contrl par une commande MPPT
Pour raliser cette simulation, nous avons utilis le logiciel MATLAB comme outil de
tests et de simulation.
UKMO 2012 49
Chapitre 4 Simulation du systme photovoltaque
UKMO 2012 50
Chapitre 4 Simulation du systme photovoltaque
-1
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(s)
23
T ens ion(V )
22
21
20
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(s)
11
P uis s anc e(w )
10
8
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps
Figure 4.3 : La puissance, la tension et le courant produit par le GPV sans MPPT
UKMO 2012 51
Chapitre 4 Simulation du systme photovoltaque
Figure 4.4: Schma bloc de la simulation d'un systme PV avec la commande MPPT
UKMO 2012 52
Chapitre 4 Simulation du systme photovoltaque
UKMO 2012 53
Chapitre 4 Simulation du systme photovoltaque
La charge 50 30 10 8 3.86
L'clairement 1000w/m^2 1000w/m^2 1000w/m^2 1000w/m^2 1000w/m^2
La puissance 74.86w 74.86w 74.86w 74.86w 74.81w
Le courant 4.35A 4.35A 4.37A 4.36A 4.40A
La tension 17.17V 17.14V 17.12V 17.16V 16.99V
UKMO 2012 54
Chapitre 4 Simulation du systme photovoltaque
Courant de la charge
1.2
Courant(A) 0.8
0.6
0.4
0.2
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(s)
Courant de GPV
4.45
4.4
Courant(A)
4.35
4.3
4.25
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(s)
50
Tension(V)
40
30
20
10
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(s)
Tension de GPV
17.5
Tension(V)
17
16.5
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(s)
UKMO 2012 55
Chapitre 4 Simulation du systme photovoltaque
La puissance de la charge
80
60
Puissance (W) 40
20
0
0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(s)
La puissanse de GPV
74.9
Puissance(W)
74.8
74.7
74.6
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(s)
UKMO 2012 56
Chapitre 4 Simulation du systme photovoltaque
10.5
P sans MPPT
10
8.5
8
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(s)
74.9
P uis s ance(w)
74.7
74.6
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(s)
Pour montrer bien l'efficacit de la commande MPPT sur la puissance maximale, nous avons
simul le systme avec diffrents clairements et une temprature fixe T=25C, les rsultats
obtenus est la suivante :
UKMO 2012 57
Chapitre 4 Simulation du systme photovoltaque
Sans MPPT
La tension de la charge pour diffrents clairement
24
G=800w/m2
22
20
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
22
20 G= 600w/m2
18
Tension(v)
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
22
20
G= 400w/m2
18
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
19.345
G= 200w/m2
19.345
19.345
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(s)
V sortie de convertisseur
100
G= 800w/m2
50
0
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
50
G= 600w/m2
Tension (v)
0
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
40
G= 400w/m2
20
0
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
25
G= 200w/m2
20
15
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(w)
UKMO 2012 58
Chapitre 4 Simulation du systme photovoltaque
Avec MPPT
La puissance max pour diffrents clairement
10
8 P pour 800w/m2
6
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
10
P pour 600w/m2
8
Puissance(w)
6
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
10
P pour 400w/m2
8
6
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
7.4845
P pour 200w/m2
7.4845
7.4845
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(s)
59.8
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
44.6
44.2
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
28.8
MPPT pour 400w/m2
28.6
28.4
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
12.426
MPPT pour 200w/m2
12.424
12.422
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Temps(s)
UKMO 2012 59
Chapitre 4 Simulation du systme photovoltaque
A travers les figures (4.11, 4.12 et 4.13, 4.14) obtenus nous avons montr que le
convertisseur DC-DC et la commande MPPT effectuent correctement leurs rles. Le
convertisseur fournit dans les conditions optimales une tension sa sortie suprieure celle
fournie par le gnrateur PV. La commande MPPT adapte le gnrateur PV la charge:
transfert de la puissance maximale. En remarque toujours que le rsultat reste le mme pour
diffrents clairements.
4.4 Conclusion
Dans ce chapitre nous avons prsent la simulation d'un systme PV avec et sans MPPT
et on conclue que un systme de conversion d'nergie, constitu par un panneau PV, et un
convertisseur DC-DC et un suiveur du point de puissance maximale (MPPT) par la mthode
de Perturbe et observe conduit un fonctionnement optimal du systme. La recherche du
point de puissance maximale est fait d'une manire stable et robuste, les rsultats obtenus
montrent que pendant une dure infrieure 2 ms et en variant dune manire automatique le
rapport cyclique du signal qui commande linterrupteur du convertisseur dnergie le systme
converge vers les conditions optimales indpendamment des conditions mtorologiques et de
la variation de la charge ces rsultats, forts intressants, montrent que lutilisation de la
commande MPPT permet damliorer dune manire considrable le rendement des
installations photovoltaques.
UKMO 2012 60
Conclusion gnrale
Conclusion gnrale
La demande mondiale en nergie volue rapidement et les ressources naturelles de
l'nergie telles que l'uranium, le gaz et le ptrole diminuent en raison d'une grande diffusion et
dveloppement de l'industrie ces dernires annes. Pour couvrir les besoins en nergie, des
recherches sont conduits l'nergie renouvelable. Une des nergies renouvelables qui peut
accomplir la demande est lnergie solaire photovoltaque, cest une nergie propre,
silencieuse, disponible et gratuite. Cest dailleurs ce explique que son utilisation connat une
croissance significative dans le monde.
L'utilisation d'nergie solaire comme source alternative d'nergie, souffre du cot lev
des cellules solaires, du faible rendement et de puissance intermittent selon la fluctuation des
conditions atmosphriques. Par consquent, n'importe quelle conception de systme
d'application d'nergie solaire, devrait prendre en compte ces inconvnients. Ce travail sagit
dune contribution l'optimisation du fonctionnement dun systme photovoltaque.
Les travaux prsents concernent plus particulirement la problmatique du couplage entre un
gnrateur photovoltaque et une charge de type continue. En effet, ce type de couplage
souffre encore du problme du transfert de puissance du gnrateur photovoltaque la
charge qui est souvent trs loin des possibilits du gnrateur.
Avant daborder la problmatique de conversion, nous avons tout dabord tudi le
contexte et la problmatique de lnergie solaire photovoltaque. Ceci savrer ncessaire pour
acqurir une connaissance approfondie du comportement dun gnrateur photovoltaque.
Nous avons tabli des programmes de calculs sous lenvironnement Matlab 7.0 permettant de
tracer la caractristique I-V pour diffrentes clairements et tempratures solaires et de
calculer la puissance maximale dlivre par la cellule photovoltaque correspondante.
Nous rappelons dans le deuxime chapitre les principales caractristiques des convertisseurs
statiques dont il faut tenir compte pour dimensionner les tages dadaptation entre le
gnrateur photovoltaque et la charge.
Pour mieux analyser les contraintes de lexploitation dun systme photovoltaque, nous
avons tudie plusieurs mthodes doptimisation.
Dans le dernier chapitre, nous avons analys dans le simulateur Matlab 7.0 le fonctionnement
optimal dun systme PV, de moyennes puissances dont la rgulation de la puissance est
effectue par une commande MPPT, suites aux variations des conditions mtorologiques
(clairements,) et de la charge.
Conclusion gnrale
Rfrences bibliographiques
[1] K. Kassmi et M. Hamdaoui et F. Olivi Conception et modlisation dun
systme photovoltaque adapt par une commande MPPT analogique, universit de
maroc ,revue des nergies renouvelables ,pp 451 462, 2007.
[2] A. T.SINGO Systme dalimentation photovoltaque avec stockage hybride
pour lhabitat nergtiquement autonome Thse de doctorats, universit Henri
Poincar, Nancy-I,2010.
[3] J. Royer et T. Djiako et E. Schiler, B. Sadasy. Le pompage photovoltaque,
universit d'Ottawa, 1998.
[4] S.ABADA Etude et optimisation d'un gnrateur photovoltaque pour la recharge
d'une batterie avec un convertisseur sepic, Mmoire Matre es Sciences , Universit
laval, 2011.
[5] M. G. Villalva and J. R. Gazoli, and E. R. Filho Comprehensive approach to
modeling and simulation of photovoltaic arrays,Jornal of power electronics, ,May
2009.
[6] A. Bilbao Ralisation d'un commande mppt numrique universit de rovira i
virgili ,2006.
[7] H. BELGHITRI et H. Modlisation, simulation et optimisation dun systme
hybride olien-photovoltaque mmoire de magister, 2009.
[8] A. labouret et M. villoz prface de jean louis bal Energie solaire
photovoltaque, livre 4 dition,2008.
[9] B. Flche - D. Delagnes Energie solaire photovoltaique.doc, juin 2007,
http://www.scribd.com.
[10] M.belhad Modlisation Dun Systme De Captage Photovoltaque Autonome,
mmoire de magister universit de bachar ,2008.
[11] S. PETIBON Nouvelles architectures distribues de gestion et de conversion
de lnergie pour les applications photovoltaques, universit de Toulouse, 2009.
[12] L. ABBASSEN Etude de la connexion au rseau lectrique d'une centrale
photovoltaque Mmoire de magister, universit mouloud Mammeri Tizi ouzzo,
2011.
[13] L. ZAROUR Etude technique dun systme dnergie hybride photovoltaque-
olien hors rseau, Thse Magister universit de constantine, 2007.
[14] Z. Bendjellouli Contribution a la modlisation d'une cellule solaire, Mmoire
de magister universit de bachar, 2009.
Rfrences bibliographiques
Annexe 1:
Programme associ avec le panneau solaire :
clc;
%% paneaux solar sp 75 module
Iscn =4.8; %Nominal short-circuit courant [A]
Vocn =21.7; %Nominal array open-circuit voltage[v]
Imp = 4.41; %Array current @ maximum power point [A]
Vmp = 17; %Array Voltage @ maximum power point [v]
Pmax_e = Vmp*Imp; %Array maximum power [W]
Kv = -0.77; %Voltage/temperature coefficient [V/K]
Ki = 2.06e-3; %Current/temperature coefficient [A/K]
Ns = 32; %Nunber of series cells
Tn = 25+273.15; %module temperateur[k]
Gn=1000; %Irradiance
%% parameters module
Npp=1;
Nss=1;
%% converter Boost
L=0.34*10^-3;
Ce=9400*10^-6;
Co=4700*10^-6;
%Vo=12;
%% Constants
%% Algorithm parameters
%Increment of Rs
Rsinc = 0.01;
%Maximum tolerable power error
tol = 0.0001;
%Maximum number of iterations for each value of "a"
nimax = 5000;
%Voltage points in each iteration
nv = 50;
%% Adjusting algorithm
% Reference values of Rs and Rp
Rs_max = (Vocn - Vmp)/ Imp;
Rp_min = Vmp/(Iscn-Imp) - Rs_max;
% Initial guesses of Rp and Rs
Rs = 0;
Rp = Rp_min;
% The model is adjusted at the nominal condition
T = Tn;
G = Gn;
Vtn = k * Tn / q; %Thermal junction voltage (nominal)
Vt = k * T / q; %Thermal junction voltage (current temperature)
Annexe
ni = 0;
ni = ni + 1;
Io = (Ipv - Vocn/Rp)/(exp(Vocn/Vt/a/Ns)-1); %%
% Increments Rs
Rs = Rs + Rsinc;
Rp_ = Rp;
Rp = Vmp*(Vmp+Imp*Rs)/(Vmp*Ipv-
Vmp*Io*exp((Vmp+Imp*Rs)/Vt/Ns/a)+Vmp*Io-Pmax_e);
% Solving the I-V equation for several (V,I) pairs
clear V
clear I
V = 0:Vocn/nv:Vocn; % Voltage vector
I = zeros(1,size(V,2)); % Current vector
g(j) = Ipv-Io*(exp((V(j)+I(j)*Rs)/Vt/Ns/a)-1)-(V(j)+I(j)*Rs)/Rp-I(j);
g(j) = Ipv-Io*(exp((V(j)+I(j)*Rs)/Vt/Ns/a)-1)-(V(j)+I(j)*Rs)/Rp-I(j);
glin(j) = -Io*Rs/Vt/Ns/a*exp((V(j)+I(j)*Rs)/Vt/Ns/a)-Rs/Rp-1;
I_(j) = I(j) - g(j)/glin(j);
I(j) = I_(j);
end
Pmax_m = max(P);
Annexe
perror = (Pmax_m-Pmax_e);
end % while (error>tol)
if (Rp<0)
Rp = Rp_
end
Ion = Io;
%% Outputs
% I-V curve
figure
subplot(2,1,1)
grid on
hold on
title('Adjusted I-V curve');
xlabel('V [V]');
ylabel('I [A]');
xlim([0 Vocn*1.1]);
ylim([0 Iscn*1.1]);
plot(V,I,'LineWidth',2,'Color','b')
plot([0 Vmp Vocn ],[Iscn Imp 0 ],'o','LineWidth',2,'MarkerSize',5,'Color','r')
%%P-V curve
subplot(2,1,2)
grid on
hold on
title('Adjusted P-V curve');
xlabel('V [V]');
ylabel('P [W]');
xlim([0 Vocn*1.1]);
ylim([0 Vmp*Imp*1.1]);
plot(V,P,'LineWidth',2,'Color','b')
disp(sprintf('Model info:\n'));
disp(sprintf(' Rp_min = %f',Rp_min));
disp(sprintf(' Rp = %f',Rp));
disp(sprintf('Rs_max = %f',Rs_max));
disp(sprintf(' Rs = %f',Rs));
disp(sprintf(' a = %f',a));
disp(sprintf(' T = %f',T-273.15));
disp(sprintf(' G = %f',G));
disp(sprintf(' Pmax,m = %f (model)',Pmax_m));
disp(sprintf(' Pmax,e = %f (experimental)',Pmax_e));
disp(sprintf(' tol = %f',tol));
disp(sprintf('P_error = %f',perror));
disp(sprintf(' Ipv = %f',Ipv));
disp(sprintf(' Isc = %f',Isc));
disp(sprintf(' Ion = %g',Ion));
disp(sprintf('\n\n'));