1 Uvod V hidroelektrarnah poganjajo generatorje vodne turbine. Kadar določena turbina, ki ni opremljena z zapornim organom, ne obratuje, pride do večjega ali manjšega puščanja vode skozi vodilne lopatice. To se dogaja pri Francisovi, Kaplanovi in drugih nadtlačnih turbinah. Ker puščanje vode pomeni izgubo energije, želimo poznati količino vode, ki uhaja skozi zaprte vodilne lopatice. Na opisani način lahko ovrednotimo izgubljeno energijo, poleg tega pa lahko primerjamo puščanje pred določenimi posegi na vodilnih lopaticah turbine in po njih. Nekatere hidroelektrarne so opremljene z ultrazvočnimi merilci pretoka, s katerimi bi načeloma bilo mogoče meriti pretok skozi turbino, ko so vodilne Prejet 21. april, 2010 Odobren 25. julij, 2010 240 Leban, Volk, Mihalič lopatice zaprte. Problem takega merjenja je, da je pretok puščanja majhen v primerjavi z nazivnim pretokom ( na splošno manj kot 1 %). Pri tako majhnih pretokih so tovrstne merilne naprave zelo nenatančne. Pretok skozi zaprte vodilne lopatice pa je mogoče določiti tudi s pomočjo barometra, če je oblika cevovoda med turbino in predturbinsko zapornico matematično opisljiva. Ker metode, ki bi določala puščanje na takšen način, ni bilo mogoče izslediti v literaturi, smo izdelali princip, ki ga opisujemo v nadaljevanju. 2 Osnovni princip Osnovni princip temelji na predpostavki, da je vsota preseka špranje, skozi katero pušča voda, med preizkusom konstantna in da velja [1]: pkQ ⋅=vodilnika , (1) pri tem sta: k…konstanta in p…razlika tlaka pred špranjo in za njo. 3 Metoda Pretok vode skozi turbino pri zaprtih vodilnih lopaticah določimo tako, da zapremo predturbinsko zapornico. Ob predpostavki, da je puščanje predturbinske zapornice manjše od puščanja turbine, se začne nivo vode v delu cevovoda pod predturbinsko zapornico nižati. Predpostavimo najprej, da predturbinska zapornica idealno tesni. Nižanje nivoja vode v delu cevovoda med predturbinsko zapornico in turbino je tako odvisno le od puščanja turbine. Ob nižanju nivoja vode ne sme nastati podtlak, ker bi bile meritve sicer napačne. Nastanek podtlaka preprečujejo prezračevalne cevi, ki se sicer uporabljajo za izpust zraka ob napolnitvi praznega cevovoda prek "bypassa" [2]. S pomočjo barometra merimo časovno odvisnost spremembe višine gladine vode. Ob poznavanju oblike cevovoda lahko tako določimo spremembo volumna v odvisnosti od časa in padca, s čimer dobimo tudi pretok skozi turbino v odvisnosti od padca. Princip merjenja prikazuje slika 1. Pretok skozi turbino v določenem trenutku je enak [3]: t V dt dV tQ ∆ ∆ ≅=)( , (2) pri čemer so: Q(t)…pretok skozi turbino, ∆V…sprememba volumna vode in ∆t…opazovano časovno obdobje. H sp.vode H vode1; t1 H vode2; t2 V habs. Predturbinska zapornica Vodilne lopatice Slika 1: Določitev spremembe volumna ∆V določimo iz spremembe gladine vode ∆H Fig. 1: Determining the volume change ∆V from a water level change ∆H. Približek t V dt dV ∆ ∆ ≅ velja tem bolj, čim krajše je opazovano časovno obdobje. Tako se določi pretok skozi zaprte vodilne lopatice turbine pri določeni gladini vode, ki je pod nivojem predturbinske zapornice, in ob določeni gladini spodnje vode. Iz tega lahko izračunamo namišljeni presek špranje zaprtih vodilnih lopatic, ki ga označimo z S [4]: .2 abshgt V v Q S ⋅⋅⋅∆ ∆ == , (3) kjer so: S…namišljeni presek špranje zaprtih vodilnih lopatic, Q…pretok skozi turbino, v…hitrost vode, g… gravitacijski pospešek, habs.…absolutni padec vode, ∆V… sprememba volumna vode in ∆t… opazovano časovno obdobje. Absolutni padec vode je enak povprečju nivoja vode v delu cevovoda v časovnem intervalu merjenja, zmanjšan za nivo spodnje vode: vodesp vodevode abs H HH h . 21 . 2 − + = , (4) pri čemer so: habs.… absolutni padec vode, Hvode1…višina vode na začetku meritve, Hvode2…višina vode na koncu meritve in Hsp.vode…višina spodnje vode. Napako zaradi predpostavke: t V dt dV ∆ ∆ = , (5) zmanjšamo tako, da se meritev spremembe volumna opravi v več točkah. Za končno vrednost namišljenega Določanje puščanja vodnih turbin in predturbinskih zapornic v hidroelektrarni Doblar 2 241 preseka S vzamemo povprečno dobljeno vrednost. Ko je namišljeni presek znan, se lahko preračuna puščanje vodilnih lopatic na dejanski nivo zgornje vode, saj nas zanima puščanje pri dvignjeni predturbinski zapornici (slika 2). Preračun temelji na enačbi: )(2 ..vodilnika vodespvodezg HHgSQ −⋅⋅⋅= , (6) kjer so: Qvodilnika…pretok skozi vodilne lopatice, S… namišljeni presek špranje zaprtih vodilnih lopatic, g… gravitacijski pospešek, Hzg.vode…višina zgornje vode in Hsp.vode…višina spodnje vode. H sp.vode H zg.vode Slika 2: Hidro agregat z dvignjeno predturbinsko zapornico Fig. 2: Hydro turbine with inlet water gate lifted 4 Puščanje predturbinske zapornice V praksi se je izkazalo, da puščanje predturbinske zapornice ni zanemarljivo. To pomeni, da izmerjena sprememba volumna vode v delu cevovoda ne izraža neposredno puščanja zaprtih vodilnih lopatic turbine, saj v opazovanem času določena količina vode priteče skozi predturbinsko zapornico. Količina vode, ki odteče skozi turbino v opazovanem času, je torej enaka vsoti spremembe izmerjenega volumna in količini vode, ki priteče v del cevovoda skozi predturbinsko zapornico. Da lahko to določimo, je treba poznati puščanje predturbinske zapornice. 4.1 Določitev puščanja predturbinske zapornice Puščanje predturbinske zapornice določimo tako, da pustimo zapornico zaprto, dokler se padanje vode v cevovodu ne ustavi. Nivo vode v delu cevovoda ne pade na nivo spodnje vode, saj pri določeni (majhni) višinski razliki zapornica pušča enako kot vodilne lopatice (slika 3). V stanju, ko gladina vode ne upada več (stacionarno stanje), velja, da je puščanje zapornice enako puščanju vodilnih lopatic turbine. H sp.vode_v_stac._st. H zg.vode_v_stac._st. hv_stac._st. Slika 3: Gladina vode pred vodilnimi lopaticami, po umiritvi upadanja. Fig. 3: Water level in front of guide vanes in the stationary state. Če je puščanje predturbinske zapornice enako puščanju vodilnih lopatic v stacionarnem stanju, potem velja: ,2 )(2 ._._ ._.__.._.__. ustac.stanj vvodilnikazapornice ststacv ststacvvodespststacvvodezg hgS HHgS QQ ⋅⋅⋅= =−⋅⋅⋅= == (7) pri čemer so: Qzapornice…pretok skozi predturbinsko zapornico, Qvodilnika v stac.stanju…pretok skozi vodilne lopatice v stacionarnem stanju, Hzg.vode_v_stac._st.…višina vode v delu cevovoda v stacionarnem stanju, Hsp.vode_v_stac._st.…višina spodnje vode v stacionarnem stanju in hv_stac._st.…absolutni vodni padec v stacionarnem stanju. Izmerjeni spremembi volumna (poglavje 3) je treba dodati še količino vode, ki se nateče skozi predturbinsko zapornico. Dobimo popravljeno vrednost spremembe volumna: S zaporniceV V Q t∆ = ∆ + ⋅∆ . (8) S pomočjo te enačbe lahko izračunamo namišljeni presek S: . .2 2 zaporniceS abs abs V Q tV S t g h t g h ∆ + ⋅∆∆ = = ∆ ⋅ ⋅ ⋅ ∆ ⋅ ⋅ ⋅ . (9) Tako dobimo dve enačbi (enačbi 9 in 7) z neznanima vrednostma S in Qzapornice. Če v enačbo 9 vstavimo enačbo 7, dobimo enačbo, s pomočjo katere lahko izračunamo namišljeni presek S, ki puščanje 242 Leban, Volk, Mihalič predturbinske zapornice upošteva, vendar vrednosti pretoka puščanja le-te ni treba poznati, če poznamo višinsko razliko med nivojem vode pred turbino in za njo ob spuščenih zapornicah in zaprtem vodilniku v stacionarnem stanju (po umiritvi padanja gladine vode v opazovanem odseku cevovoda). ( )._._. 22 ststacvabs hghgt V S ⋅⋅−⋅⋅⋅∆ ∆ = . (10) Seveda se lahko (če nas to zanima) izračuna tudi puščanje predturbinske zapornice s pomočjo enačbe 7, ko enkrat poznamo vrednost namišljenega preseka S. 5 Izračun spremembe volumna vode za cevovod, kakršen je v HE Doblar 2 Kot je bilo že omenjeno, je opisani postopek za določitev puščanja predturbinske zapornice in vodilnih lopatic turbine mogoče izvesti samo če je oblika cevovoda med predturbinsko zapornico in turbino znane oblike. 16 .1 6 17 .8 0 5.50 Slika 4: Oblika dela cevovoda v HE Doblar 2 Fig. 4: The form of the penstock section in HPP Doblar 2 Z drugimi besedami, vedeti moramo, kakšno spremembo volumna vode pomeni določena sprememba višine gladine vode. To je bilo izračunano za obliko odseka cevovoda v HE Doblar 2. Obliko cevovoda med predturbinsko zapornico in turbino v HE Doblar 2 bi lahko poimenovali "zvita" cev. Obliko cevovoda prikazuje slika 4. Del tega cevovoda je enak delu geometrijskega telesa, ki se imenuje torus. Presek torusa popisuje Cassinijeva krivulja [5], opisana s spodnjo enačbo (slika 5): ( ) ( )222222222 4 cxryxcar +⋅⋅=+++− , (11) pri tem so: a…polmer cevi, r…radij torusa in c…oddaljenost ravnine, ki preseka torus, od središča torusa [6]. Slika 5: Geometrija za izračun preseka torusa Fig. 5: Determination of a torus-cross-section area. Iz Cassinijeve enačbe izrazimo koordinato y in dobimo: ( )2222222 4 cxrxcary +⋅⋅+−−+−= . (12) Presek torusa se izračuna kot integral: ∫⋅= zb dxyS nb T 2 , (13) oziroma ( ) n 2 2 2 2 2 2 2 T b 2 4 zb S r a c x r x c dx= − + − − + ⋅ +∫ , (14) pri čemer sta meji integrala enaki: 22 hrhb zz −⋅⋅= in (15) 2)2()2(2 ahrahb nn ⋅−−⋅⋅−⋅= ; h > 2a oz. 0=nb ; h < 2a. (16) Pri tem so: arrz += , (17) arrn −= in (18) harc −+= . (19) Volumen v delu cevovoda nad gladino vode dobimo kot integral presekov torusa: T T 0 h V S dh= ∫ (20) in končno: Določanje puščanja vodnih turbin in predturbinskih zapornic v hidroelektrarni Doblar 2 243 ( ) dhdxcxrxcarV h bz 42 0 b 2222222 T n ∫ ∫ +⋅⋅+−−+−= . (21) Ker enačba 14 ni analitično rešljiva, smo volumen za opisano obliko dela cevovoda izračunali numerično, s pomočjo programa Matlab. Potek izračuna prikazuje spodnja koda: korak=0.001 V=0 for h=0:korak:16.16 c=r+a-h bz=sqrt(2*h*rz-h^2) if h>2*a bn=sqrt(2*(h-2*a)*rn-(h-2*a)^2) else bn=0 end y=0 S=0 for x=bn:korak:bz y=sqrt(-r^2+a^2-c^2-x^2+sqrt((4*r^2)*(x^2+c^2))) S=S+y*korak*2 end V=V+S*korak end V [ m 3 ] Slika 6: Odvisnost volumna od višine gladine vode za cevovod, kakršen je v HE Doblar 2 Fig. 6: Correlation between the water level depth and volume, calculated according to the penstock form of HPP Doblar 2. Za dimenzije, prikazane na sliki 4, smo dobili odvisnost volumna od oddaljenosti gladine vode od vrha predturbinske zapornice (h), kakor je prikazano na sliki 6. 6 Želeno maksimalno odstopanje meritev Če bi želeli določiti natančnost meritve, bi morali imeti rezultate referenčnih meritev, ki pa jih za sprejemljive stroške skoraj ni mogoče izvesti. Želimo si seveda čim natančnejše rezultate. Natančnost rezultatov je mogoče v določeni meri oceniti tudi na podlagi stresanja rezultatov. Oceniti želimo tudi podatek o puščanju predturbinske zapornice, saj je to merilo za morebitno menjavo tesnil ipd. Na točnost izračunane vrednosti vplivajo natančnost naprav za merjenje tlaka, natančnost dimenzij in oblike cevovoda ter natančnost zagotavljanja predpostavk med merjenjem (konstantnost gladine gornje vode - pomembne za izračun puščanja predturbinskih zapornic - in konstantnost gladine spodnje vode - pomembne za izračun puščanja vodilnika). Čeprav smo izpeljali odvisnosti spremembe volumna glede na spremembo nivoja vode za celotno cev, se je treba zavedati, da v naravi cevi niso gladko zaprte na gornjem in spodnjem koncu. Spodnji konec ni problematičen, ker tako ali tako leži pod nivojem spodnje vode in za meritve ne pride v poštev. Zgornji konec cevi pa se nadaljuje v zaporne organe, ki ne odrežejo vodnega stolpca idealno. Zato je smiselno nekoliko počakati in meritve izvesti šele takrat, ko se nivo vode v cevi spusti pod presek vhodne odprtine cevi. 7 Rezultati Opisano metodo smo uporabili na turbini v HE Doblar 2. Spremembe višine gladine vode smo merili s pomočjo barometra. Najprej smo pri odprti predturbinski zapornici odčitali tlak barometra. Iz natančnih meritev nadmorske višine gladine vode pripadajoče akumulacije smo "umerili" barometer. Tako smo preverili odvisnost med tlakom barometra in nadmorsko višino gladine vode. V naslednji fazi smo predturbinsko zapornico zaprli in približno eno uro odčitavali tlak barometra. Pri tem smo opravili približno 20 meritev. Nato smo po približno 24 urah ponovno odčitali tlak. Ta meritev nam je pokazala tlak v stacionarnem stanju, ki ga potrebujemo za eliminacijo vpliva puščanja predturbinske zapornice. Pri tem je bilo treba paziti, da v tem času nismo izvajali večjih sprememb obratovanja agregatov vzporedne elektrarne Doblar 1. To namreč vpliva na spremembe nivoja spodnje vode, kar odčitek tlaka v stacionarnem stanju popači (pravzaprav pri spremembah spodnje vode odmerek ni v stacionarnem stanju). Za vsako meritev sprememb smo izračunali puščanje vodilnih lopatic in puščanje predturbinske zapornice. Vrednosti puščanja vodilnih lopatic za različne meritve agregata Doblar 2 so se gibale med 0,40 m3/s in 0,46 m3/s s povprečjem 0,43 m3/s. Vrednosti puščanja predturbinske zapornice pa so se gibale med 0,18 m3/s in 0,21 m3/s. Povprečna vrednost puščanja predturbinske zapornice je bila 0,19 m3/s. Standardna deviacija je bila 4,6 % povprečne vrednosti. 7.1 Ovrednotenje izgubljene energije Poleg poznavanja puščanja vodilnih lopatic turbine nas zanima tudi ekvivalentna količina izgube energije. 244 Leban, Volk, Mihalič Povprečno letno število ur, ko je agregat v HE Doblar 2 v pripravljenosti in ne obratuje, je približno 3500 ur. Če upoštevamo izračunano puščanje vodilnih lopatic tega agregata, ki znaša približno 0,43 m3/s (0.41 % nazivnega pretoka), ter specifično porabo vode, ki je okrog 2,6 m3/s/MW, dobimo letno izgubljeno energijo, ki znaša skoraj 600 MWh. To je ekvivalentno izgubi približno 45.000 evrov na leto. Za primerjavo velja omeniti, da je bilo v letu 2007 v vseh sončnih elektrarnah v Sloveniji proizvedenih 600 MWh električne energije [7], kar je enako izgubam puščanja agregata v HE Doblar 2. Seveda pa tu niti ni upoštevano, da gre pri hidroagregatu za napovedljivo energijo, ki jo lahko izkoriščamo takrat, ko je to smiselno. Vrednost take energije pa je neprimerno večja od energije bolj ali manj nenapovedljive sončne energije, ki je ne moremo izkoriščati po želji. 8 Sklep Z opisano metodo je mogoče določiti puščanje vodilnih lopatic vodnih turbin tako rekoč brez investiranja v merilno opremo. Potrebujemo le razmeroma natančen barometer, ki je na voljo v skoraj vsaki hidroelektrarni, poleg tega pa moramo poznati geometrijo pretočnih organov med predturbinskimi zapornicami in vodilnikom. Meritve, s katerimi se določi pretok vode skozi turbino, pa stanejo tudi do 70.000 evrov [8]. Slabost metode je ta, da omogoča izvedbo meritev le v hidroelektrarnah, ki so opremljene s predturbinsko zapornico. Sama meritev je razmeroma hitra. Če želimo upoštevati tudi puščanje predturbinskih zapornic, je treba doseči ravnovesno stanje med puščanjem predturbinske zapornice in vodilnika, kar pa praviloma traja več ur. V tem času je merjeni agregat nerazpoložljiv, hkrati pa se ne sme spreminjati nivo spodnje vode, če želimo doseči zadostno natančnost. "Stranski produkt" tovrstnih meritev je določitev puščanja predturbinske zapornice. Opisana metoda bi lahko bila v pomoč tudi pri splošnih meritvah pretoka s pomočjo Gibsonove metode. Pri merjenju pretoka z omenjeno metodo je treba določiti pretok puščanja turbine, saj ta vpliva na končno vrednost pretoka, dobljenega s pomočjo Gibsonove metode [9, 10]. Meritve smo izvedli na novejšem agregatu hidroelektrarne Doblar 2. Uporabili smo dokaj natančen barometer, poleg tega pa so tudi dimenzije pretočnih naprav dokaj natančno definirane. Standardna deviacija je znašala približno 4,6%, kar je več kot zadovoljiva natančnost, saj je naš cilj oceniti puščanje, in ne določitev natančne vrednosti. Tako smo z že obstoječo opremo dobili želeno informacijo o puščanju vodilnih lopatic in predturbinske zapornice. Omenimo tudi, da lahko označimo nabor nam dostopne literature, ki bi to ali tej podobno metodo merjenja puščanja predturbinskih lopatic opisovala kot neobstoječega. 9 Viri [1] Frank M. White: Fluid Mechanics – 4th edition, McGraw Hill, 1998 [2] J. Giesecke, E. Mosonyi: Wasserkraftanlagen: Planung, Bau und Betrieb, Springer, 2003 [3] C. T. Crowe, D. F. Elger, B. C. Williams, J. A. Roberson: Engineering Fluid Mechanics – Ninth Edition, Wiley, 2009 [4] F. Kreith: Mechanical Engineering Handbook, CRC Press, 1999 [5] M. R. Spiegel, S. Lipschutz, J. Liu: Mathematical Handbook of Formulas and Tables – Third Edition, McGraw Hill, 2009 [6] Weisstein, Eric W. "Spiric Section." From MathWorld - A Wolfram Web Resource. http://mathworld.wolfram. com/SpiricSection.html, 2009 [7] Institut Jožef Stefan: Letni energetski pregled za leto 2007, Ljubljana, 2009 [8] EPRI: Flow Measurement at Hydro Facilities: Achieving Efficiency, Compliance, and Optimal Operation, 2002 [9] A. Adamkowski, W. Janicki, J. Kubiak, G. Urquiza, F. Sierra, J. M. Fernandez D.: Water Turbine Efficiency Measurements Using the Gibson Method Based on Special Instrumentation Installed Inside Pipelines, 6th International Conference on Innovation in Hydraulic Efficiency Measurements, 2006 [10] E. da Costa Bortoni: New Developments in Gibson’s Method for Flow Measurement in Hydro Power Plants, Flow Measurement and Instrumentation 19, Elsevier, 2008 Miha Leban je leta 2007 diplomiral na Fakulteti za elektrotehniko, Univerze v Ljubljani. Od takrat je zaposlen na Soąkih elektrarnah Nova Gorica. Od leta 2008 je tudi ątudent magistrskega ątudija elektroenergetike na isti fakulteti. Rajko Volk je diplomiral leta 1986 na univerzitetnem študiju Fakultete za elektrotehniko v Ljubljani. Trenutno je vodja sektorja proizvodnje na Soških elektrarnah Nova Gorica in sodeluje pri različnih študijah s področja optimizacije proizvodnje električne energije. Rafael Mihalič je diplomiral leta 1986, magistriral leta 1989 in doktoriral leta 1993 na Fakulteti za elektrotehniko in računalništvo v Ljubljani. Po diplomi je postal asistent na omenjeni fakulteti na Katedri za elektroenergetske sisteme in naprave. Med letoma 1988 in 1991 je bil zaposlen pri SIEMENS AG Erlangen na inštitutu za razdeljevanje električne energije in načrtovanje omrežij. Od leta 2005 je redni profesor na Fakulteti za elektrotehniko v Ljubljani. Je član CIGRE, član IEEE, predsednik ŠK B4 SLOKO CIGRE. Področje delovanja vključuje predvsem analizo elektroenergetskih sistemov in naprav FACTS.