___ _____________________________________
Temat: Obliczanie rozpływów mocy u sieci rozdzielczej
___________ _____________________________
Bąk Wojciech
EN-DI-3
L
Rok akademicki 2013/2014
23.01.2014r.
Zadanie
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 (Pi, Qi) in all nodes of electric network. Note: sign „+“ determines power supply, sign „-“ determines power comsumption,
active and reactive power flows (Pij, Qij) 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 (95_AlFe) |
0,319 | 0,428 | 2,690 | 20 |
1-3 (120_AlFe) |
0,225 | 0,411 | 2,750 | 15 |
2-3 (150_AlFe) |
0,201 | 0,403 | 2,850 | 25 |
Rozwiązanie
Na podstawie poniższych wzorów obliczmy dla danych linii rezystancje ( R ) ,
reaktancje ( X ) , susceptancję ( B ) .
R = Rk • l [Ω]
X = Xk • l [Ω]
B = Bk • l [μS]
Wyniki przedstawiono w poniższej tabeli
Linia | R [Ω] | X [Ω] | B [µS] |
---|---|---|---|
1-2 (95_AlFe) |
6,38 | 8,56 | 53,8 |
1-3 (120_AlFe) |
3,375 | 6,165 | 41,25 |
2-3 (150_AlFe) |
5,025 | 10,075 | 71,25 |
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\}$$
Yii = ||
admitancje wzajemne gdy miedzy węzłami są linie
$= - \frac{1}{R_{\text{ij}} + j \bullet X_{\text{ij}}}$
Yij = ||
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
Wartości uzyskane w wyniku obliczeń w Matlabie:
|| = 120540 V = 121 kV
α22 = ( 0, 5638)≈0, 6
=( 120540+j1186,1 ) V = ( 121+j1,2)kV
|| = 118530 V = 119 kV
α33 = (−0,4090) ≈ −0, 4
=( 118530−j846) V = ( 119−j0,8)kV ∖ n
P1 = −4, 4 MW
Q1 = −6, 32 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=Pp−Pk
Wartości mocy uzyskane w wyniku obliczeń w Matlabie:
_________________ _______________________
Moc czynną i bierną uzyskano na podstawie mocy pozornej.
Linia 12
Pp=21, 03 MW
Pk=20, 82 MW
Qp=5, 7 MVAr
Qk=6, 2 MVAr ∖ n
Pp=16, 41 MW
Pk= 16, 35 MW
Qp=0, 13 MVAr
Qk=0, 33 MVAr
Pp=28, 97 MW
Pk=28, 64 MW
Qp=9, 30 MVAr
Qk=9, 67 MVAr
__________ _____________________________
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
ΔP12 = 0, 211 MW
Linia 1-3
ΔP13 = 0, 064 MW
Linia 2-3
ΔP23 = 0, 324 MW
__________ ______________________________