Chod sítě

Chod sítě (Load Flow), neboli výpočet ustáleného chodu třífázové souměrně zatížené elektrické sítě o harmonickém napětí, je spolu se zkratovými výpočty základní výpočetní úloha nad modelem sítě, daného topologií sítě a parametry jejích větví (rezistance, reaktance, konduktance a susceptance). Cílem výpočtu je nad modelem sítě na základě zadaných komplexních výkonových injekcí v uzlech sítě stanovit komplexní napětí v uzlech sítě a potažmo komplexní výkonové toky po větvích sítě.

Výpočet chodu sítě

Náhradní jednopólové schéma vedení (výše) a transformátoru (níže)

Chod kruhové sítě

Zkonstruujme výkonovou bilanci komplexního výkonu injektovaného v $i$-tém uzlu sítě o $n$-uzlech:

${𝐒}_i=\sqrt{3}{𝐔}_i{{𝐈}_i}^{*}={𝐔}_i\sum _j^{}{{𝐘}_{ij}}^{*}{{𝐔}_j}^{*}=P_i+i{Q}_{i}i,j=1,\dots ,n$

kde $𝐔$ resp. $𝐈$ jsou fázory příslušného napětí resp. proudu a $𝐘$ je příslušný prvek admitanční matice, tj.:

$f_i(\overrightarrow{U},\overrightarrow{\phi })=U_i\sum _j^{}{Y}_{\text{ij}}{U}_j\cos ({\phi }_i-{\phi }_j-{\alpha }_{\text{ij}})=P_i$

$f_{i+n}(\overrightarrow{U},\overrightarrow{\phi })=U_i\sum _j^{}{Y}_{\text{ij}}{U}_j\sin ({\phi }_i-{\phi }_j-{\alpha }_{\text{ij}})=Q_i$

kde $\phi$ resp. $\alpha$ je úhel napětí v příslušném uzlu resp. úhel admitance mezi příslušnými uzly v polárním tvaru.

Výpočet chodu sítě pak spočívá v řešení výše uvedené soustavy nelineárních rovnic uzlových komplexních výkonových bilancí pro velikosti a úhly napětí v uzlech např. Newtonovou iterační metodou podle věty o pevném bodě, zaručující existenci jistého okolí řešení s vlastností, že leží-li v něm počáteční aproximace řešení ${\overrightarrow{x}}_0$ tvořená jmenovitými napětími a nulovými úhly, algoritmus konverguje s přesností $\overrightarrow{\epsilon }$ k řešení $\overrightarrow{x}=[\overrightarrow{U},\overrightarrow{\phi }]$.

Funkce $f_i$ můžeme totiž aproximovat pomocí prvních dvou členů Taylorova rozvoje následující linearizací:

$f_i\left({\overrightarrow{x}}_{p+1}\right)=f_i\left({\overrightarrow{x}}_p\right)+d{f}_i\left({\overrightarrow{x}}_p\right)d{f}_i\left({\overrightarrow{x}}_p\right)={\mathrm{\nabla }f}_i\left({\overrightarrow{x}}_p\right)({\overrightarrow{x}}_{p+1}-{\overrightarrow{x}}_p)i=1,\dots ,2n$

kde $p$ je počítadlo iterací ($p=0,1,2,\cdots$) a $d{f}_i$ je totální diferenciál funkce $f_i$ v daném bodě představující tečnou nadrovinu k funkci $f_i$ v daném bodě a $\mathrm{\nabla }{f}_i$ je gradient funkce $f_i$ v daném bodě představující směr maximálního růstu funkce $f_i$ v daném bodě a pro vektorovou funkci $\overrightarrow{f}$ pak můžeme zapsat pro ${\overrightarrow{x}}_{p+1}$ dostatečně blízké řešení $\overrightarrow{x}$:

$\overrightarrow{f}\left({\overrightarrow{x}}_{p+1}\right)=\overrightarrow{f}\left({\overrightarrow{x}}_p\right)+𝐉\left({\overrightarrow{x}}_p\right)({\overrightarrow{x}}_{p+1}-{\overrightarrow{x}}_p)\le \overrightarrow{\epsilon }$

kde $𝐉$ je Jacobiho matice vyjadřující obecnou derivaci vektorové funkce $\overrightarrow{f}$ v daném bodě:

$𝐉=\left[\begin{array}{cc}\frac{\partial 𝐏}{\partial \mathit{\phi }}& \frac{\partial 𝐏}{\partial 𝐔}\\ \frac{\partial 𝐐}{\partial \mathit{\phi }}& \frac{\partial 𝐐}{\partial 𝐔}\end{array}\right]$

tj. kde pro ${\psi }_{ij}={\phi }_i-{\phi }_j-{\alpha }_{ij}$ a $i,j,k=1,\dots ,n$ platí:

$\frac{\partial 𝐏}{\partial \mathit{\phi }}=\left[\begin{array}{ccc}-U_1\sum _k^{}{a}_{1k}+U_1{a}_{11}& \cdots & U_1{a}_{1n}\\ ⋮& \ddots & ⋮\\ U_n{a}_{n1}& \cdots & -U_n\sum _k^{}{a}_{nk}+U_n{a}_{nn}\end{array}\right]$

$\frac{\partial 𝐐}{\partial \mathit{\phi }}=\left[\begin{array}{ccc}{U}_1\sum _k^{}{b}_{1k}-U_1{b}_{11}& \cdots & -U_1{b}_{1n}\\ ⋮& \ddots & ⋮\\ -U_n{b}_{n1}& \cdots & U_n\sum _k^{}{b}_{nk}-U_n{b}_{nn}\end{array}\right]$

$\frac{\partial 𝐏}{\partial 𝐔}=\left[\begin{array}{ccc}\sum _k^{}{b}_{1k}+c_{11}& \cdots & c_{1n}\\ ⋮& \ddots & ⋮\\ c_{n1}& \cdots & \sum _k^{}{b}_{nk}+c_{nn}\end{array}\right]$

$\frac{\partial 𝐐}{\partial 𝐔}=\left[\begin{array}{ccc}\sum _k^{}{a}_{1k}+d_{11}& \cdots & d_{1n}\\ ⋮& \ddots & ⋮\\ d_{n1}& \cdots & \sum _k^{}{a}_{nk}+d_{nn}\end{array}\right]$

kde: $a_{ij}=U_j{Y}_{ij}\sin {\psi }_{ij}{b}_{ij}=U_j{Y}_{ij}\cos {\psi }_{ij}{c}_{ij}=U_i{Y}_{ij}\cos {\psi }_{ij}{d}_{ij}=U_i{Y}_{ij}\sin {\psi }_{ij}$

tj. Jacobiho matice výše uvedené soustavy nelineárních rovnic je singulární, regularity dosáhneme vypuštěním úhlu napětí libovolně vybraného uzlu sítě ze seznamu proměnných spolu s vypuštěním bilanční rovnice činného výkonu vybraného uzlu ze seznamu rovnic, tj. vypuštěním odpovídajícího sloupce a řádku matice s tím, že vypuštěný úhel napětí, v tzv. referenčním uzlu, pak musíme zadat, např. jako nulový, k němuž se pak všechny další úhly napětí ostatních uzlů sítě budou vztahovat, a ve vybraném uzlu pak lze dopočtem vyrovnat bilanci činného výkonu v síti.

Pozn.: Závislost činného výkonu na velikosti napětí je slabší než na úhlu napětí a závislost jalového výkonu na úhlu napětí je slabší než na velikosti napětí, tj. za cenu ztráty přesnosti výpočtu lze mimodiagonální submatice Jacobiho matice považovat za nulové, tj. obecně nesymetrická Jacobiho matice přejde v matici symetrickou, čímž se výpočet chodu sítě (Decoupled Load Flow) urychlí. Za cenu další ztráty přesnosti výpočtu pomocí dalších aproximací lze výpočet chodu sítě (Fast Decoupled Load Flow) dále urychlovat.

Výkonové toky po větvích sítě pak určíme následovně:

$P_{ij}=A_{ij}{U_i}^2-U_i{p}_{ij}{U}_j(A_{ij}\cos {\alpha }_{ij}-B_{ij}\sin {\alpha }_{ij})$

$Q_{ij}=B_{ij}{U_i}^2-U_i{p}_{ij}{U}_j(B_{ij}\cos {\alpha }_{ij}+A_{ij}\sin {\alpha }_{ij})$

$A_{ij}=\frac{R_{ij}}{Z_{ij}^2}{B}_{ij}=\frac{X_{ij}}{Z_{ij}^2}{Z}_{ij}^2=R_{ij}^2+X_{ij}^2{\alpha }_{ij}={\phi }_i-{\phi }_j{p}_{ij}=\frac{U_i^{nom}}{U_j^{nom}}{U}_i^{nom}\ge U_j^{nom}$

a kde $P_{ij}$ resp. $Q_{ij}$ je činný resp. jalový výkon tekoucí z $i$-tého uzlu do $j$-tého uzlu po větvi o parametrech $R_{ij}$, $X_{ij}$ přepočtených k $i$-tému uzlu.

Chod paprskové sítě

Zkonstruujme výkonovou bilanci komplexního výkonu tekoucího $i$-tou větví sítě propojující uzly $i$ a $i+1$:

${𝐒}_i-{𝐒}_{i+1}=\sqrt{3}({𝐔}_i-{𝐔}_{i+1}){𝐈}_i^{*}=3{𝐙}_i{\left|\frac{{𝐒}_i^{*}}{\sqrt{3}{𝐔}_i^{*}}\right|}^2={𝐙}_i\frac{P_i^2+Q_i^2}{U_i^2}$

kde $𝐔$ resp. $𝐈$ jsou fázory příslušného napětí resp. proudu a $𝐙$ je impedance příslušné větve, tj.:

$f_i^P(P_i,Q_i,U_i^2)=P_i-R_i\frac{P_i^2+Q_i^2}{U_i^2}=P_{i+1}$

$f_i^Q(P_i,Q_i,U_i^2)=Q_i-X_i\frac{P_i^2+Q_i^2}{U_i^2}=Q_{i+1}$

$f_i^U(P_i,Q_i,U_i^2)=U_i^2-2(R_i{P}_i+X_i{Q}_i)+Z_i^2\frac{P_i^2+Q_i^2}{U_i^2}=U_{i+1}^2$

Výpočet chodu sítě pak spočívá v postupném řešení výše uvedených soustav nelineárních rovnic výkonových bilancí toků činných resp. jalových výkonů po větvích sítě spolu s napětím v jejich krajních uzlech, a to např. Newtonovou iterační metodou. Derivace funkcí $f_i$ můžeme aproximovat při zavedení poměrných jednotek ($\left|{𝐔}_i\right|\cong 1$) pro $R_i,X_i\ll 1$ následovně:

$\frac{\partial f_i^P}{\partial P_i}=1-2P_i\frac{R_i}{U_i^2}\cong 1\frac{\partial f_i^P}{\partial Q_i}=-2Q_i\frac{R_i}{U_i^2}\cong 0\frac{\partial f_i^P}{\partial U_i^2}=R_i\frac{S_i^2}{U_i^4}\cong \mathrm{\Delta }{P}_i$

$\frac{\partial f_i^Q}{\partial P_i}=-2P_i\frac{X_i}{U_i^2}\cong 0\frac{\partial f_i^Q}{\partial Q_i}=1-2Q_i\frac{X_i}{U_i^2}\cong 1\frac{\partial f_i^Q}{\partial U_i^2}=X_i\frac{S_i^2}{U_i^4}\cong \mathrm{\Delta }{Q}_i$

$\frac{\partial f_i^U}{\partial P_i}=-2{(R}_i-P_i\frac{Z_i^2}{U_i^2})\cong 0\frac{\partial f_i^U}{\partial Q_i}=-2(X_i-Q_i\frac{Z_i^2}{U_i^2})\cong 0\frac{\partial f_i^U}{\partial U_i^2}=1-Z_i^2\frac{S_i^2}{U_i^4}\cong 1$

tj. determinant Jacobiho matice výše uvedené soustavy nelineárních rovnic je přibližně roven jedné a soustava je tedy jednoznačně řešitelná.

Pozn.: Chod paprskové sítě sice zanedbává příčné parametry větví (menší přesnost výpočtu), zato probíhá výrazně rychleji než chod kruhové sítě.

Související články

• Elektrická síť
• Elektrická přenosová soustava
• Elektrická distribuční soustava
• Inteligentní sítě

Externí odkazy

• Výpočty chodu sítě