UNVERSITE DES SCIENCES ET DE LA TECHNOLOGIE HAOUARI BOUMEDIENE
Faculté d’Electronique et d’Informatique/Département Electrotechnique
. MASTER ME & EI (Module : Réseaux de transport et de distribution d’énergie électrique) .
SERIE D'EXERCICES N°1 : Ecoulement de puissance
Exercice N°1 : Soit le schéma unifilaire du réseau de transport électrique suivant :
104.8 200 MW
8 MW
150
76.23 MVar200+j
1
MVar
150 2
Bus N°1 QG1
Référence
Slack bus
100.4
MWPG3= 100.8
100 MW MVarQ
G3
Branche (nœud à nœud) Impédance sérieZs (pu) admittance shunt (Y /2) (pu)
1-2 0.02 + j 0.15 j 0.05
1-3 0.01 + j 0.20 j 0.10
2-3 0.01 + j 0.10 j 0.05
N° Nœud Tensions des nœuds Puissances générées Puissances consommées
Magnitude (pu) Angle (deg) P (MW) Q (MVAr) P (MW) Q(MVAr)
1 1.05 0 - - 0 0
2 - - 0 0 200 150
3 1.04 - 100 - 0 0
On prend 100 MVA comme puissance base et une précision égale à 0.01.
1°)- Calculer les admittances séries des différents branches du réseau ?
2°)- Donner l’algorithme de formation des éléments de la matrice admittance Ybus ?
3°)- Former la matrice admittance Ybus ?
4°)- En utilisant la méthode de Gauss-Seidel. Déterminer les tensions aux nœuds 2 et 3 ?
5°)- Trouver les puissances active et réactive du nœud balancier (slack bus) ?
6°)- Calculer la puissance réactive au nœud machine (3), PV bus ?
7°)- En déduire les puissances transitées et les pertes le long des différentes branches ?
.
[Link] &[Link] 2019-2020
UNVERSITE DES SCIENCES ET DE LA TECHNOLOGIE HAOUARI BOUMEDIENE
Faculté d’Electronique et d’Informatique/Département Electrotechnique
. MASTER ME & EI (Module : Réseaux de transport et de distribution d’énergie électrique) .
Exercice N°2 : Soit le schéma unifilaire du réseau de distribution électrique suivant :
4
3
Branche (nœud à nœud) Impédance série Zs (pu) admittance shunt (Y /2) (pu)
1-2 0.1 + j 0.3 j 0.10 au lieu de 1.0
2-3 0.07 + j 0.20 j 0.00
2-4 0.05 + j 0.15 j 0.00
N° Nœud Tensions des nœuds Puissances consommées
Magnitude (pu) Angle (deg) P (MW) Q(MVAr)
1 1.05 0 0 0
2 - - 5 2
3 - - 4 1.5
4 - - 6 3
On prend 100 MVA comme puissance base et une précision égale à 0.01.
1°)- Calculer les admittances séries des différents branches du réseau ?
2°)- Former la matrice admittance Ybus ?
3°)- En utilisant la méthode de Gauss-Seidel. Déterminer les tensions aux nœuds 2, 3 et 4 ?
5°)- Trouver les puissances active et réactive du nœud balancier (slack bus) ?
6°)- Calculer les puissances transitées et les pertes dans les différentes branches ?
7°)- Calculer les pertes active et réactive totales dans le réseau ?
.
[Link] &[Link] 2019-2020
UNVERSITE DES SCIENCES ET DE LA TECHNOLOGIE HAOUARI BOUMEDIENE
Faculté d’Electronique et d’Informatique/Département Electrotechnique
. MASTER ME & EI (Module : Réseaux de transport et de distribution d’énergie électrique) .
Solution de l’exercice N°1 :
1) Admittances séries des branches :
1
𝑦𝑠𝑖𝑗 = 𝑧
𝑠𝑖𝑗
1 1
𝑦𝑠12 = 𝑧 = 0.02+𝑗 0.15 = 0.8734 – j 6.5502 pu
𝑠12
1 1
𝑦𝑠13 = 𝑧 = 0.01+𝑗 0.2 = 0.2494 – j 4.9875 pu
𝑠13
1 1
𝑦𝑠23 = 𝑧 = 0.01+𝑗 0.1 = 0.9901 – j 9.9010 pu
𝑠23
2) Algorithme de formation des éléments de la matrice admittance Ybus :
𝑌𝑖𝑗 = −𝑦𝑠𝑖𝑗
𝑁
𝑌𝑖𝑖 = 𝑗 =1(𝑦𝑠𝑖𝑗 + 𝑦𝑠𝑖𝑗 )
𝑗 ≠𝑖
3) Formation de la matrice admittance Ybus :
Eléments diagonaux :
3
𝑌11 = 𝑗 =1(𝑦𝑠1𝑗 + 𝑦𝑠1𝑗 ) =𝑦𝑠12 + 𝑦𝑠12 +𝑦𝑠13 + 𝑦𝑠13
𝑗 ≠1
𝑌11 = (0.0+j0.05)+(0.8734–j 6.5502)+(0 + j0.10)+(0.2494 – j4.9875)
𝑌11 = 1.1228 -j11.3877 pu
3
𝑌22 = 𝑗 =1(𝑦𝑠2𝑗 + 𝑦𝑠2𝑗 ) =𝑦𝑠21 + 𝑦𝑠21 +𝑦𝑠23 + 𝑦𝑠23
𝑗 ≠2
𝑌22 = (0.0+j0.05)+(0.8734–j 6.5502)+(0 + j0.05)+( 0.9901 – j 9.9010)
𝑌22 =1.8635 -j16.3512 pu
3
𝑌33 = 𝑗 =1(𝑦𝑠3𝑗 + 𝑦𝑠3𝑗 ) =𝑦𝑠31 + 𝑦𝑠31 +𝑦𝑠32 + 𝑦𝑠32
𝑗 ≠3
𝑌33 = (0 + j0.10)+(0.2494 – j4.9875) +(0 + j0.05)+( 0.9901 – j 9.9010)
𝑌33 =1.2395 -14.7385 pu
Eléments non diagonaux :
𝑌12 = 𝑌21 = −𝑦𝑠12 =-(0.8734 – j 6.5502)= - 0.8734 + j 6.5502 pu
𝑌13 = 𝑌31 = −𝑦𝑠13 =-(0.2494 – j 4.9875)= - 0.2494 + j 4.9875 pu
𝑌23 = 𝑌32 − 𝑦𝑠23 =-(0.9901 – j 9.9010)= - 0.9901 + j 9.9010 pu
La matrice admittance du réseau est :
+ 1.1228 − j11.3877 − 0.8734 + j 6.5502 − 0.2494 + j 4.9875
𝑌 = − 0.8734 + j 6.5502 + 1.8635 − j16.3512 − 0.9901 + j 9.9010
− 0.2494 + j 4.9875 − 0.9901 + j 9.9010 + 1.2395 − 14.7385
.
[Link] &[Link] 2019-2020
UNVERSITE DES SCIENCES ET DE LA TECHNOLOGIE HAOUARI BOUMEDIENE
Faculté d’Electronique et d’Informatique/Département Electrotechnique
. MASTER ME & EI (Module : Réseaux de transport et de distribution d’énergie électrique) .
4) Détermination des tensions aux nœuds 2 et 3 par la méthode Gauss-Seidel :
(𝑘+1) 1 𝑃𝑖 −𝑗 𝑄𝑖 𝑖−1 𝑘+1 𝑁 𝑘
𝑉𝑖 = ∗(𝑘) −( 𝑗=1(𝑌𝑖𝑗 𝑉𝑗 + 𝑗=𝑖+1(𝑌𝑖𝑗 𝑉𝑗 )
𝑌𝑖𝑖 𝑉𝑖
(𝑘+1) 1 𝑃𝑖 −𝑗 𝑄𝑖 (𝑘+1) (𝑘+1) (𝑘) (𝑘)
𝑉𝑖 = 𝑌𝑖𝑖 ∗(𝑘) − (𝑌𝑖1 𝑉1 + ⋯ + 𝑌𝑖𝑖 −1 𝑉𝑖−1 + 𝑌𝑖𝑖 +1 𝑉𝑖+1 + ⋯ + 𝑌𝑖𝑁 𝑉𝑁 )
𝑉𝑖
Avant de commencer les calculs des tensions, il faut :
Classifier les nœuds, donner une solution initiale et calculer les puissances injectées :
(0) (𝑓𝑖𝑛𝑎𝑙𝑒 )
Nœud N°1 : Nœud de référence, 𝑉1 = 𝑉1 = 1.05 + 𝑗 0.0 = 1.05 𝑒 𝑗 0 pu
(0) 200+𝑗 150
Nœud N°2 : Nœud PQ, 𝑉2 = 1 + 𝑗0.0 = 1 𝑒 𝑗 0 pu, 𝑆𝐼2 = 𝑆𝐺2 − 𝑆𝐶2 = 0 − = −2 − 𝑗 1.5 pu
100
(0)
Nœud N°3 : Nœud PV, 𝑉3 = 1.04 + 𝑗0.0 = 1.04 𝑒 𝑗 0 pu, 𝑃𝐼3 = 𝑃𝐺3 − 𝑃𝐶3 = 100/100 − 0 = 1 pu
Calcul des tensions nodales à l’itération N°1 :
Nœud N°1 : La tension ne change pas (Référence)
Nœud N°2 : PQ
(1) 1 𝑃𝐼2 −𝑗 𝑄𝐼2 2−1 1 3 0
𝑉2 = ∗(0) −( 𝑗=1 (𝑌2𝑗 𝑉𝑗 + 𝑗=2+1(𝑌2𝑗 𝑉𝑗 )
𝑌22 𝑉2
(1) 1 𝑃𝐼2 −𝑗 𝑄𝐼2 (1) (0)
𝑉2 = 𝑌22 ∗(0) − (𝑌21 𝑉1 + 𝑌23 𝑉3 )
𝑉2
(1) 1 −2−𝑗(−1.5)
𝑉2 = − − 0.8734 + j 6.5502 1.05 + − 0.9901 + j 9.9010 (1.04))
1.8635 −j16.3512 1
(1)
𝑉2 = 0.9460 – j 0.1111 = 0.9525 𝑒 𝑗 (−0.1169 𝑟𝑑 ) 𝑝𝑢
Nœud N°3 : PV
(1)
(1) 1 𝑃𝐼3 −𝑗 𝑄𝐼3 3−1 1
𝑉3 = 𝑌33 ∗(0) −( 𝑗=1 (𝑌3𝑗 𝑉𝑗 + 0)
𝑉3
(1)
(1) 1 𝑃𝐼3 −𝑗 𝑄𝐼3 (1) (1)
𝑉3 = 𝑌33 ∗(0) − (𝑌31 𝑉1 + 𝑌32 𝑉2 )
𝑉3
On n’a pas la puissance réactive injectée, il faut la calculer en premier lieu par :
(𝑘+1) ∗ 𝑘 𝑖−1 𝑘+1 𝑁 𝑘
𝑄𝐼𝑖 = −𝐼𝑚 𝑉𝑖 𝑗=1(𝑌𝑖𝑗 𝑉𝑗 + 𝑗=𝑖(𝑌𝑖𝑗 𝑉𝑗
(1) ∗ 0 2 1 3 0
𝑄𝐼3 = −𝐼𝑚 𝑉3 𝑗=1(𝑌3𝑗 𝑉𝑗 + 𝑗=3(𝑌3𝑗 𝑉𝑗
(1) ∗ 0 1 1 0
𝑄𝐼3 = −𝐼𝑚 𝑉3 𝑌31 𝑉1 + 𝑌32 𝑉2 + 𝑌33 𝑉3
𝑄3 = −𝐼𝑚 (1.04) −0.2494 + 𝑗 4.9875 1.05 + −0.9901 + 𝑗 9.9010 0.9460 – j 0.1111 + 1.2395 − 14.7385 (1.04)
𝑄3 = 0.6397 𝑝𝑢
(1) 1 1−𝑗 0.6397
𝑉3 = − ((− 0.2494 + j 4.9875)(1.05) + (− 0.9901 + j 9.9010)(0.9525 ej (−0.1169 rd ) ))
1.2395 −14.7385 1.04
(1)
𝑉3 = 1.0387 – 𝑗 0.0154 𝑝𝑢 = 1.0388 𝑒𝑗 (−0.0148 𝑟𝑑) 𝑝𝑢
.
[Link] &[Link] 2019-2020
UNVERSITE DES SCIENCES ET DE LA TECHNOLOGIE HAOUARI BOUMEDIENE
Faculté d’Electronique et d’Informatique/Département Electrotechnique
. MASTER ME & EI (Module : Réseaux de transport et de distribution d’énergie électrique) .
(𝟏)
Correction du module de la tension 𝐕𝟑 :
(1)
(1) 𝑠𝑝 é𝑐𝑖𝑓𝑖 é 𝑉3
𝑉3 = 𝑉3 1 = 1.04 𝑒 𝑗 (−0.0148 𝑟𝑑 ) 𝑝𝑢 = 1.0399 − 0.0154 𝑝𝑢
𝑎𝑏𝑠 (𝑉3 )
Test de convergence :
𝑑𝑉1 = 0 (Nœud de référence)
(1) (0)
𝑑𝑉2 = 𝑉2 − 𝑉2 = 0.9460 – j 0.1111 − 1 = −0.054 − 𝑗 0.1111 𝑝𝑢
(1) (0)
𝑑𝑉3 = 𝑉3 − 𝑉3 = 1.0387 – 𝑗 0.0154 − 1.04 = −0.0013 – 𝑗 0.0154 𝑝𝑢
Erreur=max (0, 0.054, 0.1111, 0.0013, 0.0154)=0.1111 > Précision=0.01 Passer à l’itération N°2.
Calcul des tensions nodales à l’itération N°2 :
Nœud N°1 : La tension ne change pas (Référence)
Nœud N°2 : PQ
(2) 1 𝑃𝐼2 −𝑗 𝑄𝐼2 2−1 2 3 1
𝑉2 = ∗(1) −( 𝑗=1 (𝑌2𝑗 𝑉𝑗 + 𝑗=2+1(𝑌2𝑗 𝑉𝑗 )
𝑌22 𝑉2
(2) 1 𝑃𝐼2 −𝑗 𝑄𝐼2 (2) (1)
𝑉2 = 𝑌22 ∗(1) − (𝑌21 𝑉1 + 𝑌23 𝑉3 )
𝑉2
(2) 1 −2−𝑗 (−1.5)
𝑉2 = − − 0.8734 + j 6.5502 1.05 + − 0.9901 + j 9.9010 (1.04 𝑒 𝑗 (−0.0148 𝑟𝑑 ) ))
1.8635 −j16.3512 0.9525 𝑒 𝑗 (+0.1169 𝑟𝑑 )
(1)
Le signe positif à cause du conjugué de 𝑉2 Module après correction
(2)
𝑉2 = 0.9271 – j 0.1123 = 0.9339 𝑒 𝑗 (−0.1206 𝑟𝑑 ) 𝑝𝑢
Nœud N°3 : PV
(2)
(2) 1 𝑃 𝐼3 −𝑗 𝑄𝐼3 3−1 2
𝑉3 = ∗(1) −( 𝑗=1 (𝑌3𝑗 𝑉𝑗 + 0)
𝑌33 𝑉3
(2)
(2) 1 𝑃 𝐼3 −𝑗 𝑄𝐼3 (2) (2)
𝑉3 = ∗(1) − (𝑌31 𝑉1 + 𝑌32 𝑉2 )
𝑌33 𝑉3
On n’a pas la puissance réactive injectée, il faut la calculer en premier lieu par :
(𝑘+1) 𝑘 𝑘+1 𝑘
𝑄𝐼𝑖 = −𝐼𝑚 𝑉𝑖∗ 𝑖−1
𝑗 =1(𝑌𝑖𝑗 𝑉𝑗 + 𝑁
𝑗 =𝑖 (𝑌𝑖𝑗 𝑉𝑗
(2) 1 2 2 3 1
𝑄𝐼3 = −𝐼𝑚 𝑉3∗ 𝑗 =1(𝑌3𝑗 𝑉𝑗 + 𝑗 =3 (𝑌3𝑗 𝑉𝑗
(2) 1 2 2 1
𝑄𝐼3 = −𝐼𝑚 𝑉3∗ 𝑌31 𝑉1 + 𝑌32 𝑉2 + 𝑌33 𝑉3
𝑄3 = −𝐼𝑚 (1.04 𝑒 𝑗 (+0.0148 𝑟𝑑 ) ) −0.2494 + 𝑗 4.9875 1.05 + −0.9901 + 𝑗 9.9010 0.9339 𝑒 𝑗 (−0.1206 𝑟𝑑 ) +
1.2395 −14.7385(1.04 𝑒𝑗 −0.0148 𝑟𝑑))
𝑄3 = 0.8357 𝑝𝑢
(2) 1 1−𝑗 0.8357
𝑉3 = − ((− 0.2494 + j 4.9875)(1.05) + (− 0.9901 + j 9.9010)(0.9339 𝑒 𝑗 (−0.1206 𝑟𝑑 ) ))
1.2395 −14.7385 1.04 𝑒 𝑗 (+0.0148 𝑟𝑑 )
(2)
𝑉3 = 1.0396 – 𝑗 0.0184 𝑝𝑢 = 1.0398 𝑒 𝑗 (−0.0177 𝑟𝑑 ) 𝑝𝑢
(𝟐)
Correction du module de la tension 𝑽𝟑 :
(1)
(2) 𝑠𝑝 é𝑐𝑖𝑓𝑖 é 𝑉3
𝑉3 = 𝑉3 1 = 1.04 𝑒 𝑗 (−0.0177 𝑟𝑑 ) 𝑝𝑢
𝑎𝑏𝑠 (𝑉3 )
.
[Link] &[Link] 2019-2020
UNVERSITE DES SCIENCES ET DE LA TECHNOLOGIE HAOUARI BOUMEDIENE
Faculté d’Electronique et d’Informatique/Département Electrotechnique
. MASTER ME & EI (Module : Réseaux de transport et de distribution d’énergie électrique) .
Test de convergence :
𝑑𝑉1 = 0 + 𝑗 0 (Nœud de référence)
(2) (1)
𝑑𝑉2 = 𝑉2 − 𝑉2 = 0.9271 – j 0.1123 − 0.9460 – j 0.1111 = −0.0189 – 𝑗 0.0013 𝑝𝑢
(2) (1)
𝑑𝑉3 = 𝑉3 − 𝑉3 = 1.0396 – 𝑗 0.0184 − 1.0399 − 0.0154 = −0.0003 – 𝑗 0.0030 𝑝𝑢
Erreur=max (0, 0.0189, 0.0013, 0.0003, 0.003)=0.0189 > Précision=0.010 Passer à l’itération N°3.
Calcul des tensions nodales à l’itération N°3 :
Nœud N°1 : La tension ne change pas (Référence)
Nœud N°2 : PQ
(3) 1 𝑃 𝐼2 −𝑗 𝑄𝐼2 2−1 3 3 2
𝑉2 = ∗(2) −( 𝑗 =1 (𝑌2𝑗 𝑉𝑗 + 𝑗 =2+1 (𝑌2𝑗 𝑉𝑗 )
𝑌22 𝑉2
(3) 1 𝑃 𝐼2 −𝑗 𝑄𝐼2 (3) (2)
𝑉2 = ∗(2) − (𝑌21 𝑉1 + 𝑌23 𝑉3 )
𝑌22 𝑉2
(3)
𝑉2 =
1 −2−𝑗 (−1.5)
− − 0.8734 + j 6.5502 1.05 +
1.8635 −j16.3512 0.9339 𝑒 𝑗 (+0.1206 𝑟𝑑 )
− 0.9901 + j 9.9010 (1.04 𝑒𝑗 −0.0177 𝑟𝑑)))
(3)
𝑉2 = 0.9248 – j 0.1157 = 0.932 𝑒 𝑗 (−0.1245 𝑟𝑑 ) 𝑝𝑢
Nœud N°3 : PV
(3)
(3) 1 𝑃 𝐼3 −𝑗 𝑄𝐼3 3−1 3
𝑉3 = ∗(2) −( 𝑗 =1 (𝑌3𝑗 𝑉𝑗 + 0)
𝑌33 𝑉3
(3)
(3) 1 𝑃 𝐼3 −𝑗 𝑄𝐼3 (3) (3)
𝑉3 = ∗(2) − (𝑌31 𝑉1 + 𝑌32 𝑉2 )
𝑌33 𝑉3
On n’a pas la puissance réactive injectée, il faut la calculer en premier lieu par :
(𝑘+1) 𝑘 𝑘 +1 𝑘
𝑄𝐼𝑖 = −𝐼𝑚 𝑉𝑖∗ 𝑖−1
𝑗 =1 (𝑌𝑖𝑗 𝑉𝑗 + 𝑁
𝑗 =𝑖 (𝑌𝑖𝑗 𝑉𝑗
(3) 2 2 3 3 2
𝑄𝐼3 = −𝐼𝑚 𝑉3∗ 𝑗 =1(𝑌3𝑗 𝑉𝑗 + 𝑗 =3(𝑌3𝑗 𝑉𝑗
(3) 2 3 3 2
𝑄𝐼3 = −𝐼𝑚 𝑉3∗ 𝑌31 𝑉1 + 𝑌32 𝑉2 + 𝑌33 𝑉3
(3)
𝑄𝐼3 =
−𝐼𝑚 (1.04 𝑒 𝑗 (+0.0177 𝑟𝑑 ) ) −0.2494 + 𝑗 4.9875 1.05 + −0.9901 + 𝑗 9.9010 0.932 𝑒 𝑗 (−0.1245 𝑟𝑑 ) +
1.2395 −14.7385(1.04 𝑒𝑗 −0.0177 𝑟𝑑))
(3)
𝑄𝐼3 = 0.8562 𝑝𝑢
(3)
𝑉3 =
1 1−𝑗 0.8562
− ((− 0.2494 + j 4.9875)(1.05) + (− 0.9901 +
1.2395 −14.7385 1.04 𝑒 𝑗 (+0.0177 𝑟𝑑 )
j 9.9010)(0.932 𝑒𝑗 −0.1245 𝑟𝑑)))
(3)
𝑉3 = 1.0396 – 𝑗 0.021 𝑝𝑢 = 1.0398 𝑒 𝑗 (−0.0202 𝑟𝑑 ) 𝑝𝑢
(𝟑)
Correction du module de la tension 𝑽𝟑 :
.
[Link] &[Link] 2019-2020
UNVERSITE DES SCIENCES ET DE LA TECHNOLOGIE HAOUARI BOUMEDIENE
Faculté d’Electronique et d’Informatique/Département Electrotechnique
. MASTER ME & EI (Module : Réseaux de transport et de distribution d’énergie électrique) .
(2)
(3) 𝑠𝑝 é𝑐𝑖𝑓𝑖 é 𝑉3
𝑉3 = 𝑉3 2 = 1.04 𝑒 𝑗 (−0.0202 𝑟𝑑 ) 𝑝𝑢
𝑎𝑏𝑠 (𝑉3 )
Test de convergence :
𝑑𝑉1 = 0 (Nœud de référence)
(3) (2)
𝑑𝑉2 = 𝑉2 − 𝑉2 = 0.9248 – j 0.1157 − 0.9271 – j 0.1123 = −0.0023 – 𝑗 0.0034 𝑝𝑢
(3) (2)
𝑑𝑉3 = 𝑉3 − 𝑉3 = 1.0396 – 𝑗 0.0210 − 1.0398 – 𝑗 0.0184 = −0.0003 − 𝑗 0.0026 𝑝𝑢
Erreur=max (0, 0.0023, 0.0034, 0.0003, 0.0026)=0.0034 < Précision=0.01 Convergence à l’itération N°3.
La solution de tension est :
(3)
𝑉1 = 1.05 𝑝𝑢 𝑒 𝑗 (0 𝑟𝑑 )
(3)
𝑉2 = 0.9248 – j 0.1157 = 0.932 𝑝𝑢 𝑒 𝑗 (−0.1245 𝑟𝑑 ) = 0.932 𝑝𝑢 𝑒 𝑗 (−7.1325°)
(3)
𝑉3 = 1.0396 – 𝑗 0.021 𝑝𝑢 = 1.0398 𝑝𝑢 𝑒 𝑗 (−0.0202 𝑟𝑑 ) = 1.0398 𝑝𝑢 𝑒 𝑗 (−1.159°)
5) Puissances active et réactive du nœud balancier (Slack Bus) :
3 3 3 3 3 3 3
𝑃𝐼1 = 𝑃𝐺1 = 𝑅𝑒𝑎𝑙 𝑉1∗ 𝑗 =1 (𝑌𝑖𝑗 𝑉𝑗 = 𝑅𝑒𝑎𝑙 𝑉1∗ 𝑌11 𝑉1 + 𝑌12 𝑉2 + 𝑌13 𝑉3
𝑘 3 3 3 3 3 3
𝑄𝐼1 = 𝑄𝐺1 = −𝐼𝑚 𝑉1∗ 𝑗 =1 (𝑌1𝑗 𝑉𝑗 = −𝐼𝑚 𝑉1∗ 𝑌11 𝑉1 + 𝑌12 𝑉2 + 𝑌13 𝑉3
𝑃𝐼1 = 𝑃𝐺1 = 𝑅𝑒𝑎𝑙 ((1.05) 1.1227 – 𝑗 11.3877 1.05 + −0.8734 + 𝑗 6.5502 0.932 𝑒 𝑗 (−0.1245 𝑟𝑑 ) +
−0.2494 + 𝑗 4.9875(1.0398 𝑒𝑗 −0.0202 𝑟𝑑)))
𝑃𝐼1 = 𝑃𝐺1 = 𝑅𝑒𝑎𝑙 1.0236 − 𝑗 0.6389 = 1.0236 𝑝𝑢 = +123.6 𝑀𝑊
𝑄𝐼1 = 𝑄𝐺1 = −𝐼𝑚 ((1.05) 1.1227 – 𝑗 11.3877 1.05 + −0.8734 + 𝑗 6.5502 0.932 𝑒 𝑗 (−0.1245 𝑟𝑑 ) +
−0.2494 + 𝑗 4.9875(1.0398 𝑒𝑗 −0.0202 𝑟𝑑)))
𝑄𝐼1 = 𝑄𝐺1 = −𝐼𝑚 1.0236 − 𝑗 0.6389 = −0.6389 𝑝𝑢 = +63.89 𝑀𝑉𝑎𝑟
6) Puissance réactive au nœud générateur 3 (PV Bus) :
3 3 3 3
𝑄𝐼3 = 𝑄𝐺3 = −𝐼𝑚 𝑉3∗ 𝑌31 𝑉1 + 𝑌32 𝑉2 + 𝑌33 𝑉3
𝑄𝐺3 =
−𝐼𝑚 (1.0398 𝑒 𝑗 (+0.0202 𝑟𝑑 ) ) −0.2494 + 𝑗 4.9875 1.05 +
−0.9901 +𝑗 9.90100.932 𝑒𝑗 −0.1245 𝑟𝑑)+ 1.2395 −14.7385(1.0398 𝑒𝑗 −0.0202 𝑟𝑑))
𝑄𝐼3 = 𝑄𝐺3 = −𝐼𝑚 1.0019 − 𝑗 0.8535 = −0.8535 𝑝𝑢 = + 85.35 𝑀𝑉𝑎𝑟
7) Puissances transitées et les pertes le long des différentes branches :
Branche i-j :
𝐼𝑖𝑗 = 𝐼𝑠é𝑟𝑖𝑒 −𝑖𝑗 + 𝐼𝑠𝑢𝑛𝑡 −𝑖𝑗 = 𝑉𝑖 − 𝑉𝑗 𝑦𝑠𝑖𝑗 + 𝑉𝑖 𝑦𝑠𝑖𝑗
𝑆𝑖𝑗 = 𝑃𝑖𝑗 + 𝑗𝑄𝑖𝑗 = 𝑉𝑖 𝐼𝑖𝑗∗
.
[Link] &[Link] 2019-2020
UNVERSITE DES SCIENCES ET DE LA TECHNOLOGIE HAOUARI BOUMEDIENE
Faculté d’Electronique et d’Informatique/Département Electrotechnique
. MASTER ME & EI (Module : Réseaux de transport et de distribution d’énergie électrique) .
Branche 1-2 :
𝐼12 = 𝐼𝑠é𝑟𝑖𝑒 −12 + 𝐼𝑠𝑢𝑛𝑡 −12 = 𝑉1 − 𝑉2 𝑦𝑠12 + 𝑉1 𝑦𝑠12
𝐼12 = 1.05 − (0.9248 – j 0.1157) 0.8734 − 6.5502 + 1.05 𝑗 0.05 = 0.8674 − 𝑗 0.6667 𝑝𝑢
∗
𝑆12 = 𝑃12 + 𝑗𝑄12 = 𝑉1 𝐼12
𝑆12 = 1.05 0.8674 + 0.6667 = 0.9107 + 𝑗 0.7 𝑝𝑢 = +91.07 𝑀𝑊 + 𝑗 70 𝑀𝑉𝑎𝑟
Branche 2-1 :
𝐼21 = 𝐼𝑠é𝑟𝑖𝑒 −21 + 𝐼𝑠𝑢𝑛𝑡 −21 = 𝑉2 − 𝑉1 𝑦𝑠12 + 𝑉2 𝑦𝑠12 = −0.8616 + 𝑗 0.7654 𝑝𝑢
𝑆21 − 0.8853 – 𝑗 0.6082 = −88.53 𝑀𝑊 − 𝑗 60.82 𝑀𝑉𝑎𝑟
Perte dans la branche 12 :
PertesBranche12 = 𝑆12 + 𝑆21 = +91.07 𝑀𝑊 + 𝑗 70 𝑀𝑉𝑎𝑟 + (−88.53 𝑀𝑊 − 𝑗 60.82 𝑀𝑉𝑎𝑟)
PertesBranche12 = 2.54 𝑀𝑊 + 9.19 𝑀𝑉𝑎𝑟
Branche 1-3 :
𝐼13 = 𝐼𝑠é𝑟𝑖𝑒 −13 + 𝐼𝑠𝑢𝑛 𝑡−13 = 𝑉1 − 𝑉3 𝑦𝑠13 + 𝑉1 𝑦𝑠13
𝐼13 = 1.05 − 1.0396 – 𝑗 0.021 0.2494 – 𝑗 4.9875 + 1.05 𝑗 0.1 = 0.1075 𝑝𝑢 + 𝑗 0.0582 𝑝𝑢
∗
𝑆13 = 𝑃13 + 𝑗𝑄13 = 𝑉1 𝐼13
𝑆13 = 1.05 0.1075 − 𝑗 0.0582 = 0.1129 𝑝𝑢 – 𝑗 0.0611 𝑝𝑢 = 11.29 𝑀𝑊 − 𝑗 6.11 𝑀𝑉𝑎𝑟
Branche 3-1 :
𝐼31 = 𝐼𝑠é𝑟𝑖𝑒 −31 + 𝐼𝑠𝑢𝑛𝑡 −31 = 𝑉3 − 𝑉1 𝑦𝑠13 + 𝑉3 𝑦𝑠13
𝐼31 = 1.0396 – 𝑗 0.021 − 1.05 0.2494 – 𝑗 4.9875 + 1.0396 – 𝑗 0.021 𝑗 0.1 = −0.1054 + 𝑗 0.1507
∗
𝑆31 = 𝑃31 + 𝑗𝑄31 = 𝑉3 𝐼31
𝑆13 = 1.0383 − 𝑗 0.0224 −0.1054 − 0.1507 = −0.1127𝑝𝑢 − 0.1545 𝑝𝑢 = −11.27 𝑀𝑊 − 𝑗 15.45 𝑀𝑉𝑎𝑟
Perte dans la branche 13 :
PertesBranche13 = 𝑆13 + 𝑆31 = 0.1129 𝑝𝑢 – 𝑗 0.0611 𝑝𝑢 + −0.1127 𝑝𝑢 − 0.1545 𝑝𝑢
PertesBranche13 = 0.0001 𝑝𝑢 – 𝑗 0.2156 𝑝𝑢 = 0.01 𝑀𝑊 − 21.56 𝑀𝑉𝑎𝑟
Branche 2-3 :
𝐼23 = 𝐼𝑠é𝑟𝑖𝑒 −23 + 𝐼𝑠𝑢𝑛𝑡 −23 = 𝑉2 − 𝑉3 𝑦𝑠23 + 𝑉2 𝑦𝑠23
𝐼23 = 0.9248 – 𝑗 0.1157 − 1.0396 – 𝑗 0.021 0.9901 – 𝑗 9.9010 + 1.0396 – 𝑗 0.021 𝑗 0.05 = −1.0454 + 𝑗1.0891
∗
𝑆23 = 𝑃23 + 𝑗𝑄23 = 𝑉2 𝐼23
𝑆23 = 0.9248 – j 0.1157 −1.0454 − 𝑗1.0891 = −1.0928 − 𝑗 0.8862 𝑝𝑢 = −109.28 𝑀𝑊 − 𝑗 88.62 𝑀𝑉𝑎𝑟
Branche 3-2 :
𝐼32 = 𝐼𝑠é𝑟𝑖𝑒 −32 + 𝐼𝑠𝑢𝑛𝑡 −32 = 𝑉3 − 𝑉2 𝑦𝑠23 + 𝑉3 𝑦𝑠23
𝐼32 = 1.0396 – 𝑗 0.021 − 0.9248 – j 0.1157 0.9901 − 𝑗9.901 + 1.0396 – 𝑗 0.021 𝑗 0.05 = 1.0522 − 𝑗 0.9909
∗
𝑆32 = 𝑃32 + 𝑗𝑄32 = 𝑉3 𝐼32
𝑆32 = 1.0396 – 𝑗 0.021 1.0522 + 𝑗 0.9909 = 1.1147𝑝𝑢 + 𝑗 1.008𝑝𝑢 = 111.47 𝑀𝑊 + 𝑗 100.8 𝑀𝑉𝑎𝑟
Perte dans la branche 23 :
.
[Link] &[Link] 2019-2020
UNVERSITE DES SCIENCES ET DE LA TECHNOLOGIE HAOUARI BOUMEDIENE
Faculté d’Electronique et d’Informatique/Département Electrotechnique
. MASTER ME & EI (Module : Réseaux de transport et de distribution d’énergie électrique) .
PertesBranche23 = 𝑆23 + 𝑆32 = −109.28 𝑀𝑊 − 𝑗 88.62 𝑀𝑉𝑎𝑟 + 111.47 𝑀𝑊 + 𝑗 100.8 𝑀𝑉𝑎𝑟
PertesBranche23 = 2.19 MW + j 12.18 MVar
Pertes dans le réseau :
PertesRéseau = PertesBranche = PertesBranche12 + PertesBranche13 + PertesBranche23
PertesRéseau = 2.54 𝑀𝑊 + 9.19 𝑀𝑉𝑎𝑟 + 0.01 𝑀𝑊 − 21.56 𝑀𝑉𝑎𝑟 + (2.19 MW + j 12.18 MVar )
PertesRéseau = 4.75 MW + j 2.0 MVar
Schématisation de l’écoulement de puissance :
104.48 MW 200 MW
-76.23 MVar 150 MVar
G
Bus N°1 1.05 𝑝𝑢 𝑒 𝑗 0 𝑑𝑒𝑔 0.932 𝑝𝑢 𝑒 𝑗 −7.13 °
Bus N°2 PQ
Référence
Pertes 1-2
-109.3 -88.62
+11.29 -6.11 +91.07 MW +2.54 MW -88.53 MW
MW MVar
MW MVar +70.00 MVar +9.19 MVar -60.82 MVar
Pertes 1-3 Pertes 2-3
+0.01 MW +2.19 MW
-21.56 MVar +12.18 MVar
-11.27 -15.45 +111.5 +100.8
MW MVar MW MVar
Bus N°3 PV
1.0398 𝑝𝑢 𝑒 𝑗 −1.159°
G
100.4 MW
100.86 MVar
.
[Link] &[Link] 2019-2020
