Correction d’exercice sur le load flow
On considère le réseau suivant (7 bus) :
On cherche à déterminer l’état de celui-ci quant à la répartition de puissance.
Données des bus :
Bus Type Grandeurs connues Grandeurs
inconnues
1 Générateur contrôlé V1 = 1,03 pu – P1 = 2.2 pu Q1 – 1
2 Générateur contrôlé V2 = 1,01 pi – P2 = 1.2 pu Q2 –
3 Bus balancier V3 = 1,05 pu - = 0° P3 – Q3
4 Bus charge P4 = 0 – Q4 = 0 V4 -
5 Bus charge P4 = 0 – Q4 = 0 V5 -
6 Bus charge P4 = 0 – Q4 = 0 V6 -
7 Bus charge P7 = 3,8 pu – Q7 = 2,2 pu V7 -
Données des lignes :
Ligne Admittance
4-5 -j10
4-6 -j6
5-7 -j8
6-7 -j12
Données des transformateurs :
Transformateur Admittance
1 12,31-8,46j
2 8,8-104,52j
3 5,73-66,17j
Matrice d’admittances
(i,j) 1 2 3 4 5 6 7
1 12,31-98,46j 0 0 -12,31+98,46j 0 0 0
2 0 8,8-104,52j 0 0 -8,8+104,52j 0 0
3 0 0 5,73-66,17j 0 0 -5,73+66,17j 0
4 -12,31+98,46j 0 0 12,31-114,46j 10j 6j 0
5 0 -8,8+104,52j 0 10j 8,8-122,52j 0 8j
6 0 0 -5,73+66,17j 6j 0 5,73-84,17j 12j
7 0 0 0 0 8j 12j -
20j
Résultats du calcul
Bus Type Grandeurs connues Grandeurs calculées
1 Générateur contrôlé V1 = 1,03 pu – P1 = 2.2 pu Q1 = 0.5287 – 1= 13.49°
2 Générateur contrôlé V2 = 1,01 pi – P2 = 0.2 pu Q2 = 1.2804 – =8.1659 °
3 Bus balancier V3 = 1,05 pu - = 0° P3=0.4151 – Q3=2.3345
4 Bus charge P4 = 0 – Q4 = 0 V4=1.0224 - =12.329°
5 Bus charge P5 = 0 – Q5 = 0 V5= 0.9971 - =7.5756°
6 Bus charge P6 = 0 – Q6 = 0 V6= 1.0161 - =-0.1714°
7 Bus charge P7 = 3,8 pu – Q7 = 2,2 pu V7= 0.8522 - =-9.909°
Important :
1- Les bus 4, 5 et 6 présentent des puissances nulles puisqu’ils ne sont pas des bus de
charge (Load bus).
Les calculs ont généré les valeurs suivantes :
S4 = -9.4476e-06 + 1.1630e-06i
S5 = -9.5483e-06 + 8.2839e-07i
S6 = -4.1026e-07 - 2.1129e-08i
2- Le nombre d’itérations pour atteindre une stabilité de toutes les grandeurs se situe après
297 (le choix a été fait pour 300)
Programme :
Remarque : J’ai utilisé le nombre complexe « i » au lieu de « j »
%Matrice d'admittances
Y= [(12.308 - 98.462i) 0 0 (-12.308 + 98.462i) 0 0 0;
0 (8.8-104.52i) 0 0 (-8.8+104.52i) 0 0;
0 0 (5.7347 - 66.1697i) 0 0 (-5.7347 + 66.1697i) 0;
(-12.308 + 98.462i) 0 0 (12.308 - 114.462i) 10i 6i 0;
0 (-8.8018 + 104.5219i) 0 10i (8.8018 - 122.5219i) 0 8i;
0 0 (-5.7347 + 66.1697i) 6i 0 (5.7347 - 84.1697i) 12i;
0 0 0 0 8i 12i -20i];
%Données constantes
V1 = 1.03;
V1cor = 1.03;
P1 = 2.2;
V2 = 1.01;
V2cor = 1.01;
P2 = 1.2;
V3 = 1.05;
S4 = 0;
S5 = 0;
S6 = 0;
% S7 est comptée négativement puisque S7=Sg-Sd=0-(3.8+2.2i)= -3.8-2.2i
S7 = -3.8-2.2i;
%Initialisation
V4 = 1;
V5 = 1;
V6 = 1;
V7 = 1;
Q1 = 0;
Q2 = 0;
%Calcul
for m = 1:300
disp(m)
V4 = 1/Y(4,4)*(conj(S4/V4)-(Y(4,1)*V1cor+Y(4,2)*V2cor+Y(4,3)*V3+Y(4,5)*V5+Y(4,6)*V6+Y(4,7)*V7))
V5 = 1/Y(5,5)*(conj(S5/V5)-(Y(5,1)*V1cor+Y(5,2)*V2cor+Y(5,3)*V3+Y(5,4)*V4+Y(5,6)*V6+Y(5,7)*V7))
V6 = 1/Y(6,6)*(conj(S6/V6)-(Y(6,1)*V1cor+Y(6,2)*V2cor+Y(6,3)*V3+Y(6,4)*V4+Y(6,5)*V5+Y(6,7)*V7))
V7 = 1/Y(7,7)*(conj(S7/V7)-(Y(7,1)*V1cor+Y(7,2)*V2cor+Y(7,3)*V3+Y(7,4)*V4+Y(7,5)*V5+Y(7,6)*V6))
Q1 = imag(V1cor*conj(Y(1,1)*V1cor+Y(1,2)*V2cor+Y(1,3)*V3+Y(1,4)*V4+Y(1,5)*V5+Y(1,6)*V6+Y(1,7)*V7))
S1 = P1 + Q1*i
V1 = 1/Y(1,1)*(conj(S1/V1cor)-(Y(1,2)*V2cor+Y(1,3)*V3+Y(1,4)*V4+Y(1,5)*V5+Y(1,6)*V6+Y(1,7)*V7))
V1cor = 1.03*(cos(angle(V1))+i*sin(angle(V1)))
Q2 = imag(V2*conj(Y(2,1)*V1cor+Y(2,2)*V2cor+Y(2,3)*V3+Y(2,4)*V4+Y(2,5)*V5+Y(2,6)*V6+Y(2,7)*V7))
S2 = P2 + Q2*i
V2 = 1/Y(2,2)*(conj(S2/V2cor)-(Y(2,1)*V1cor+Y(2,3)*V3+Y(2,4)*V4+Y(2,5)*V5+Y(2,6)*V6+Y(2,7)*V7))
V2cor = 1.01*(cos(angle(V2))+i*sin(angle(V2)))
end
S3=V3*conj(Y(3,1)*V1cor+Y(3,2)*V2cor+Y(3,3)*V3+Y(3,4)*V4+Y(3,5)*V5+Y(3,6)*V6+Y(3,7)*V7)
S4=V4*conj(Y(4,1)*V1cor+Y(4,2)*V2cor+Y(4,3)*V3+Y(4,4)*V4+Y(4,5)*V5+Y(4,6)*V6+Y(4,7)*V7)
S5=V5*conj(Y(5,1)*V1cor+Y(5,2)*V2cor+Y(5,3)*V3+Y(5,4)*V4+Y(5,5)*V5+Y(5,6)*V6+Y(5,7)*V7)
S6=V6*conj(Y(6,1)*V1cor+Y(6,2)*V2cor+Y(6,3)*V3+Y(6,4)*V4+Y(6,5)*V5+Y(6,6)*V6+Y(6,7)*V7)
S7=V7*conj(Y(7,1)*V1cor+Y(7,2)*V2cor+Y(7,3)*V3+Y(7,4)*V4+Y(7,5)*V5+Y(7,6)*V6+Y(7,7)*V7)
%Ligne 1-4, Transfo 1
S14 = V1*conj((V4-V1)*Y(1,4))
S41 = V4*conj((V1-V4)*Y(1,4))
ST1 = S14+S41
disp(' ');
%Ligne 2-5, Transfo 2
S25 = V2*conj((V5-V2)*Y(5,2))
S52 = V5*conj((V2-V5)*Y(5,2))
ST2 = S25+S52
disp(' ')
%Ligne 3-6, Transfo 3
S36 = V3*conj((V6-V3)*Y(3,6))
S63 = V6*conj((V3-V6)*Y(3,6))
ST3 = S36 + S63
disp(' ');
%Ligne 4-5
S45 = V4*conj((V5-V4)*Y(5,4))
S54 = V5*conj((V4-V5)*Y(5,4))
S454 = S45+S54
disp(' ');
%Ligne 4-6
S46 = V4*conj((V6-V4)*Y(4,6))
S64 = V6*conj((V4-V6)*Y(4,6))
S464 = S46 + S64
disp(' ');
%Ligne 5-7
S57 = V5*conj((V7 - V5)*Y(7,5))
S75 = V7*conj((V5 - V7)*Y(7,5))
S575 = S75 + S57
disp(' ');
%Ligne 6 - 7
S67 = V6*conj((V7 - V6)*Y(7,6))
S76 = V7*conj((V6 - V7)*Y(7,6))
S676 = S76 + S67
disp(' ');
%Bilan
I1 = Y(1,1)*V1cor+Y(1,2)*V2cor+Y(1,3)*V3+Y(1,4)*V4+Y(1,5)*V5+Y(1,6)*V6+Y(1,7)*V7
S1 = V1cor*conj(I1)
I2 = Y(2,1)*V1cor+Y(2,2)*V2cor+Y(2,3)*V3+Y(2,4)*V4+Y(2,5)*V5+Y(2,6)*V6+Y(2,7)*V7
S2 = V2cor*conj(I2)
I3 = Y(3,1)*V1cor+Y(3,2)*V2cor+Y(3,3)*V3+Y(3,4)*V4+Y(3,5)*V5+Y(3,6)*V6+Y(3,7)*V7
S3 = V3*conj(I3)
I4 = Y(4,1)*V1cor+Y(4,2)*V2cor+Y(4,3)*V3+Y(4,4)*V4+Y(4,5)*V5+Y(4,6)*V6+Y(4,7)*V7
S4 = V4*conj(I4)
I5 = Y(5,1)*V1cor+Y(5,2)*V2cor+Y(5,3)*V3+Y(5,4)*V4+Y(5,5)*V5+Y(5,6)*V6+Y(5,7)*V7
S5 = V5*conj(I5)
I6 = Y(6,1)*V1cor+Y(6,2)*V2cor+Y(6,3)*V3+Y(6,4)*V4+Y(6,5)*V5+Y(6,6)*V6+Y(6,7)*V7
S6 = V6*conj(I6)
I7 = Y(7,1)*V1cor+Y(7,2)*V2cor+Y(7,3)*V3+Y(7,4)*V4+Y(7,5)*V5+Y(7,6)*V6+Y(7,7)*V7
S7 = V7*conj(I7)
% Puissances générées
Sg = S1 + S2 + S3
% Pertes en lignes
Sl = S454 + S464 + S676 +S575
% Pertes Transformateurs
ST=ST1+ST2+ST3
% Pertes totales (lignes + transformateurs)
Sp = Sl + ST
% Puissance demandée (Pertes totales + charge)
Sd = Sp - S7
% Bilan
S = Sg - Sd
disp(' ');
% Tensions de chaque bus
disp([abs(V1cor) angle(V1cor)*180/pi])
disp([abs(V2cor) angle(V2cor)*180/pi])
disp([abs(V3) angle(V3)*180/pi])
disp([abs(V4) angle(V4)*180/pi])
disp([abs(V5) angle(V5)*180/pi])
disp([abs(V6) angle(V6)*180/pi])
disp([abs(V7) angle(V7)*180/pi])
