EXAMENSARBETE INOM ELEKTROTEKNIK, GRUNDNIVÅ STOCKHOLM 2014
KTH SKOLAN FÖR TEKNIK OCH HÄLSA www.kth.se
TRITA-STH-2014:44
Justerbar modell av
transmissionsledning för elkraftsöverföring
ALEXANDER SVENSSON MARCIAL
ANDREAS GATU
Justerbar modell av transmissionsledning för elkraftsöverföring
An adjustable model of a transmission line for power transmission
Alexander Svensson Marcial Andreas Gatu
Examensarbete inom Elektroteknik, Grundnivå, 15 hp
Handledare på KTH: Anna Josefsson Examinator: Thomas Lindh
TRITA-STH 2014:44
KTH
Skolan för Teknik och Hälsa
136 40 Handen, Sverige
Sammanfattning
Detta examensarbete har utförts på uppdrag av Terco. Tercos PST 2220 Transmission Line and Distribution Module fungerar som en fysisk modell av ett verkligt transmissions- eller distributionsnät där fem olika typer av nät med avseende på längd, spänningsnivå och skenbar effekt är möjliga. Idag finns ett behov av en modell där användaren kan ställa in dessa parametrar så att modellen mer precist kan spegla vilka egenskaper det specifika nätet har. Här undersöks hur längden och dess inverkan på en ledning kan varieras i en modell.
En presentation av hur transmissions- och distributionsnät fungerar och beskrivs teoretiskt lägger grunden till den modell och de två approximationer som kan beskriva ett helt näts egenskaper.
Då R, L och C komponenterna behöver kunna varieras för att fysiskt kunna realisera denna teoretiska modell undersöks vilka metoder som detta kan genomföras på. För detta undersöks två tillvägagångssätt, kaskadkopp- lad pi-modell och variabel aktiv-passiv reaktans.
Flera aspekter som utrymme, kostnad och variabilitet gör att varia- bel aktiv-passiv reaktans är att föredra. Dess funktion som en varierbar spänningskälla, uppbyggd av switchar styrda med reglerteknik och puls- breddsmodulering gör att komponenterna R, L och C och dess egenskaper och inverkan på en transmissionslinje kan åstadkommas. Resultatet är att de nödvändiga R,L,C komponenterna går att variera i storlek för att kunna ingå i en varierbar transmissionsledningsmodell.
Resultatet och målen säkerställs med simuleringar där variabel aktiv- passiv reaktans visas kunna vidareutvecklas och praktiskt testas för att mo- dellera transmissions- och distributionsnät med olika längd.
Nyckelord. Variabilitet, Inverterare, Impedans, DC-AC, Pulsbreddsmodu-
lering, Övertoner, Transmissionsledning, Spänningsfall, Reaktiv effekt.
Abstract
english This diploma work has been carried out on behalf of Terco. Tercos PST 2220 Transmission Line and Distribution Module works as a physical model of a real transmission and distribution grid where five different types of networks based on length, voltage and apparent effect are available. There is today a need of a model where the user self can adjust these parameters so that the model more precisely can reflect the characteristics that the specific grid has. Here it’s investigated how the length and its impact on a line can be varied in a model.
A presentation of how the transmission and distribution grid works and are described theoretically provides the basics for the different models that can describe a whole network and its properties.
Since the R, L and C components needs to be able to be varied to be able to physically realize this theoretical model, the different methods that this can be realized through are investigated. Two approaches are investigated, the cascaded pi-model and variable active-passive reactance (VAPAR).
A number of aspects like space, cost and variability makes the variable active-passive reactance the most suited solution. Its function as a variable voltage source, made out of an four switches, operated with control techno- logy and pulse width modulation, makes it possible to imitate R, L and Cs properties and effect on a transmission line. The result is that the necessary R,L,C components are made adjustable in order to be incorporated in a adjustable transmission lin model.
The result and the goal are verified with simulations where variable active-passive reactance is proved able for further development and practical tests to model transmission and distribution lines with different length.
Keywords. Variability, Inverter, Impedance, DC-AC, Pulse width modula-
tion, Harmonics, H-bridge, Transmission line, Voltage drop, Reactive effect.
Förkortningar
R - Resistans L - Induktans C - Kapacitans P - Aktiv effekt Q - Reaktiv effekt S - Skenbar effekt Z - Impedans X - Reaktans
VAr - Volt ampere reaktiv PWM - Pulsbreddsmodulering
IGBT - Insulated-gate bipolar transistor
FACTS - Flexible alternating current transmission system VAPAR - Variable active-passive reactance
CBCR - Converter based controllable reactance BVI - Bootstrap variable inductance
AVI - Active variable inductance
THD - Total harmonic distortion
Innehåll
Figurer 1
1 Inledning 5
1.1 Bakgrund . . . . 5
1.2 Problemformulering . . . . 5
1.3 Målsättning . . . . 6
1.4 Avgränsningar . . . . 7
2 Teori och bakgrund 9 2.1 Resistans . . . . 9
2.2 Induktans . . . . 10
2.3 Kapacitans . . . . 10
2.4 Modellering av transmissionledningar . . . . 11
2.4.1 Modell 1 . . . . 11
2.4.2 Modell 2 . . . . 13
2.4.3 Modell 3 . . . . 14
2.4.4 Modellernas noggrannhet vid olika ledningslängder . . 14
3 Metoder och tillämpningar 19 3.1 Kaskadkopplade pi-länkar . . . . 19
3.2 Varierbar impedans . . . . 21
3.2.1 Variabel spänningskälla för impedanskompensering . . 21
3.2.2 Converter based controllable reactance (CBCR) . . . . 22
3.2.3 Uppbyggnad av den variabla impedansen . . . . 23
3.2.4 Varierbar induktiv impedans . . . . 25
3.2.5 Varierbar kapacitiv impedans . . . . 29
4 Realisering av variabel reaktans 33
4.1 Varierbar spänningskälla . . . . 33
5 Resultat 37
5.1 Varierbar induktiv impedans . . . . 37 5.2 Simulering av variabel kapacitans . . . . 39 5.3 Simulering av variabel transmissionsledning . . . . 39
6 Analys och diskussion 43
7 Slutsatser 47
Källförteckning 49
Bilagor 50
A Induktor i Simulink 51
Figurer
1.1 PST2220: Transmissionsmodul . . . . 6
2.1 Placering av transmissionsledningar . . . . 10
2.2 Impedanser fördelade längs transmissionsledningen . . . . 12
2.3 Ekvivalent Pi-modell . . . . 13
2.4 Modell 2: Enfasig ekvivalent över en transmissionsledning . . . . 14
2.5 Modell 3: Enfasig ekvivalent över en transmissionsledning . . . . 15
2.6 Zsom funktion av längden för de tre modellerna. Övre figur visar beloppet, nedre argumentet (grader). . . . 16
2.7 U R vid olika längder och last. a) transmissionslängd från 5 upp till 100 km. b) en last som går från kapacitiv till induktiv . . . . 17
2.8 U R vid olika längder och last. a) transmissionslängd från 5 upp till 100 km. b) en last som går från kapacitiv till induktiv . . . . 17
3.1 k stycken pi-länkar kaskadkopplade . . . . 19
3.2 Visardiagram och krets med induktans . . . . 22
3.3 Visardiagram och krets med insatt spänningskälla U i . . . . 23
3.4 vaparkrets . . . . 24
3.5 Blockschema över kretsen för den variabla impedansen . . . . 25
3.6 Blockschema med PI-regulator för resistiv-induktiv impedans . . 26
3.7 Bodediagram över simulerad impedans jämfört med verklig . . . 29
3.8 Bodediagram över filtrets överföringsfunktion . . . . 29
3.9 Blockschema för kapacitans med PI-regulator . . . . 30
3.10 Bodediagram över den varierbara kapacitansen, med och utan filterkapacitans. . . . 32
4.1 Inverterare/konverterare . . . . 34
4.2 Unipolar PWM . . . . 35
Figurer
4.3 Switchad spänning . . . . 35
5.1 Jämförelse mellan verklig impedans och simulerad . . . . 38
5.2 Jämförelse mellan verklig kapacitans/km och simulerad . . . . . 40
5.3 Trefasig transmissionsledning . . . . 40
A.1 Effektelektroniken som utgör DC-AC konverteraren i den variab- la induktorn . . . . 52
A.2 Reglerkretsen för den variabla induktorn . . . . 53
A.3 Pulsbreddsmoduleringen . . . . 53
2
Kapitel 1
Inledning
1.1 Bakgrund
Terco är ett företag som utvecklar och tillverkar laborationsutrustning för skola och utbildning på gymnasial och eftergymnasialnivå. En av Tercos produkter är en elkraftssimulator som simulerar förloppet från kraftgenere- ringen till slutförbrukningen. Simulatorn är en fysisk modell där syftet är att studenter med inkoppling av en generator ska kunna mäta, avläsa och skapa sig en förståelse för hur ett elkraftsnät fungerar. Simulatorn går under namnet PST 2200 och är uppbyggd av fem stycken moduler som kan fungera fristående. En av modulerna i PST 2200 heter PST 2220 Transmission Line
& Distribution Module. Denna modul simulerar en transmissionsledning och är föremål för detta arbete.
1.2 Problemformulering
Ibland nöjer sig inte kunder med PST 2220 standardutförande, de önskar att den skall kunna simulera transmissionsledningar enligt egna specifikationer.
Detta innebär att Terco måste designa en ny modell, något som skulle kunna undvikas i en varierbar modell.
Terco har därför en önskan om att uppgradera modulen PST 2220 till en
simulator där användaren själv kan bestämma vilken längd transmissions-
ledningen skall ha och därmed reducera de befintliga fem transmissionsled-
ningar till en modell. En sådan modell skulle också ge användaren större
flexibilitet och anpassningsbarhet till olika experiment och det lokala elnä-
KAPITEL 1. INLEDNING
Figur 1.1: PST2220: Transmissionsmodul
tets utformning.
Det finns flera teoretiska modeller över en transmissionsledningar. De har alla gemensamt att transmissionledningen modelleras med impedanser.
Dessa impedanser varierar med ledningens längd. Detta innebär att, för att realisera en transmissionsledning av varierbar längd så måste dessa impedan- ser kunna varieras för märkspänningen och -strömmen 400V/2A trefas.
1.3 Målsättning
Målet med examensarbetet är att föreslå en konstruktion som kan simule- ra en transmissionsledning där transmissionsledningens längd kan varieras.
Den modell som föreslås i detta arbete skall kunna fungera i ett 400V/2A trefassystem.
Med en jämförelse av de matematiska modeller som beskriver elkraftsö- verföring, väljs den bäst lämpade modellen för att uppfylla modellens syfte.
Två stycken metoder för att realisera detta undersöks. Den första metoden bygger på att kaskadkoppla mindre segment av en transmissionsledningsmo- dell där varje segment består av fasta värden på impedanserna. Den andra metoden baseras på effektelektronik och reglerteknik i syfte att skapa vari-
4
1.4. AVGRÄNSNINGAR
abla impedanser.
Målet är att den föreslagna konstruktionen också skall kunna realiseras som ett system som är fysiskt mindre än det befintliga PST2220 systemet.
Lösningen skall vidare vara ekonomiskt försvarbar då det kommer till pro- duktion.
De som använder detta system använder extern mätutrusning så som multimeter och oscilloskop. Därför är en målsättning med den färdiga lös- ningen att den ström och spänning som mäts skall vara minimalt påverkad av störningar som lösningen kan medföra, ett oscilloskop skall, om än med störningar tydligt visa att det är en sinusvåg.
Det ska vara möjligt för Terco då de vill planera en uppgradering av PST2220 till en variabel modell att kunna utgå från den lösning som pre- senteras i detta arbete.
1.4 Avgränsningar
Den tilltänkta lösningen kommer att gälla en trefas kraftledning med sta- tionär växelspänning och symmetrisk impedanslast.
PST 2220 har transmissionsledningar av längderna 5,20,50 och 100 km. Där-
för begränsas detta arbete att ta fram en dynamisk modell där längden
är varierbar i intervallet 0 ≤ ≤ 100 [km]. Lösningen baseras på Terco’s
modell över en transmissionsledning för 230 kV.
Kapitel 2
Teori och bakgrund
Högspänd växelspänning transporteras längs luftledningar och markkablar.
Överföringen sker trefasigt vilket innebär att tre stycken fasledare fordras.
Detta kapitel kommer att beskriva tre kraftledningsparametrar (R,L,C), som sedan används för att ta fram en teoretisk modell över kraftledare vid sym- metriska förhållanden. Denna modell är sedan grunden för detta arbetes mål och syfte: att ge förslag på en justerbar fysisk modell.
2.1 Resistans
En ledare har en resistans R vilken kan beräknas som[10]
r T = ρ T
A [Ω/m] (2.1)
R T = r T · [Ω] (2.2)
I ekvation 2.1 är A [m 2 ] ledarens tvärsnittsarea och ρ T resistiviteten vid temperaturen T. I ekvation 2.2 är [m] ledarens längd.
På grund av denna resistans kommer en effekt att utvecklas i ledaren då en ström flyter genom den. Denna effekt definieras som en förlusteffekt och är lika med produkten av resistansen och strömkvadraten:
P f = RI 2 (2.3)
KAPITEL 2. TEORI OCH BAKGRUND
1
3
2
D 31
D 12
D 23
Figur 2.1: Placering av transmissionsledningar
2.2 Induktans
Induktans defineras som förhållandet mellan magnetflödet λ och strömmen I genom en ledare, d.v.s L = λ I . Om flera strömförande ledare ligger intill varandra påverkar ledarnas respektive magnetfält varandra och ger upphov till en en induktans som består av självinduktans och en ömsesidig induk- tans.
Induktansen för varje ledare i en trefasig transmissionsledning kan beräknas från nedan samband [10]
L = 2 · 10 −7 · ln
3√ D 12 D 23 D 31
D s (2.4)
Ur ekvation 2.4 är D 12 , D 23 , D 31 avstånden mellan fasledarna (figur 2.1). D s
är en storhet som uttrycker ledarens tvärsnittsradie.
Ekvation 2.4 förutsätter att trefasledarna är skruvade, det vill säga att de fysiskt byter position mellan varandra. Detta görs för att varje ledares in- duktans i genomsnitt skall vara lika varandra.
2.3 Kapacitans
Det är det elektriska fältet, spänningen, mellan ledarna och mellan ledar- na och marken som bestämmer kapacitansen. Kapacitansen kan generellt beräknas som [2].
8
2.4. MODELLERING AV TRANSMISSIONLEDNINGAR
C = ε · · a
d (2.5)
där
ε = dielektricitetskonstanten
= ledarlängd (m)
a = ledarens diameter (m)
d = avståndet mellan ledarna och ledare och jord (m)
Spänningen som uppstår mellan ledare modelleras med en kapacitans C o , benämnd ömsesidig kapacitans. Spänningen mellan respektive ledare och jord modelleras med kapacitansen C e , den så kallade egenkapacitansen.
I de modeller över transmissionsledningar som presenteras i avsnitten 2.4.1 och 2.4.2 används driftkapacitansen C d = C e + 3C o [2]
2.4 Modellering av transmissionledningar
En lednings resistans, kapacitans och induktans som beskrevs i de föregåen- de avsnitten är angivna per längdenhet. En transmissionsledning har dessa storheter jämt fördelade längs transmissionsledningens hela längd. Det finns tre stycken modeller beskrivna i litteraturen för att modellera transmis- sionledningar vid symmetrisk drift [10], [14] Dessa beskrivs nedan.
Tabell 2.1: Kapacitans och induktansvärden för den fysiska modellen
Spänningsnivå (kV) C (μF/km) L (mH/km) R (Ω/km
230 0.0250 0.31 0.04
2.4.1 Modell 1
Figur 2.2 illustrerar en ledning som är Δx m lång, vilket medför att im- pedansernas storlekar är zΔx [Ω] och yΔx [S]. Där z och y är lika med
z = r + jωl [Ω/m] (2.6)
KAPITEL 2. TEORI OCH BAKGRUND
y = jωc [S/m] (2.7)
−
+
U (x + Δx)
I(x + Δx)
yΔx
zΔx I(x)
+
− U (x)
Figur 2.2: Impedanser fördelade längs transmissionsledningen Genom att tillämpa Kirchoffs ström- och spänningslag erhålls samban- den
U (x + Δx) − U (x)
Δx = zI(x) (2.8)
I(x + Δx) − I(x)
Δx = y · U (x + Δx) (2.9)
Lösningarna till differentialekvationerna i (2.8) och (2.9) presenteras och sammanfattas i (2.10) där U (x) och I(x) är spänning respektive ström som funktion av ledningens längd [10].
U (x) I(x)
=
cosh(γx) Z c · sinh(γx)
Z 1
csinh(γx) cosh(γx)
U R
I R
(2.10) En modell som används för att modellera hur en transmissionsledning beter sig vid symmetrisk drift är den så kallade Pi-modellen (figur 2.3). Denna
10
2.4. MODELLERING AV TRANSMISSIONLEDNINGAR
modell består av en serieimpedans Z och två stycken shuntelement 1 2 Y pla- cerade på varsin sida av Z.
−
+
U S
I S
Y 2
I 1
Z I
Y 2
I 2
I R
+
− U R
Figur 2.3: Ekvivalent Pi-modell Ur kretsen i figur 2.3 kan följande samband ställas upp
U s
I s
=
1 + Y Z 2 Z Y 1 + Y Z 4 1 + Y Z 2
U R
I R
(2.11) Genom att utgå från (2.10) kan Y och Z i (2.11) bestämmas till
Z = Z c sinh(γx) (2.12)
Y
2 = tanh γx 2 Z c
(2.13) Med kännedom om ledarnas fysiska utförande och hur de är placerade i rummet i förhållande till varandra kan resistans, induktans och kapacitansen bestämmas. Genom att veta hur lång transmissionsledningen är kan man nu genom 2.11 beräkna spänning/ström i mottagaränden av ledningen samt hur effektflödet ser ut.
2.4.2 Modell 2
En annan förekommande modell är en approximation av den ovan nämnda
modellen. Här försummas effekten av kontinuerligt distribuerade längspara-
KAPITEL 2. TEORI OCH BAKGRUND
metrar och istället görs antagandet att de går att bunta ihop enligt nedan
Z = (R + jωL) · l [Ω] (2.14)
Y = (jωC) · l [S] (2.15)
Sambanden mellan inspänning/-ström och utspänning/-ström blir samma som i (2.11) men med Y och Z enligt (2.14) och (2.15)
−
+
U S
I S
1 2 c · l z · l
1 2 c · l I R
+
− U R
Figur 2.4: Modell 2: Enfasig ekvivalent över en transmissionsledning
2.4.3 Modell 3
I denna modell har ytterligare förenklingar gjorts, shuntelementen har för- summats, detta på grund av dess underordnade betydelse vid låga spän- ningsnivåer och korta transmissionsledningar.(figur 2.5).
2.4.4 Modellernas noggrannhet vid olika ledningslängder
I samtliga tre modeller är impedanserna funktioner av ledningens längd.
Nedanstående graf (figur 2.6) visar hur serieimpedansen Z varierar då trans- missionsledningens längd ökar. Här har ett intervall på 0 ≤ l ≤ 3000 km valts för att tydligt visa skillnaden.
Impedanserna i Figur 2.4 kan konverteras till impedanserna i Figur 2.2 genom att multiplicera dessa med korrigeringsfaktorer M 1 och M 2 som är lika med (2.16) och (2.17) [10]
12
2.4. MODELLERING AV TRANSMISSIONLEDNINGAR
−
+
U S
I S z · l I R
+
− U R
Figur 2.5: Modell 3: Enfasig ekvivalent över en transmissionsledning Tabell 2.2: Simulerade längsparametrar
R/km L/km C/km
12,73 mΩ 0.9337 mH 12,74 nF
M 1 = sinh(γ · l)
γ · l (2.16)
M 2 = tanh( γ·l 2 )
γ·l 2
(2.17) Om serieimpedansen i (figur 2.4) är Z blir den för (figur 2.2) Z · M 1 och för Y gäller Y · M 2 , dessa är komplexa storheter. Vid små ledningslängder blir Z = Z · M 1 ≈ Z och Y = Y · M 2 ≈ Y .
Den längsta längden som studeras i detta arbete är 100 km. Vid 100 km blir Z=1.2730 + j29.3331 Ω och Z’=1.2680 + j29.2758 Ω Vilket ger ett fel på 0.2% för att vid 5 km ha ett fel på 0.0068%. Detta visar att ju kortare en ledning är desto mer noggrann är approximationen av modell 2.
Vid lägre belastningar kommer spänningen för Figur 2.5 att avvika allt
mer. Vilken modell som används beror på vilka längsparametrar och vilken
längd som ska studeras. Av graferna att döma så kan modell 3, vid aktuella
längsparametrar, användas då ledningen är förhållandevis kort. Då det är
KAPITEL 2. TEORI OCH BAKGRUND
Figur 2.6: Z som funktion av längden för de tre modellerna. Övre figur visar beloppet, nedre argumentet (grader).
fråga om en längre ledning så måste modell 2 eller 3 tillgripas för att få ett mer nogrannt resultat. Det finns ingen absolut gräns för när modell 2 och 3 slutar att ge tillfredställande resultat, detta är något som får avgöras från fall till fall.
14
2.4. MODELLERING AV TRANSMISSIONLEDNINGAR
(a) Z last =110 Ω (b) |Z last | =110 med varierat argument Figur 2.7: U R vid olika längder och last. a) transmissionslängd från 5 upp till 100 km. b) en last som går från kapacitiv till induktiv
(a) Z last =1100 Ω (b) |Z last | =1100 med varierat argument
Figur 2.8: U R vid olika längder och last. a) transmissionslängd från 5 upp
till 100 km. b) en last som går från kapacitiv till induktiv
Kapitel 3
Metoder och tillämpningar
Detta kapitel behandlar de tekniker som studerats för att kunna realisera en modell av en transmissionsledning med variabel längd. Den modell som skall realiseras är modell 1 (Figur 2.2), detta då det är en bättre approxi- mation än modell 2 och 3 för aktuellt längdintervall. Modell 2 och 3 är en approximationer av modell 1.
3.1 Kaskadkopplade pi-länkar
Om en pi-länk (Figur 2.4) som är k km lång kopplas samman med k stycken identiska pi-länkar av längden k km så att den totala längden blir km så erhålls en mer förfinad modell än den i Figur 2.4 och blir en än mer nogrann approximation av den i Figur 2.2.
Pi-länk 1
Pi-länk 2
Pi-länk k
Figur 3.1: k stycken pi-länkar kaskadkopplade
Genom att utnyttja detta faktum är det möjligt att konstruera k stycken
segment som sedan kopplas in och ut beroende på önskad transmissions-
längd. Ju större k som väljs för systemet desto fler steg skapas i intervallet
KAPITEL 3. METODER OCH TILLÄMPNINGAR
0 ≤ ≤ 100 km, och därmed desto flera pi länkar krävs det. Detta betyder i praktikten 2k kapacitanser och k st resistanser och induktanser per fas.
Om den maximala längden som ska modelleras är = 100 km och k = 10 väljs blir antalet pi-länkar 10 st där varje pi-länk motsvarar 10 km. Detta resulterar i 60 st kapacitanser, 30 st resistanser och 30 st induktanser för ett trefassystem. I ett sådant system kan användaren ställa in en 0, 10, 20, 30, 40, 50, 60, 70, 80, 90 och 100 km lång transmissionsledning.
Vid simulering av denna krets fås värdena i tabell 3.1, och de använda RLC värdena är presenterade i tabell 3.2.
Tabell 3.1: Spänning och ström vid sändar och mottagarsida vid olika led- ningslängd
Länkar U S [V ] I S [A] U R [V ] I R [A]
0 325, 2651e −j0,0450 1, 9713e −j0,2632 325, 2651e −j0,0450 1, 9713e −j0,2632 1 325, 2651e −j0,0450 1, 9709e j3,4783 324, 5336e −j0,3869 1, 9669e −j0.6051 2 325, 2651e −j0.0451 1, 9804e j7,1856 323, 9295e −j0.7372 1, 9632e −j0.9554 3 325, 2651e −j0.0451 1, 9996e j10,8223 323, 4511e −j1.0963 1, 9603e −j1,3145 4 325.2651e −j0.0452 2, 0281e j14,3557 323, 0968e −j1,4643 1, 9581e −j1,6825 5 325, 2651e −j0.0452 2, 0656e j17,7578 322, 8654e −j1,8416 1, 9567e −j2,0598 6 325, 2651e −j0.0453 2, 1117e j21,0065 322, 7556e −j2,2284 1, 9561e −j2,4465 7 325, 2651e −j0.0453 2, 1658e j24,0857 322, 7666e −j2,6249 1, 9561e −j2,8431 8 325, 2651e −j0.0454 2, 2273e j26,5893 322, 8976e −j3,0315 1, 9569e −j3,2497 9 325, 2651e −j0.0454 2, 2959e j29,7001 323, 1479e −j3,4485 1, 9585e −j3,6667 10 325, 2651e −j0.0454 2, 3710e j32,2296 323, 5173e −j3,8763 1, 9607e −j4,0944
Tabell 3.2: Grundläggande längsparametrar
R/km L/km C/km
0,04 Ω 0,31 mH 0.1375 μF
Lasten som används i detta exempel består av 2oo Ω och 2 mH per fas.
18
3.2. VARIERBAR IMPEDANS
Tabell 3.3: Pi-länk parametrar för länkar om 10 km. Där C/länk fördelas i två lika stora delar vid sändarsida respektive mottagarsida.
R/länk L/länk C/länk 0,4 mΩ 3,1 mH 1,375 μF
3.2 Varierbar impedans
Då skenbar effekt överförs i transmissionsledningar ökar de aktiva effektför- lusterna med stigande överförd reaktiv effekt enligt
P f = R P 1 2 + Q 2 res
U 1 2 (3.1)
Vidare medför serieimpedansen att ett spänningsfall uppträder på ledningen och att maximal överförd effekt minskas. Detta har bidragit till att metoder utvecklats i syfte att motverka dessa problem. En enkel lösning för att mot- verka aktiva effektförluster är att ansluta en shuntkondensator som minskar den överförda reaktiva effekten. En liknande metod finns i att istället anslu- ta en kondensator i serie med ledningen och på så vis minskas ledningens reaktans och resulterar i ett lägre spänningsfall.
Det finns idag skäl att kunna variera reaktanser för att kunna möta de varierade fenomen som uppstår i elkraftssystem. Flexible alternating current transmission system (FACTS) är ett samlingsnamn för olika tekniker vars syfte är att dynamiskt kunna kontrollera olika transmissionsparametrar[1].
Gemensamt för dessa tekniker är att de nyttjar kraftelektronik, mer specifikt är det tekniker där switching används.
De lösningar som testats i detta arbete bygger på teknikerna beskrivna av FACTS. Istället för att minska en lednings reaktans skall istället dessa metoder användas till att skapa reaktans. Den metod som kommer användas i detta arbete för att skapa reaktans finns beskriven i [5][7][12][6][15] och be- nämns av sina upphovsmän till variable active-passive reactance (VAPAR).
I följande avsnitt beskrivs denna metod.
3.2.1 Variabel spänningskälla för impedanskompensering
Spänningen över induktorn i figur 2.4 kan tecknas som
KAPITEL 3. METODER OCH TILLÄMPNINGAR
U l = jωL L · I = jX l · I (3.2)
U l
U z
I
U i
(a) visardiagram
−
+
U i
L l
Z
+ U l −
+
− U z
(b) enfasig krets med induktor
Figur 3.2: Visardiagram och krets med induktans
Denna induktor kan ersättas med en spänningskälla med spänningen
U c = U l = jX l · I (3.3)
Ur ström och spänningsperspektiv är kretsen i figur 3.2b och 3.3b oföränd- rade. Detta kan även skådas ur figurerna 3.2a och 3.3a.
Om man i kretsen i figur 3.2b byter aktuell induktor mot en annan med ett annat värde kommer detta att återspeglas som en spänningsförändring i kretsen i figur 3.3b. Detta innebär att varierar man spänningen U c så kommer kretsen att se det som en förändring i induktans [8].
3.2.2 Converter based controllable reactance (CBCR)
Det har utvecklats tekniker vars syfte är att emulera reaktans, såväl ka- pacitiv som induktiv. Några av dessa är variable active- passive reactance (VAPAR) [5], boostrap variable inductance (BVI) [9] [3] och active variable inductance (AVI) [13]. De är alla uppbyggda på samma sätt men skiljer åt då det kommer till reglering av spänning och ström. Dessa går alla under namnet converter based controllable reactance (CBCR) [4]. VAPAR och AVI
20
3.2. VARIERBAR IMPEDANS
U c
U z
I
U i
(a) visardiagram
−
+
U i Z
+ U c −
+
− U z
(b) enfasig krets med induktor Figur 3.3: Visardiagram och krets med insatt spänningskälla U i
klarar att emulera en induktans och kapacitans medan BVI är strikt en vari- abel fiktiv induktans. CBCR har utvecklats i syfte att möta de problem som uppstår vid elkraftsöverföring. Det som gör de variabla impedanserna intres- santa, och då speciellt den induktiva impedansen, är att de har möjlighet att skapa en negativ induktans och kan på så sätt seriekoppla negativa in- duktanser på ledningar för att minska den effektiva serieimpedansen. CBCR är den andra metod som använts i detta arbete i syfte att skapa de vari- abla impedanserna som fodras för en justerbar transmissionsledning. Även om de olika teknikerna i CBCR är snarlika är det VAPAR som denna lös- ning bygger på, detta då den metoden även klarar att emulera en variabel kapacitans.
3.2.3 Uppbyggnad av den variabla impedansen
I avsnitt 3.2.1 beskrevs hur en spänningskälla kan härma beteendet för en impedans som kan vara både kapacitiv och induktiv genom att styra spän- ningen enligt (3.2). På samma metod bygger de ovan nämnda teknikerna i CBCR. I detta avsnitt skall grunden för uppbyggnaden av en variabel im- pedans presenteras i form av kretsscheman och ett ekvivalent blockschema.
Ur (figur 3.4) visas uppbyggnaden av den variabla impedansen. Ur figuren är V q den tidigare nämnda varierbara spänningskällan, L a ,L f ,C f och R f
ingår i ett lågpassfilter och Z(s) är en godtycklig yttre impedans och U s är
KAPITEL 3. METODER OCH TILLÄMPNINGAR
en yttre spänningskälla.
Lågpassfiltret är nödvändigt då den varierbara spänningskällan kommer att ge upphov till högfrekventa övertoner som ger störningar på spänning och ström i systemet, detta framkommer i avsnitt 4.1. Tillsammans utgör komponenterna, L a , L f , C f , R f och spänningskällan V q den variabla im- pedansen och kommer enligt ohms lag vara lika med Z i = U I
ii
.
U s
Z(s) I i L a
C f
I c
R f
L f I f
V q +
− U i
Figur 3.4: vaparkrets
Kretsens motsvarande blockschema kan ses ur (figur 3.5) där samban- den G 1 och G sys är framtagna med tillämpning av kirchoffs ström- och spänningslagar på kretsen i (figur 3.4). Detta är gjort med definitionen av strömmen I i (s) som utsignal och spänningarna V q (s) och U s (s) som insig- naler.
Den beskrivna kretsen med blockschema är generell och kan göras till både en kapacitiv impedans samt en induktiv impedans, detta beroende av hur spänningen V q styrs.
G 1 (s) = 1 + s 2 L f C f
1 + sR f C f
(3.4)
G sys (s) = Z(s) + s(L a + L f ) + (Z(s) + sL a )s 2 L f C f
1 + sR f C f
(3.5)
22
3.2. VARIERBAR IMPEDANS
1 G
sys(s)
G 1 (s)
I i
U s
V q e
Figur 3.5: Blockschema över kretsen för den variabla impedansen
3.2.4 Varierbar induktiv impedans
I detta avsnitt presenteras hur kretsen i föregående avsnitt skall regleras för att uppnå en induktiv impedans. Denna impedans skall sedan användas för att motsvara transmissionsledningen serieimpedans, vilken består av en resistans R och en induktans L.
Om spänningen U i (s) och strömmen I i (s) i kretsen (figur 3.4) följer sam- bandet
U i (s) = I i (s)(R + sL) (3.6)
Då kommer kretsen bete sig som en resistiv-induktiv belastning med im- pedansen Z iL = R + sL. Detta kan åstadkommas genom att med en PI- regulator reglera spänningen V q baserat på differensen av ärvärdet I i (s) och börvärdet I iref
Ärvärdet är den ström som går i kretsen och börvärdet är den önskade strömmen. Börvärdet fås genom att mäta spänningen U i (s) som sedan di- videras med R + sL, det vill säga genom att kasta om uttrycket i ekvation (3.6).
Enligt (figur 3.4) gäller sambandet U i (s) = U s − Z(s)I i (s) med vilken ett komplett blockschema kan konstrueras (figur 3.6).
Med blockschemat i (figur 3.6) kan en överföringsfunktion från utsignalen
I i till insignalen U s tecknas:
KAPITEL 3. METODER OCH TILLÄMPNINGAR
sL+R 1 G pi (s) G
sys1 (s)
Z(s)
G 1 (s)
U i I ref V c I i
− −
U s
U s
Figur 3.6: Blockschema med PI-regulator för resistiv-induktiv impedans
I i
U s
= (3.7)
= Z i s
2L
fC
f1+sR
fC
f+ Z i + G pi
Z i s
3L
aL
fC
f1+sR
fC
f+ sZ i (L a + L f ) + Z i G pi + Z(s) G pi + Z i + Z i s
2L
fC
f1+sR
fC
fG pi = K 1 + τ s
τ s (3.8)
Z i = sL + R (3.9)
Genom att utnyttja sambandet U i (s) = U s − Z(s)I i kan överförings- funktionen U I
ii
erhållas genom (3.7), som är den impedans som skall kunna varieras och ska enligt (3.6) skall vara lika med Z iL = R + sL
U i
I i = Z iL = (sL + R)H(s) (3.10) H(s) i uttryck (3.10) är en restfaktor som tillkommer på grund av filtret och regulatorn, denna är lika med
H(s) =
s
3L
aL
fC
f1+sR
fC
f+ s(L a + L f ) + G pi
(sL + R) 1 + 1+sR s
2L
fC
ff
C
f+ G pi
(3.11)
24
3.2. VARIERBAR IMPEDANS
För att Z iL skall kunna emulera impedansen sL+R så måste |H(ω)| = 1 och H(ω) = 0 kring området ω = 100π som är systemets frekvens.
För att underlätta analysen av H(s) kommer filtret att approximeras med den filtrering som induktansen L e = L a + L f gör. Denna approximation är giltig fram till filtrets resonansfrekvens. Denna frekvens bestäms utifrån filtrets parametrar, vilka är valda så resonansfrekvensen hamnar mycket högre än 50 Hz. Vidare så underlättar denna approximation konfigurationen av PI-regulatorn.
Med nämnd approximation blir nu (3.10) lika med:
Z iL = (sL + R) sL e + K 1+τs τ s
(sL + R) + K 1+τs τ s = (sL + R) · H (s) (3.12) Egenskaperna för H (s) undersöks här i ett bodediagram vars amplitudkurva tecknas
|H (ω)| dB = 20 log 10
1 + ω 4
2tn
(d 2 t − 1 2 )ω 2 + ω ω
44tn
1/2
1 + ω 4
2nn
(d 2 n − 1 2 )ω 2 + ω ω
44nn
1/2 (3.13)
I (3.13) kan det ses hur |H(ω)| dB är lika med 0 vid låga frekvenser.
Om täljaruttrycket studeras så är denna lika med enhetsförstärkning (=0 dB) vid alla frekvenser fram till f c1 = 2π 1 ω tn , vid denna frekvens börjar täljaren att öka med +40 dB/dekad. Det motsatta gäller nämnaruttrycket som istället avtar med -40 dB/dekad vid f c2 = 2π 1 ω tn . Detta innebär att beroende på hur valet av de två frekvenserna görs så kommer |H(ω)| dB att sluta på ett värde större eller mindre än noll. I detta arbete är frekvenserna valda så att f c2 < f c1 . Då kommer |H(ω)| dB gå mot en negativ förstärkning med ökande frekvens. Översatt till värden på induktanser innebär detta att L e ≤ L vilket framgår av följande två ekvationer:
ω tn =
K
L e τ (3.14)
ω nn =
K
Lτ (3.15)
KAPITEL 3. METODER OCH TILLÄMPNINGAR
I detta arbete är som tidigare nämnt 1.55 mH ≤ L ≤ 31 mH vilket leder till att valet av L e väljs till 1.55 mH.
Vidare så måste f c2 > 50 Hz, detta sätter följande restriktion på valet av K och τ
K
τ > (100π) 2 · L max ≈ 3060 (3.16) Konfigurationen av PI-regulatorn är gjord i Simulink med den inbyggda funktionen PID-tune. De värden på K och τ som erhölls där var K=72 och τ = 1/812.
Approximationen av filtret till L e gäller som tidigare nämnt fram till filt- rets resonansfrekvens, för vilken filtret böjer av och får en högre dämpning.
Filtrets resonansfrekvens är f resonans och är lika med (3.17).
ω resonans =
L a + L f
L a L f C f
(3.17)
Genom olika fördelningar av värdet på L e på L a och L f går det att tillsammans med kapacitansen C f att styra var filtrets resonansfrekvens skall hamna. De frekvenskomponenter som behöver filtreras bort i detta arbete ligger kring 40 kHz, vilket framgår i senare kapitel då implementationen av V q beskrivs.
Som störst blir L a · L f då värdet på L e fördelas lika mellan dem. Detta innebär att kondensatorstorleken kan göras minimal för en given resonans- frekvens. Vid placeringen av ω resonans måste denna ligga mellan systemets frekvens (ω ≈ 314 rad/s) och mindre än de frekvenser som filtret skall dämpa. Den frekvens som är vald i denna tillämpning och som uppfyller nämnt villkor är ω resonans = 1000rad/s vilket innebär att kondensatorn kan beräknas till C f =2500μF.
Med storheterna L a , L f och C f kan värdet på resistansen R d bestämmas.
Resistorn är till för att dämpa den höga förstärkning som inträffar vid reso- nansfrekvensen (figur 3.8). Värdet på resistansen är bestämt till R f = 0.5Ω baserat på jämförelse av bodediagram för olika värden på resistansen.
26
KAPITEL 3. METODER OCH TILLÄMPNINGAR
U i = 1
sC I i (3.18)
Skillnaden här jämfört med den induktiva impedansen är att det är ström som integreras istället för spänning. Från blockschemat (figur 3.9) kan sy- stemets totala överföringsfunktion tecknas som
I i
U s
= sC 1 + K 1+τs τ s + 1+sR s
2L
fC
ff
C
fs 2 C(L a + L f ) + sC s 1+sR
3L
fL
faC C
ff+ K 1+τs τ s + Z(s) · sC
1 + 1+sR s
2L
fC
ff
C
f+ K 1+τs τ s
(3.19) Från denna överföringsfunktion kan som tidigare impedansen tas fram till
Z c (s) = 1
sC H c (s) (3.20)
där H c (s) är restfaktorn som är lika med
H c (s) = s 2 C(L a + L f ) + sC s 1+sR
3L
fC
ffC L
fa+ K 1+τs τ s 1 + K 1+τs τ s + 1+sR s
2L
fC
ff
C
f(3.21)
G P I G sys
Z(s) sC 1
U i I i
U iref
U s
Figur 3.9: Blockschema för kapacitans med PI-regulator
28
3.2. VARIERBAR IMPEDANS
Som i fallet för den varierbara induktiva impedansen kommer analysen av H c (s) göras med filterapproximationen där kapacitansen försummas, detta leder till ett förenklat uttryck av den kapacitiva impedansen (3.22), denna förenkling görs i syfte att undersöka parametrarna K, L e och τ inverkan på hur väl (3.20) efterliknar en kapacitans. Genom att låta H c (s) övergå till H c (ω) så blir amplitudkurvan för H c (ω) lika med (3.23)
1 sC
s 2 CL a + K 1+τs τ s 1 + K 1+τs τ s = 1
sC H c (s) (3.22)
|H c (ω)| dB = 20 log 10
1 + ω 2 τ 2 (1 − ω 2 T 2 ) 2 1/2
1 + ω 2 τ 2 (1 + K 1 ) 2
1/2 (3.23)
T =
CL e
K (3.24)
Täljaruttrycket har två stycken brytpunkter ω t1 = 1 τ och ω resonans = T 1 . Nämnaruttrycket har en brytpunkt vid ω n = τ (1+ 1
1k