Réseaux électriques

Introduction

L’objectif est de maintenir la tension du réseau de transport dans des limites données.
Or ces tensions varient avec la puissance consommée. Des moyens de réglage télécommandés
de la puissance réactive maintiennent la tension dans une plage donnée.
ex:+ ou - 10% en BT aérien.
En bureau d’étude, on étudie la variation de tension en fonction de la variation de puissance.

=> Élaboration d’une stratégie.

En temps réel, si variation de production ou de consommation, on calcule les variations

prévisibles de tension

=> Réglage qui doit tenir compte de la capacité de transport.

Load Flow - durée 2h- G. Clerc et P. Auriol 1

 Réseaux électriques

Définition du problème
Déterminer complètement l’état du réseau dans des hypothèses données de production et de
consommation. bus

jeu de barre
bus de référence
lignes
PG0
(slack bus) 
1 PG2
G G
QG0 QG2
V2
2 3 Pc4, Qc4
Pc0, Qc0
V0 
G PG3
QG3
4
Pc3, Qc3
 V1
V3
Pc1, Qc1
PG1 consommation
G
QG1

Load Flow - durée 2h- G. Clerc et P. Auriol 2

 Réseaux électriques

Les paramètres

A chaque sommet sont définis :

• P puissance active injectée ou soutirée
• Q puissance réactive injectée ou soutirée
• V module de la tension
• q déphasage de la tension

On distingue trois types de sommets :

• Le sommet producteur (noté PV) pour lequel on connaît P et V. Les inconnues sont Q et q.
• Le sommet consommateur (noté PQ) pour lequel on connaît P et Q. Les inconnues sont V et q.
• Le sommet bilan ou balancier où on fixe V et q. Les inconnues sont P et Q.

Dans chaque branche, on détermine les puissances qui circulent et les pertes

Load Flow - durée 2h- G. Clerc et P. Auriol 3

 Réseaux électriques Qij E j  Vje
jq j Qji
Ei  Vi e jqi
Ei Yij Ej
Pij Pji
Mise en équations
Au sommet i Pi  jQi  Ei Ii * yij yji
 
avec I i  yii Ei   I ij  yii Ei    yij Ei   Yij Ei  E j 
  
j i  j i j i 
   
Pi  jQi  Ei I i *  Ei  Ei*  yii*   yij *     Ei  E j  Yij* 
*
Donc  
 * j i  j i  Attention aux
  2
Donc
 j i  j i

Pi  jQi   yii   yij  Vi   Yij *Vi e jqi Vi e  jqi  V j e j
 jq
 notations en yij et
Yij
 
 ii  ij   Gii  jH ii
Posons Yij  Gij  jH ij et Yii   y  y 
 j i 
Donc
  2
Pi  jQi   Gii  jH ii   Gij  jH ij Vi   Gij  jH ij ViV j cos(q i  q j )  j sin(q i  q j ) 

 j i  j i

Donc
  2
Pi  jQi   Gii   Gij Vi  ViV j Gij cos(q i  q j )  H ij sin(q i  q j ) 

 j i  j i

  2
 j  H ii   H ij Vi  j ViV j H ij cos(q i  q j )  Gij sin(q i  q j ) 

 j i  j i

Load Flow - durée 2h- G. Clerc et P. Auriol 4

 Réseaux électriques
Qij E j  Vje
jq j Qji
Ei  Vi e jqi
Ei Yij Ej
Pij Pji

yij yji

  2
Pi  jQi   Gii   Gij Vi  ViV j Gij cos(q i  q j )  H ij sin(q i  q j ) 

 j i  j i

  2
 j  H ii   H ij Vi  j ViV j H ij cos(q i  q j )  Gij sin(q i  q j ) 

 j i  j i

On en déduit

  2
Pi   Gii   Gij Vi  V V G i j ij cos(q i  q j )  H ij sin(q i  q j ) 
 j i  j i

  2
Qi   H ii   H ij Vi  ViV j H ij cos(q i  q j )  Gij sin(q i  q j ) 

 j i  j i

Load Flow - durée 2h- G. Clerc et P. Auriol 5

 Réseaux électriques

Un exemple simple

1
Gij  0 H ij   yij  y ji  0
X
On a donc :
1 1
 P  VV g sin q  P  VV g sin q
X X
1 1
Q  V 2  VVg cos(q )
X X
Vg est la tension du noeud balancier, V celle du nœud de charge, q sa phase/noeud balancier et
P+jQ sa consommation.

Load Flow - durée 2h- G. Clerc et P. Auriol 6

 Réseaux électriques

En éliminant q
V   2QX  V V
2 2
g
2 2
 
 X 2 P2  Q2  0

Pour avoir une solution (en V²), le discriminant doit être positif d’où la courbe dans le plan P,Q.
2
 PX 
 2   QX  0.25  0
V  V 2
 g  g

Load Flow - durée 2h - G. Clerc et P. Auriol

Load Flow - durée 2h- G. Clerc et P. Auriol 7

 Réseaux électriques

On obtient alors la solution (en posant y=V²) :

2 2
Vg Vg
V  y   QX   X 2 P 2  XQVg
2

2 4

PX
q   arcsin
VV g

Load Flow - durée 2h- G. Clerc et P. Auriol 8

 Réseaux électriques

Résolution
S’il y a s sommets, alors il y a 2s équations fondamentales à 2s inconnues (équations non
linéaires)
Supposons que tous les nœuds soient consommateurs.
Donc pour tout i, Pi et Qi sont donnés, Vi et qi sont inconnus (en fait on suppose connus les
puissances actives et réactives aux nœuds producteurs).
On cherche Vi et qi en chaque nœud pour satisfaire aux équations précédentes.
Le problème est non linéaire => méthode itérative.
1. Initialisation
On suppose Vi (0) et qi (0) connus (on se le fixe) en chaque sommet.
2. Calcul de l’erreur
On calcule (P0)i et (Q0)i en chaque sommet. Les écarts sont donnés par
DPi=Pi-(P0)i
DQi=Qi-(Q0)i
Load Flow - durée 2h- G. Clerc et P. Auriol 9
 Réseaux électriques

3. Correction
On a Pi  Pi (V1,..., Vs ,q1,..., qs )
Qi  Qi (V1,..., Vs ,q1,..., qs )
Pi P
Donc DPi   DV j   i Dq j
j V j j q j
Q Q
DQi   i DV j   i Dq j
j V j j q j
Pi
Pour j≠i  Vi Gij cos(q i  q j )  H ij sin(q i  q j ) 
V j
Pi
 ViV j Gij sin(q i  q j )  H ij cos(q i  q j ) 
q j
Qi
 Vi H ij cos(q i  q j )  Gij sin(q i  q j ) 
V j
Qi
 ViV j  Gij cos(q i  q j )  H ij sin(q i  q j ) 
q j

Load Flow - durée 2h- G. Clerc et P. Auriol 10

 Réseaux électriques

Pi  
D’où  2 Gii   Gij Vi  V j Gij cos(q i  q j )  H ij sin(q i  q j ) 
Vi  j i  j
Pi
  ViV j  Gij sin(q i  q j )  H ij cos(q i  q j ) 
(j=i)
q i j

Qi  
 2 H ii   H ij Vi  V j H ij cos(q i  q j )  Gij sin(q i  q j ) 
Vi  j i  j

Qi
  ViV j  H ij sin(q i  q j )  Gij cos(q i  q j ) 
q i j

Soit le vecteur correction tension DV  DV1,..., DVs t

Soit le vecteur correction de phase Dq  Dq1,..., Dqs t

DP  DP1,..., DPs t
Soit le vecteur correction de puissance
active
Soit le vecteur correction de puissance DQ  DQ1,..., DQs 
t

réactive
Load Flow - durée 2h- G. Clerc et P. Auriol 11
 Réseaux électriques

 Pi Pi 
On obtient alors le système linéaire  DP   q j V j   Dq 
 
DQ    Q Qi  DV 
  i
 q V j 
 j
En résolvant le système linéaire, on obtient : DV  DV1,..., DVs t

Dq  Dq1,..., Dqs t

Ce qui nous donne : q1  q1( 0)  Dq1

...
qs  qs ( 0)  Dqs
V1  V1( 0)  DV1
...
Vs  Vs ( 0)  DVs

Load Flow - durée 2h- G. Clerc et P. Auriol 12

 Réseaux électriques

On pose q i (0)  q i ,..., qs (0)  qs

Vi ( 0)  Vi ,..., Vs ( 0)  Vs

On recommence alors le processus en calculant Pi Pi Qi Qi

, , ,
les nouvelles valeurs de Vi q i Vi q i

 Pi Pi 
On résout le système  DP   q j V j   Dq 

DQ    Q Qi  DV 
  i
 q V j 
 j

Load Flow - durée 2h- G. Clerc et P. Auriol 13

 Réseaux électriques

Particularité du système
En général qi  q j  0

Gij  Hij
Pi P
 i
V j q j

Qi Q
 i
q j V j

Donc  Pi 
0 
 DP   q j Dq
  
DQ    Qi  DV 
 
 0 V j 

Le nœuds ne sont pas tous reliés entre eux généralement. Gij et Hij sont nuls pour i et j non
adjacents.
Peu de liaison => matrice très creuses => méthodes numériques spéciales
Load Flow - durée 2h- G. Clerc et P. Auriol 14
 Réseaux électriques

Nœuds producteurs
Aux nœuds producteurs, Qi est inconnue et Vi fixée.
[DP]=[ ][Dq] et [DP]=[ ][DV]

DVi=0 (Vi fixée) => on supprime la ième colonne

  2
Qi   H ii   H ij Vi  ViV j H ij cos(q i  q j )  Gij sin(q i  q j )  => on supprime la ième ligne

 j i  j i

On réduit le système à un ordre inférieur (= nb nœuds consommateurs) .

Load Flow - durée 2h- G. Clerc et P. Auriol 15

 Réseaux électriques

Nœud bilan ou balancier

Deux remarques :
1. On ne peut spécifier les puissances Pi et Qi à tous les nœuds. Le bilan de puissance complexe
du réseau s’écrit :
n

P  p
i 1
i

Q
i 1
i q

Où p et q représentent respectivement les pertes actives et réactives du réseau. Donner tous les Pi
et Qi reviendrait à spécifier les pertes. Or celles-ci sont fonction des courants dans les branches et
des tensions aux nœuds lesquelles ne sont a priori pas connues
N 1
Au nœud balancier
N 
P  P p
i 1
i

où les différents termes de la somme sont spécifiés dans les données, tandis que, p n’est connu
qu’`a l’issue du calcul
2. Par ailleurs il convient de calculer les déphasages par rapport à un nœud de référence.
On choisit le sommet où la puissance disponible est très forte : on fixe V et q.
Ce nœud fournit les pertes en ligne. Il n’y a pas d’inconnue à déterminer à ce nœud.

Load Flow - durée 2h- G. Clerc et P. Auriol 16

 Réseaux électriques

 
y ij *   Gij  jHij  Gij  jHij
Qij Qji
On choisit le sommet où la puissance disponible Ei yij Ej
Pij
est très forte : on fixe V et q. Pji
Ce nœud fournit les pertes en ligne.
y'ij y'ji
Calcul des tensions en ligne
Remarque
(Voir cours sur Mise en équation Cij
Yij  y ij réseau) y 'ij *   j   jhij
Bilan des puissances 2
Puissance transitant de i vers j Pertes parallèles à la ligne au nœud i

  
SijL  Ei y ij Ei  E j * pij// L  y 'ij Ei
2

Puissance transitant de j vers i Pertes parallèles à la ligne au nœud j

  
S ji L  E j y ji Ei  E j * p ji // L  y ' ji E j
2

    
On a donc Sij  Ei y ij Ei  E j *  y 'ij Ei * Ei Et    
S ji  E j y ji Ei  E j *  y ' ji E j * E j
Les Pertes globales sont données par : DS  Sij  S ji
Load Flow - durée 2h- G. Clerc et P. Auriol 17
 Réseaux électriques

Contraintes réelles
Qr
Pour les générateurs et les récepteurs
Courbe
Qmin(P )  Q  Qmax(P ) Q1 moyenne

Q2
Pr

P1
Zone de charge d’un récepteur
Quelques soient les nœuds, les variations de tension limitée (en V et q)
Certaines données sont incertaines ou inaccessibles
Il faut prendre en compte les contraintes économiques, les disponibilités, …

Problème d’optimisation. Compromis technico-économique

Load Flow - durée 2h- G. Clerc et P. Auriol 18

 Réseaux électriques

Bibliographie

[1] P. Auriol, Cours réseau

[2] W. D. Stevenson, Elements of power system analysis
[3] [Link]

Load Flow - durée 2h- G. Clerc et P. Auriol 19