POLITECHNIKA RZESZOWSKA

im. Ignacego Łukasiewicza

PRZESYŁ ENERGII ELEKTRYCZNEJ

Projekt 2

___ _____________________________________

Temat: Obliczanie rozpływów mocy u sieci rozdzielczej

___________ _____________________________

Biagio Marciano

EN-DI-3

L7

Rok akademicki 2013/2014

23.01.2014r.

Zadanie

Task 2:

For given simple power network calculate:

• voltages U (magnitude and phase angle of the phasor) in all nodes of electric network,

• active and reactive power (P_i, Q_i) in all nodes of electric network. Note: sign „+“ determines power supply, sign „-“ determines power comsumption,

• active and reactive power flows (P_ij, Q_ij) on the lines (values and directions),

• active power losses in power network.

Use the Newton-Rhapson method or Gauss-Seidl method. You can solve this task manually or using software MATLAB. In case of using MATLAB, please, print the script of the program.

line	Rk [Ω.km^-1]	Xk [Ω.km^-1]	Bk [µS.km^-1]	l [km]
1-2 (350_AlFe)	0,085	0,394	3,050	35
1-3 (210_AlFe)	0,130	0,400	2,900	40
2-3 (240_AlFe)	0,125	0,403	2,869	25

Rozwiązanie

Na podstawie poniższych wzorów obliczmy dla danych linii rezystancje ( R ) ,
reaktancje ( X ) , susceptancję ( B ) .

R = R_k • l  [Ω] 

X = X_k • l  [Ω] 

B = B_k • l  [μS] 

Wyniki przedstawiono w poniższej tabeli

Linia	R [Ω]	X [Ω]	B [µS]
1-2 (350_AlFe)	2,975	13,79	106,75
1-3 (210_AlFe)	5,2	16	116
2-3 (240_AlFe)	3,125	10,075	71,725

Schemat zastępczy sieci rozdzielczej

Dla naszych danych obliczamy admitancje według poniższych wzorów

• admitancje własne węzłów ( do węzłów wchodzą tylko linie )

$$= = \sum_{\text{jϵ}N_{i}}^{}\frac{1}{R_{\text{ij}} + j \bullet X_{\text{ij}}} + j \bullet \frac{B_{\text{ij}}}{2}\text{\ \ \ }\text{gdzie}\ {\ N}_{i}\ \in \{ 1,2,3\}$$

Y_ii = ||

• admitancje wzajemne gdy miedzy węzłami są linie

$= - \frac{1}{R_{\text{ij}} + j \bullet X_{\text{ij}}}$

Y_ij = ||

Teoria:

Celem naszego projektu jest wyznaczenie określenie wszystkich niewiadomych wielkości ||,δ, P,  Q w każdym węźle i na ich podstawie obliczenie rozpływów mocy . Zależności między mocami węzłowymi P ,  Q , a napięciami węzłowymi są następujące ( wartości zespolone ) :

$$P_{i} = \sum_{j = 1}^{n}{U_{i}{\bullet U}_{j} \bullet Y_{\text{ij}} \bullet \cos{(\delta_{i} - \delta_{j} - \alpha_{\text{ij}})}}\ \ \ \ \ \ i = 2\ ,3\ldots n\ $$

$$Q_{i} = \sum_{j = 1}^{n}{U_{i}{\bullet U}_{j} \bullet Y_{\text{ij}} \bullet \sin{(\delta_{i} - \delta_{j} - \alpha_{\text{ij}})}}\ \ \ \ \ \ i = 2\ ,3\ldots n\ $$

Zlinearyzowane równania mocowo-napięciowe przedstawione powyżej można przedstawić dla wszystkich węzłów sieci w postaci macierzowej :

$\begin{bmatrix} \begin{matrix} P_{1} \\ \vdots \\ P_{n} \\ \end{matrix} \\ Q_{1} \\ \vdots \\ Q_{m} \\ \end{bmatrix}\ \ = \ \begin{bmatrix} \left\lbrack \frac{\partial P}{\partial\delta} \right\rbrack & \left\lbrack \frac{\partial P}{\partial U} \right\rbrack \\ \left\lbrack \frac{\partial Q}{\partial\delta} \right\rbrack & \left\lbrack \frac{\partial Q}{\partial U} \right\rbrack \\ \end{bmatrix}\ \begin{bmatrix} \begin{matrix} \delta_{1} \\ \vdots \\ \delta_{n} \\ \end{matrix} \\ U_{1} \\ \vdots \\ U_{m} \\ \end{bmatrix}$

 Y          =             J                  X

+ algorytm możesz wrzucić z jego materiałów

Wartości uzyskane w wyniku obliczeń w Matlabie:

|| = 113000 V = 113 kV

α₂₂ = (0, 4498)≈0, 45

=( 11300+j887,01 ) V = ( 113+j0,887)kV

|| = 114020 V = 114 kV

α₃₃ = ( − 1.9682)≈2

=(  113950−j3916) V = ( 114−j3,9 )kV ∖ n

P₁ = 24, 05 MW

Q₁ = 14, 56 MVAr

Wzory wykorzystane do obliczenia przepływów mocy oraz strat mocy czynnej:

• prąd w linii
  ${}_{\mathbf{L}}\mathbf{=}\frac{{}_{\mathbf{p}}\mathbf{-}{}_{\mathbf{k}}}{\mathbf{R}_{\mathbf{L}}\mathbf{+ j}\mathbf{X}_{\mathbf{L}}}$
  

• prądy ładowania linii
  

$\mathbf{= j}\frac{\mathbf{B}_{\mathbf{L}}}{\mathbf{2}}\mathbf{\bullet}{}_{\mathbf{p}}$
$\mathbf{\ }\mathbf{= j}\frac{\mathbf{B}_{\mathbf{L}}}{\mathbf{2}}\mathbf{\bullet}_{\mathbf{k}}$

• prąd i moc wpływające do linii
  

_p=_L+

_p=_pp^*

• prąd i moc wypływające z linii
  
  _k=_L−
  _k=_kk^*
  

• straty mocy czynnej
  
   ΔP=P_p−P_k
  

Wartości mocy uzyskane w wyniku obliczeń w Matlabie:

_________________ _______________________

Moc czynną i bierną uzyskano na podstawie mocy pozornej.

• Linia 12

P_p=3, 69 MW

P_k=3, 62 MW

Q_p=17, 84 MVAr

Q_k=16, 78 MVAr ∖ n

P_p=27, 67 MW

P_k=27, 37 MW

Q_p=2, 23 MVAr

Q_k=1, 64 MVAr

P_p=46, 3 MW

P_k=45, 63 MW

Q_p=26, 36  MVAr

Q_k=25, 11MVAr

__________ _____________________________

Na podstawie tych wartości przedstawiono na poniższym rysunku rozpływ mocy w liniach .

Rozpływy mocy

Wartości strat mocy uzyskane w wyniku obliczeń w Matlabie

_________ _______________________________

Linia 1-2

ΔP₁₂ =  0, 072 MW

Linia 1-3
ΔP₁₃ = 0,  302MW

Linia 2-3
ΔP₂₃ = 0, 674 MW

__________ ______________________________