Hydromodelovanie vážskej kaskády: Rozdiel medzi revíziami

Z Kiwiki
Skočit na navigaci Skočit na vyhledávání
d
 
(11 medziľahlých úprav od 2 ďalších používateľov nie je zobrazených)
Riadok 1: Riadok 1:
 
[[Kategória:Študentské práce]][[Kategória:Diplomové práce]][[Kategória:Modelovanie]]
 
[[Kategória:Študentské práce]][[Kategória:Diplomové práce]][[Kategória:Modelovanie]]
{{Praca_uvod|2|Návrh a realizácia hydromodelovania vodných elektrární vážskej kaskády|Základné teoretické rovnice prúdenia v otvorených korytách|Hydromodelovanie vážskej kaskády}}
+
{{Praca_uvod|3|Návrh a realizácia hydromodelovania vodných elektrární vážskej kaskády|Základné teoretické rovnice prúdenia v otvorených korytách|Opis modelovaného systému a jednotlivých častí hydromodelu|Hydromodelovanie vážskej kaskády}}
 
__TOC__
 
__TOC__
 +
 
= =
 
= =
==Úloha hydromodelu a opis jeho jednotlivých častí==
+
Modelovanie a simulácia je jednou z hlavných metód ako overovať chovanie analyzovaných systémov. Matematický model je abstrakcia pretransformovaná do počítačového programu, ktorý je podobný pôvodnému objektu, na ktorom vykonávame experimenty. Modelový experiment je riešenie vhodne formulovanej matematickej úlohy s rešpektovaním vzťahov opisujúcich model pre konkrétne podmienky. Komunikácia užívateľa s programom prebieha cez príkazy a cez zadanie počiatočných a koncových podmienok konkrétnych parametrov. Výsledky sú z dôvodu prehľadnosti zobrazené v grafickej aj číselnej forme. Optimalizovaným modelovaním hľadáme také riešenie úlohy, ktoré spĺňa naše predpoklady výsledku, prípadne nám umožní sa vyhnúť problémovým situáciám v modelovanom systéme. Našou úlohou, ktorú budeme ďalej rozoberať, je optimalizované modelovanie vážskej kaskády. V druhej kapitole boli uvedené všetky časti a moduly súvisiace so samotným hydromodelovaním. Postupnosť krokov modelovania je nasledovná:
Tento programový nástroj je matematicko-logickým opisom reálneho objektu, ktorý sa skladá zo skupín, ktorým patria objekty a nakoniec prvky na výrobu elektrickej energie.
+
*Vstup pre zadanie počiatočných a koncových podmienok;
*Skupiny – Stupeň vážskej kaskády
+
*Hydromodelovanie;
*Objekty – Vodná elektráreň
+
*Analýza výsledkov hydromodelovania.
*Prvky – Turbogenerátor
 
  
Hydromodelovanie vážskej kaskády bude v istých prístupoch obsahovať zjednodušenia ako napr.
+
==Vstup pre zadanie počiatočných a koncových podmienok==
*Horná hladina vodnej nádrže bude pre ľubovolný prietok príhaťovou VE pre túto hydrostatická, to znamená že nebude vyhodnocovaná hydraulická strata. Bude to z dôvodu, že nedisponujeme relevantne nameranými hodnotami hydraulickej straty pri rôznych prietokoch a rôznych prevádzkových hladinách vodnej nádrže.
+
Hydromodel je komplexný nástroj pozostávajúci z jednotlivých modulov, ktoré naprogramované v prostredí VBA. Práca s nástrojom si však vyžaduje hlbšie znalosti v danej problematike. Pre zadávanie vstupov, počiatočných a koncových podmienok, nám slúži súbor s programom Hydromodel.xls, prostredníctvom ktorého je možné editovať vstupy pre celý model. Vstupy sú kontrolované proti zadaniu nekorektných údajov a zároveň v prípade, že zadaná hodnota presahuje obmedzujúce parametre, model príslušné vstupy vyznačí (podfarbí pozadie a zmení farbu písma) a zablokuje spustenie hydromodelovania. Takéto zadávanie vstupov môžeme nazvať - zadávanie vstupov naslepo.
*Takisto nebudú vyhodnocované straty v tlakových privádzačoch na VE, nakoľko nie dostupné konštrukčné parametre týchto objektov. Chyba spôsobená zanedbaním tejto hydraulickej straty bude malá a kompenzovaná hydraulickými stratami vyhodnotenými v kanáloch.
 
*Nebudú uvažované zanesenia hrubých hrablíc na haťových objektoch pri vstupe do prívodných kanálov pre VE a takisto jemných hrablíc priamo na prívodnom objekte VE (vstup do špirály).
 
 
Súčasťou hydromodelu bude hydrologický, optimalizačný, hydraulický a transformačný model.
 
*Hydrologický model predstavuje systém medzipovodí, do ktorých vstupujú parametre prítoku z vyššie postavených nádrží, bočné prítoky a odtoky z medzipovodí s energetickým alebo neenergetickým využitím. Maximálne a minimálne prevádzkové hladiny s možnosťou obmedzenia hladín v špecifických prevádzkových prípadoch. Takisto aj obmedzenia v dostupnosti energetických zariadení (GO, BO, REK, EZ, IV). Súčasťou hydrologického modelu budú mapy limitných prietokov pre kaskády a elektrárne prislúchajúcich vodných nádrží, ktoré zodpovedajú hladinám na vodných nádržiach a kanáloch, raster s krokom 0,05m pre celý rozsah prevádzkových hladín. Tieto prietokové mapy budú vyhodnotené od povodňových vôd až po bežné prietokové režimy (hltnosť vodnej elektrárne, derivačného kanála), ktoré dokáže daný stupeň spracovať. Súčasťou budú takisto dotokové časy vody pre kanál a jalové vody pre staré korytá medzi jednotlivými stupňami. Výstupom tohto modelu budú údaje o prietokoch cez jednotlivé stupne vážskej kaskády a vytvarovanie priebehu hladiny na jednotlivých nádržiach počas 24 hodín v hodinovom rastri (max. 25 hodín).
 
*Optimalizačný model má za úlohu na základe vstupných podmienok a obmedzení z hydrologického modelu navrhnúť prietokové pomery, rešpektujúc okrajové podmienky hladín a prietokov. Ďalej nasadenie jednotlivých turbogenerátorov do prevádzky s obmedzeniami prietoku v závislosti od pozície hladiny v nádrži a rezervu prietoku pre regulačnú službu počas 24 hodín v hodinovom rastri (max. 25 hodín).
 
*Hydraulický model na základe výsledkov predchádzajúcich dvoch modelov navrhne hydraulické nastavenie hladín kanálov jednotlivých stupňov vážskej kaskády. Táto hydraulická závislosť je vyhodnocovaná aj na kanáloch pod poslednou VE danej kaskády, kde odpadný kanál zaúsťuje do nasledujúcej vodnej nádrže prípadne koryta rieky. Horná hladina vodnej nádrže bude hydrostatická bez ovplyvnenia drsností nádrže pri rôznych prietokoch, ale hladina na konci odpadného kanála (poslednej VE) ústiaceho do nádrže, bude vyhodnotená cez minimálny potenciál z DH odpadného kanála pod VE. Výstupom tohto modelu sú teda spádové pomery na jednotlivých elektrárňach.
 
*Transformačný model na základe prietoku, spádu, účinnosti a dostupnosti technologického zariadenia bude dávať výsledok v podobe výkonového zaťaženia pre danú elektráreň a hodinu s presným počtom nasadených turbogenerátorov.
 
  
Úzke previazanie bude prevládať medzi hydrologickým a optimalizačným modelom, ktorých výstupy budú ďalej spracované súčastne hydraulickým a transformačným modelom.
 
  
==Prehľad stupňov vážskej kaskády==
+
Pri spustení (otvorení súboru Hydromodel.xls) sa zobrazí priestor s výsledkami tvorený číselným zobrazením nasadeného výkonu jednotlivých objektov VE, ktorých výkony sumarizuje skupina SVK. Nakoniec súčet všetkých skupín je vyjadrený v položke DDZ Váh a graficky zobrazený spolu s položkou Ocenenie DDZ (Obr. 3.1). Suma ocenenia DDZ vyjadruje prácu a veľkosti hodnôt ocenenia vyjadrujú pomer nasadenia v priebehu celého dňa.
Hydromodel bude pozostávať z nasledujúcich stupňov vážskej kaskády a k nim prislúchajúcimi elektrárňami spolu s počtom turbogenerátorov:
 
*VN Orava – VE Orava (OR) – 2xTG;
 
*VN Tvrdošín – VE Tvrdošín (TV) – 3xTG;
 
*VN Liptovská Mara – VE Liptovská Mara (LM) – 4xTG;
 
*VN Bešeňová – VE Bešeňová (BE) – 2xTG;
 
*VN Krpeľany – VE Krpeľany (KR) – 3xTG, VE Sučany (SU) – 3xTG, VE Lipovec (LI) – 3xTG;
 
*VN Žilina – VE Žilina (ZI) – 2xTG;
 
*VN Hričov – VE Hričov (HR) – 3xTG, VE Mikšová (MI) – 3xTG, VE Považská Bystrica (PB) – 3xTG;
 
*VN Nosice – VE Nosice (NO) – 3xTG;
 
*VN Dolné Kočkovce – VE Ladce (LA) – 2xTG, VE Ilava (IL) – 2xTG, VE Dubnica (DU) – 2xTG, VE Trenčín (TN) – 2xTG;
 
*VN Trenčianske Biskupice – VE Kostolná (KO) – 2xTG,        VE Nové Mesto (NM) – 2xTG, VE Horná Streda (HS) – 2xTG;
 
*VN Sĺňava – VE Madunice (MA) – 3xTG;
 
*VN Kráľová – VE Kráľová (KA) – 2xTG.
 
  
[[Súbor:rb_dp2.1.jpg|framed|center|Obr. 2.1 Schéma vodných nádrží a stupňov vážskej kaskády]]
+
[[Súbor:Narh_3.1.png|framed|center|Obr. 3.1 Ovládací panel]]
  
==Rozdelenie VD na základné skupiny==
+
V súčasnej dobe nástroj nedisponuje databázou vzhľadom k rozsahu aktuálne riešenej úlohy. Rozšírenie bude vyžadovať analýzu a následne tvorbu vhodne nastaveného manažmentu dát, ktorý zabezpečí rozumné riešenie pre použitie rozšírenia na týždennú, mesačnú a ročnú prípravu prevádzky a väzbu na riadiaci systém. V súčasnej dobe databázu nahrádza interná pamäť, ktorá zabezpečuje potrebný manažment údajov v rámci D-1, D a D+1. Plánovanie nasadenia výkonu počas dňa pomocou modelu vieme realizovať v rozsahu 23 až 25 hodín.
Spôsob prevádzky na jednotlivých vodných dielach z hľadiska hospodárenia s vodou rozdeľujeme na vodné diela
 
;Akumulačné:sú prevádzkované cez hladinou vo vyrovnávacej nádrži, z ktorej je odtok vody riadený prietokovou reguláciou
 
;Regulačné:sú prevádzkované v medziach prevádzkových hladín danej nádrže a odtoky sú riadené výkonovou reguláciou (v tomto prípade je dôležitý výkon aj napriek tomu, že výkon sa riadi cez prietok)
 
  
Medzi akumulačné VD patrí VN Orava a Liptovská Mara. Perióda prázdnenia a plnenia akumulačnej vodnej nádrže je jeden rok a vyrovnávacej vodnej nádrže, ktorá sa nachádza pod nimi je táto perióda jeden deň.
 
  
[[Súbor:rb_dp2.2.png|framed|center|Obr. 2.2  Schéma vodného diela akumulačného typu]]
+
Do jednotlivých medzipovodí sú prítoky vstupom, ktorý kryje časť výroby. Model sám nedokáže na požadovanú výrobu samočinne nastaviť odtoky z vyrovnávacích nádrží TV a BE tak, aby kryl ostatnú časť výroby. Na tento účel je modul, súčasť hydromodelu, ktorý umožňuje aproximovať výrobu s presnosťou, ktorú rozoberieme neskôr. Modul sprehľadňuje zadávanie vybraných vstupov a v okamihu zadania dáva odozvu v podobe odhadu výroby. Modul sa volá zo súboru Hydromodel.xls. Pri tomto volaní si odovzdajú informácie o aktuálnych pomeroch prietoku a výkonu. V prípade veľkých vôd, keď dochádza k odstaveniu skupín SVK a prietokové, výkonové pomery sa rovnajú nule, modul s názvom Aproximácia_výroby_Váh.xls, ktorý aproximuje výrobu, použije prednastavené technické parametre TG jednotlivých objektov VE. Prostredie modulu je na obrázku 3.2. Pri odovzdávaní informácii o pomeroch modul aproximácie výroby preberie z riadiacej časti aj prítoky za D-1 z jednotlivých medzipovodí a v rámci výpočtu prostredníctvom aktuálnych prítokov z D-1 a priemerných hodnôt prietoku cez skupiny SVK počas dňa D stanový najlepší odhad výroby pre hydromodel. Tak ako riadiaca časť programu aj modul aproximácie používa časť hydraulickej časti pre vodné nádrže a modul pre dotok vody. Modul aproximácie integruje prvky prehľadného nastavenia celého riadiaceho súboru Hydromodel.xls.
  
Medzi regulačné VD patrí VN Tvrdošín, Bešeňová, Krpeľany, Žilina, Hričov, Dolné Kočkovce, Trenčianske Biskupice, Madunice, tu je perióda jeden deň a VN Nosice, Kráľová tu je perióda prázdnenia a plnenie jeden týždeň.
+
[[Súbor:Narh_3.2.png|framed|center|Obr. 3.2 Panel aproximácie výroby]]
  
Čas periód plnenia a prázdnenia je závislý od veľkosti vodnej nádrže a od významu pre sústavu celej vážskej kaskády. Nádrž z ročným vyrovnaním hladín sa nachádza nad všetkými VN vážskej kaskády a jej význam je v čase prebytku vody (z dažďov, topenia snehov) túto akumulovať a v čase nedostatku vody (letné mesiace, zimné mesiace) v jednotlivých medzipovodiach dotovať zvýšenými odtokmi. Nádrž z týždenným vyrovnaním (ako je VD NO) zabezpečuje vodu pre nižšie postavené stupne vážskej kaskády na pološpičkovú (priebežnú až špičkovú) alebo špičkovú (nasadenie je realizované v čase dňa) prevádzku a zároveň súži na zachytenie vody zo špičkovej prevádzky vyššie postaveného stupňa. Zároveň zachytáva vodu z prívalových letných búrok a tým zmierňuje ekonomické dopady na nevyužitý hydroenergetický potenciál v tejto časti Váhu. Prevádzka nádrží z denným vyrovnaním hladín je závislá od momentálnej situácie bočných prítokov do jednotlivých medzipovodí, kde sa tieto nádrže nachádzajú. V čase kedy je možné túto vodu zachytiť sa nasadzovanie orientuje do denných hodín. Inak sa na danom stupni realizuje pološpičková v špecifických prípadoch priebežná prevádzka (povodne, prietokové lebo hladinové režimy).  
+
Jednotlivé položky sú:
 +
#Minimálny počet nasadenia TG skupiny SVK do prevádzky;
 +
#Maximálny počet nasadenia TG skupiny SVK do prevádzky;
 +
#Maximálny prietok cez skupiny SVK;
 +
#Hodnotu počiatočnej hladiny VN;
 +
#Hodnotu plánovanej koncovej hladiny VN;
 +
#Prítoky vody do jednotlivých medzipovodí;
 +
#Biologické odtoky z medzipovodí do koryta;
 +
#Priemerný denný prietok cez SVK;
 +
#Stanovenie špičkovej alebo priebežnej prevádzky skupín SVK;
 +
#CheckBox tlačidlá stanovujeme prepúšťanie Qj do hltnosti kanálov cez kanál alebo hať;
 +
#Prepínač umožňujúci prenášať obmedzenia oblastí 1, 2, 3;.
 +
#CheckBox povolenia priebežného výpočtu po zadaní hodnôt v oblastí 4, 5, 6, 7.
  
[[Súbor:rb_dp2.3.png|framed|center|Obr. 2.3  Schémy vodných diel (príhaťová, kanálová) regulačného typu]]
 
  
Z pohľadu prevádzky VE môžeme tieto rozdeliť na:
+
Celý nástroj zabezpečuje integritu dát, nekorektné vstupy výstražne označí a v prípade odoslania (napr. záporný priemerný prietok SVK) blokuje odoslanie. Samotná presnosť aproximácie je závislá od oblasti prietokov, v ktorej sa modeluje. Pri prevádzke do hltnosti kanálov a mierne nad túto hltnosť je odhad uspokojivý. V oblasti prietokov veľkých vôd, kde nastáva aj odstavenie SVK z prevádzky a pomerové hodnoty prietoku a výkonu sú načítané z technických parametrov TG objektov VE, je odhad slabý. V tomto prípade sa využíva prednosť v jednoduchosti zadávania údajov, nie presnosť odhadu.
;Príhaťové elektrárne:nachádzajú sa v tesnej blízkosti vodnej nádrže a v rámci hydromodelovania horná hladina na týchto VE sa bude považovať za hydrostatickú. Pod týmito VE sa môže nachádzať buď odpadný kanál alebo staré koryto.
 
;Kanálové elektrárne:nachádzajú sa na prívodnom kanály z vodnej nádrže alebo vyššie postavenej VE a vodu odvádzajú odpadným kanálom k nasledujúcej VE, prípadne do nádrže, alebo starého koryta.
 
  
==Hydrologický model==
 
Tento model poskytuje údaje o prietokoch cez jednotlivé stupne vážskej kaskády. Poskytuje údaje o disponibilite plnenia a prázdnenia vodných nádrží, a tiež o všetkých prietokoch do nich vstupujúcich a z nich vystupujúcich. Limitné hodnoty prietokov, ktoré je možné previesť derivačnými kanálmi a spracovať tak cez VE boli vyhodnotené do prietokových máp prostredníctvom vytvoreného nástroja hydraulicko – transformačného modelu. Úlohou týchto máp je popísať limitné prietoky pre optimalizačnú metódu ešte pred samotným odoslaním do hydraulicko – transformačného modelu, aby sa nestal stav, ktorý by mal za následok vrátenie príliš veľkého množstva vody (prebytok alebo deficit) späť na prehodnotenie a úpravu obmedzujúcich prietokov do optimalizačnej metódy. Ak by sme s limitnými prietokmi neuvažovali, pri samotnom spustení optimalizovania a riešenia celej úlohy, mohli by sme navodiť niekoľko stavov, ktoré by spôsobili enormné predĺženie času spracovania alebo velmi skreslené výsledky. Tým by sa stal celý nástroj nepoužiteľný!
 
  
Väzby v hydrologickom modely vážskej kaskády môžeme rozdeliť do bodov:
+
Príklad presností odhadov:
#Prietoky vstupujúce a vystupujúce v danom profile (medzipovodí);
+
Zaujímavosť zistená pri analýze odpúšťania vody z VN TV na krytie výkonu Váhu je v tom, že plné výkonové nasýtenie Váh dosiahne približne za štyri dni. Hladiny sú pritom modelované tak, aby počiatočná sa rovnala koncovej hladine. Dotoky z dňa D model po ukončení modelovania SVK pripraví pre D+1. V modelovom príklade nebudeme uvažovať jalovú vodu starým korytom. Všetky medzipovodia majú nastavený bočný prítok rovnajúci sa biologickému odtoku a odtok z vyrovnávacích nádrží TV, BE je nastavený na minimálny odtok.
#Dotokové časy vody na úseku medzi jednotlivými stupňami v bežnej prevádzke VE;
 
#Dotokové časy jalovej vody cez hať a staré koryto do nasledujúceho profilu (medzipovodia).
 
  
===Prietoky vstupujúce a vystupujúce v danom profile (medzipovodí)===
+
'''Tab. 3.1 Výkonová nasýtenosť VE Váh'''
Riadenie prietoku cez VE na akumulačnom VD je závislé od odtoku z vyrovnávacej vodnej nádrže a tiež od objemu plnenia alebo prázdnenia tejto nádrže.  
 
[[Súbor:rb_dp2.4.png|framed|center|Obr. 2.4 Grafické znázornenia jednotlivých prietokov pre akumulačné VD]]
 
Základnú rovnicu pre riadenie vody typu akumulačného VD môžeme vyjadriť ako
 
  
{{vzorec|
+
{|border="2" cellspacing="0" cellpadding="4" width="100%"
<math>\overline{{{Q}_{VE}}}=\frac{\frac{\left( {{V}_{kon\operatorname{cov}\acute{y}}}-{{V}_{po\check{c}iato\check{c}n\acute{y}}} \right)}{3600}+\sum\limits_{krok=1}^{t}{{{Q}_{OdtokVN\left( krok \right)}}}-\sum\limits_{krok=1}^{t}{{{Q}_{MedziPovodie\left( krok \right)}}}-\sum\limits_{krok=1}^{t}{{{Q}_{Ha\check{t}\left( krok \right)}}}}{{{t}_{hodiny}}}</math>|2.1}}
+
|align = "center"|'''P.č.'''
 +
|align = "center"|'''Aproximácia'''<br><nowiki>[</nowiki>MWh<nowiki>]</nowiki>
 +
|align = "center"|'''Hydromodel'''<br><nowiki>[</nowiki>MWh<nowiki>]</nowiki>
 +
|align = "center"|'''Rozdiel'''<br><nowiki>[</nowiki>MWh<nowiki>]</nowiki><br>
 +
|align = "center"|'''Čas spracovania úlohy'''<br><nowiki>[</nowiki>sek.<nowiki>]</nowiki>
  
kde:
 
{|
 
 
|-
 
|-
|t
+
|1.
|časový interval výpočtu [hod];
+
|align = "center"|1517,6
 +
|align = "center"|1693,7
 +
|align = "center"|176,1
 +
|align = "center"|16,55
 +
 
 
|-
 
|-
|krok
+
|2.
|krok výpočtu [hod];
+
|align = "center"|3280,4
 +
|align = "center"|3383,7
 +
|align = "center"|103,3
 +
|align = "center"|25,14
 +
 
 
|-
 
|-
|<math>\overline{{{Q}_{VE}}}</math>
+
|3.
|priemerný prietok vody cez VE z akumulačnej nádrže za interval t [m3.s-1];
+
|align = "center"|4699,2
 +
|align = "center"|4827,9
 +
|align = "center"|128,7
 +
|align = "center"|17,94
 +
 
 
|-
 
|-
|<math>V_{Počiatočný}</math>
+
|4.
|objem vyrovnávacej nádrže na začiatku intervalu t [<math>m^3</math>];
+
|align = "center"|4928,4
 +
|align = "center"|5066,7
 +
|align = "center"|138,3
 +
|align = "center"|18,88
 +
 
 
|-
 
|-
|<math>V_{Koncový}</math>
+
|5.
|objem vyrovnávacej nádrže na konci intervalu t [<math>m^3</math>];
+
|align = "center"|4946,4
|-
+
|align = "center"|5110,8
|<math>ΣQ_{OdtokVN}</math>
+
|align = "center"|164,4
|suma plánovaných odtokov z vyrovnávacej nádrže pre nadlepšovanie výroby za interval t [<math>m^3.s_{-1}</math>];
+
|align = "center"|20,59
|-
 
|<math>ΣQ_{MedziPovodie}</math>
 
|suma predikcie bočného prítoku vody do profilu vyrovnávacej vodnej nádrže za interval t [<math>m^3.s_{-1}</math>];
 
|-
 
|<math>ΣQ_{Hať}</math>
 
|suma jalových odtokov z akumulačnej nádrže cez hať za interval t [<math>m^3.s_{-1}</math>];
 
|}
 
 
 
[[Súbor:rb_dp2.5.png|framed|center|Obr. 2.5  Grafické znázornenia jednotlivých prietokov pre regulačné VD]]
 
 
 
===Dotokové časy vody na úseku medzi jednotlivými stupňami v bežnej prevádzke VE===
 
Tento faktor v podstatnej miere ovplyvňuje časovú závislosť nasadenia VE medzi jednotlivými stupňami vážskej kaskády, ktoré majú priamu, ale aj nepriamou hydraulickou väzbou. Pri nasadení VE to znamená významné ovplyvnenie chovania sa hladín vodnej nádrže a tým aj prevádzky VE v čase, kým voda je ešte len na ceste z vyššieho stupňa do daného profilu. Úseky derivačných kanálov nie je potrebné zaťažovať týmto časom (napriek tomu, že dotokové časy sú aj v kanáloch) z dôvodu rovnakého prietoku v čase cez prvú a poslednou VE daného stupňa vážskej kaskády.
 
 
 
Medzi úseky, ktoré sú týmto významne ovplyvnené patria:
 
*VE Tvrdošín - VN Krpeľany;
 
*VE Bešeňová - VN Krpeľany;
 
*VE Lipovec - VN Žilina;
 
*VE Žilina - VN Hričov;
 
*VE Považská Bystrica - VN Nosice;
 
*VE Nosice - VN Dolné Kočkovce;
 
*VE Trenčín - VN Trenčianske Biskupice;
 
*VE Horná Streda - VN Sĺňava;
 
*VE Madunice - VN Kráľová.
 
 
 
[[Súbor:rb_dp2.6.png|framed|center|Obr. 2.6  Grafické znázornenie úseku medzi stupňami]]
 
 
 
Riešenie problému dotokových časov a transformáciu prietokovej vlny budeme realizovať pomocou vrstvových posunov v diskrétnom tvare. Tieto vrstvy budú vypočítané funkciou lineárnej interpolácie z bodov maximálneho a minimálneho prietoku a ich časovej závislosti dotoku z vyššie postavenej VE do nasledujúcej vodnej nádrže.
 
 
 
Pre priblíženie problému uvediem príklad, v ktorom bude stanovený odtok v každej hodine dňa z VE Madunice. Následne je každý odtok vyhodnotený funkciou
 
<math>{{Q}_{t\left( krok \right)}}={{t}_{Dotok}}\left( {{Q}_{\Pr itok\left( krok \right)}} \right)</math> , ktorá rozdelí odtok do časových vrstiev, v ktorých sa budú pohybovať k stanovenému cieľu. Výsledkom transformácie odtoku nám vznikne prítok s posunom jednotlivých časových vrstvách a teda aj prietokov pre jednotlivé hodiny dňa. Celý príklad ilustruje obrázok 2.7 a 2.8.
 
 
 
[[Súbor:rb_dp2.7.png|framed|center|Obr. 2.7 Príklad rozloženia prietoku cez VE do časových vrstiev a jeho transformácia do nasledujúcej VN]]
 
  
Veľmi dôležitý parameter bude pre nás suma odtoku a suma prítoku, ktoré sa musia rovnať za rovnaký časový interval v opačnom prípade vnesieme chybu do celého výpočtu. Pomocou spomínaného riešenia sa dokážeme priblížiť k reálnemu priebehu prietokovej vlny, ktorej priebeh a časový posun je transformovaný do tvaru ako je to uvedené na obr. 2.8.
 
 
[[Súbor:rb_dp2.8.png|framed|center|Obr. 2.8 Grafické znázornenie – odtok z VE, transformácia v čase a čas prítoku do nasledujúcej VN]]
 
 
Vážska kaskáda má rôzne parametre dotokov vody medzi jednotlivými objektmi a stupňami, preto budú uvedené v tabuľke 2.1.
 
 
{|class=wikitable
 
|+Tab. 2.1  Hodnoty odtokov a&nbsp;k&nbsp;nim prislúchajúci čas dotoku vody medzi objektmi
 
|-
 
|'''Úsek medzi objektom a stupňom'''
 
|'''Prietok <nowiki>[</nowiki>m'''<sup>'''3'''</sup>.s<sup>-1</sup>'''<nowiki>]</nowiki>'''
 
|'''Čas dotoku t <nowiki>[</nowiki>hod<nowiki>]</nowiki>'''
 
|&nbsp;
 
|'''Prietok <nowiki>[</nowiki>m'''<sup>'''3'''</sup>.s<sup>-1</sup>'''<nowiki>]</nowiki>'''
 
|'''Čas dotoku t <nowiki>[</nowiki>hod<nowiki>]</nowiki>'''
 
|-
 
|VE Orava - VE Tvrdošín
 
|800
 
|0
 
|&nbsp;
 
|15
 
|0
 
|-
 
|VE Liptovská Mara - VE Bešeňová
 
|800
 
|0
 
|&nbsp;
 
|50
 
|0
 
|-
 
|Lipovec - Žilina
 
|210
 
|1,4
 
|&nbsp;
 
|30
 
|3,5
 
|-
 
|Žilina - Hričov
 
|720
 
|0,3
 
|&nbsp;
 
|60
 
|0,75
 
|-
 
|Považská Bystrica - Nosice
 
|500
 
|0,3
 
|&nbsp;
 
|50
 
|0,75
 
|-
 
|Nosice - Kočkovce
 
|390
 
|0,3
 
|&nbsp;
 
|40
 
|0,75
 
|-
 
|Trenčín - Trenčianske Biskupice
 
|160
 
|0,3
 
|&nbsp;
 
|40
 
|0,75
 
|-
 
|Horná Streda - Dráhovce
 
|180
 
|0,6
 
|&nbsp;
 
|40
 
|1,5
 
|-
 
|Madunice - Kráľová
 
|300
 
|2,8
 
|&nbsp;
 
|40
 
|7
 
 
|}
 
|}
  
===Dotokové časy jalovej vody cez hať a staré koryto do nasledujúceho profilu (medzipovodia)===
 
V čase veľkých vôd, kedy sú už kapacity energetického využitia vyčerpané je potrebné ostatnú vodu previesť neenergeticky cez hať do koryta rieky. Má podobne časový charakter ako pri čisto energetickej prevádzke VE. Riešený bude však úsek medzi jednotlivými stupňami vážskej kaskády. Prevádzka VE je pri dosiahnutí určitých prietokov cez hať do koryta rieky po časovom intervale dotoku vody ďalej obmedzená a je to spôsobené vzdutím dolnej hladiny pod VE, kedy kontrolované parametre spádových pomerov nie je možné dodržať. V takomto prípade bude hydraulický model znižovať prietokovú kapacitu turbogenerátorov na udržanie predpísaných spádových pomerov na jednotlivých VE.
 
  
Medzi úseky, ktoré sú týmto významne ovplyvnené patria:
+
V tabuľke 3.1 prvý bod je počiatočný stav Váhu s minimálnymi prítokmi do medzipovodí a odtokmi z vyrovnávacích nádrží. V druhom bode je zmenený iba odtok z VN TV z 4 m3.s-1 skokovo na 65 m3.s-1. Je vidieť postupné výkonové nasycovanie na Váhu a plné nasýtenie nastáva na štvrtý deň. Pri aproximácii výroby je v tomto príklade rozdiel v rozsahu od 100 MWh do 180 MWh, čo môžeme považovať za uspokojivý výsledok. Chyba je spôsobená spriemerovaním odtokov z SVK počas dňa, pričom hydromodel nasadzuje prietok do špičiek a väčšina objemu ostáva v dni D a nie je presunutá do dňa D+1. Pokiaľ by modul aproximácie nevyužíval model dotokov vody, hodnota rozdielu odhadu voči hodnote z hydromodelu by mohla presiahnuť 1500 MWh v čase väčších zmien prítokov, čo je neprípustné!
*VD Tvrdošín - VN Krpeľany;
 
*VD Bešeňová - VN Krpeľany;
 
*VD Krpeľany - VN Žilina;
 
*VD Žilina - VN Hričov;
 
*VD Hričov - VN Nosice;
 
*VD Nosice - VN Dolné Kočkovce;
 
*VD Dolné Kočkovce - VN Trenčianske Biskupice;
 
*VD Trenčianske Biskupice - VN Sĺňava;
 
*VD Dráhovce - VN Kráľová.
 
  
[[Súbor:rb_dp2.9.png|framed|center|Obr. 2.9  Grafické znázornenie úseku medzi vodnými nádržami stupňov vážskej kaskády]]
 
  
Riešenie problému dotokových časov a transformácia prietokovej vlny bude realizované rovnakou metódou ako v kap. 2.3.2 pomocou vrstvových posunov v diskrétnom tvare.
+
Odľahčenie výroby namodelujeme analogicky. VN TV nastavíme pôvodný odtok 4 m3.s-1. Tabuľka 3.2 opisuje tak, ako v prípade nasýtenia výkonu v prvom bode, počiatočný stav Váhu a v nasledujúcom bode zníženie odtoku z VN TV. Rozdiel odhadu je pri odľahčovaní v rozsahu od 160 MWh do 220 MWh, čo môžeme opäť považovať za uspokojivé. Opäť dostávame štyri dni do úplného odľahčenia výroby VE.
Nasledujúci príklad je počítaný medzi stupňami s jalovým prepadom cez haťový objekt VD DR (naväzuje na predchádzajúci príklad). Výsledkom transformácie odtoku nám vznikne prítok s posunom jednotlivých časových vrstiev a teda aj prietokov pre jednotlivé hodiny dňa. Celý príklad ilustruje obrázok 2.10 a 2.11.
 
  
[[Súbor:rb_dp2.10.png|framed|center|Obr. 2.10 Príklad rozloženia prietoku cez hať do časových vrstiev a jeho transformácia do nasledujúcej VN]]
+
'''Tab. 3.2  Odľahčenie výroby VE Váh'''
  
Rovnako ako v kapitole 2.4.2 suma odtoku a suma prítoku sa musí rovnať za rovnaký časový interval v opačnom prípade vnesieme chybu do celého výpočtu.
+
{|border="2" cellspacing="0" cellpadding="4" width="100%"
 +
|align = "center"|'''P.č.'''
 +
|align = "center"|'''Aproximácia'''<br><nowiki>[</nowiki>MWh<nowiki>]</nowiki>
 +
|align = "center"|'''Hydromodel'''<br><nowiki>[</nowiki>MWh<nowiki>]</nowiki>
 +
|align = "center"|'''Rozdiel'''<br><nowiki>[</nowiki>MWh<nowiki>]</nowiki><br>
 +
|align = "center"|'''Čas spracovania úlohy'''<br><nowiki>[</nowiki>sek.<nowiki>]</nowiki>
  
[[Súbor:rb_dp2.11.png|framed|center|Obr. 2.11  Grafické znázornenie – odtok z hate, transformácia v čase a čas prítoku do nasledujúcej VN]]
 
 
Vážska kaskáda má rôzne parametre dotokov vody medzi jednotlivými stupňami, preto budú uvedené v tabuľke 2.2.
 
 
{|class=wikitable
 
|+Tab. 2.2  Hodnoty odtokov a k nim prislúchajúci čas dotoku vody medzi stupňami
 
 
|-
 
|-
!Úsek medzi objektom a stupňom
+
|1.
!Prietok <nowiki>[</nowiki>m<sup>3</sup>.s<sup>-1</sup><nowiki>]</nowiki>
+
|align = "center"|4946,4
!Čas dotoku t <nowiki>[</nowiki>hod<nowiki>]</nowiki>
+
|align = "center"|5110,8
!&nbsp;
+
|align = "center"|164,4
!Prietok <nowiki>[</nowiki>m<sup>3</sup>.s<sup>-1</sup><nowiki>]</nowiki>
+
|align = "center"|20,59
!Čas dotoku t <nowiki>[</nowiki>hod<nowiki>]</nowiki>
 
|-
 
|VD Tvrdošín - VN Krpeľany
 
|800
 
|4,5
 
|&nbsp;
 
|5
 
|11
 
  
 
|-
 
|-
|VD Bešeňová - VN Krpeľany
+
|2.
|700
+
|align = "center"|3195,2
|2,5
+
|align = "center"|3412,9
|&nbsp;
+
|align = "center"|217,7
|15
+
|align = "center"|29,38
|6,1
 
  
 
|-
 
|-
|VD Krpeľany - VN Žilina
+
|3.
|700
+
|align = "center"|1828,5
|3,4
+
|align = "center"|2031,7
|&nbsp;
+
|align = "center"|203,2
|20
+
|align = "center"|20,41
|8,3
 
  
 
|-
 
|-
|VD Žilina - VN Hričov
+
|4.
|700
+
|align = "center"|1526,1
|0,3
+
|align = "center"|1709,3
|&nbsp;
+
|align = "center"|183,2
|20
+
|align = "center"|17,41
|0,8
 
  
 
|-
 
|-
|VD Hričov - VN Nosice
+
|5.
|700
+
|align = "center"|1514,6
|3,5
+
|align = "center"|1693,7
|&nbsp;
+
|align = "center"|179,1
|20
+
|align = "center"|16,98
|8,6
 
  
|-
 
|VD Nosice - VN Dolné Kočkovce
 
|700
 
|0,6
 
|&nbsp;
 
|20
 
|1,5
 
 
|-
 
|VD D. Kočkovce - VN T. Biskupice
 
|700
 
|3,5
 
|&nbsp;
 
|20
 
|8,6
 
 
|-
 
|VD T. Biskupice - VN Sĺňava
 
|700
 
|4,5
 
|&nbsp;
 
|20
 
|11
 
 
|-
 
|VD Dráhovce - VN Kráľová
 
|700
 
|4
 
|&nbsp;
 
|20
 
|9,8
 
|}
 
 
Ak by sme spojili oba príklady do jedného vznikol by nám priebeh zobrazený na obrázku 2.12. Na celej vážskej kaskáde sú derivačné kanále a staré korytá Váhu zaústené velmi blízko cieľa dotoku vody do vodnej nádrže okrem zaústenia odpadného kanála pod VE Lipovec a VE Madunice. Pre zjednodušenie tu bude prijatý rovnaký prístup ako pri ostatných stupňoch.
 
 
[[Súbor:rb_dp2.12.png|framed|center|Obr. 2.12  Grafické znázornenie – odtok z VE a hate, transformácia v čase a čas prítoku do nasledujúcej VN]]
 
 
==Hydraulický model==
 
 
Údaje vystupujúce z hydrologického modelu sú vstupnými pre hydraulický model. Úlohou hydraulického modelu je výpočet spádových pomerov a prietokových obmedzení, ktoré sú v súčinnosti s transformačným modelom ďalej vyhodnotené a doplnené do limitných prietokových máp. Pomocou limitných prietokových máp bude zabezpečená na úrovni optimalizovania korektnosť výpočtu maximálneho prietoku cez VE alebo derivačný kanál. Bude tým zabezpečená rýchla priechodnosť a konečnosť výpočtu. Hydraulický model bude mať za úlohu:
 
*Nastaviť pozíciu hladiny v nádrži;
 
*Určiť hydraulické väzby medzi jednotlivými VE;
 
*Vyhodnotiť hladinové režimy v derivačných kanáloch;
 
*Stanoviť maximá energeticky využiteľného prietoku pomocou predchádzajúcich bodov (mapa limitných prietokov).
 
 
Hydraulický model môžeme rozdeliť na nasledujúce úseky:
 
#Nádrž;
 
#Úsek nádrž – VE;
 
#Úsek VE – VE;
 
#Úsek VE – nádrž;
 
#Úsek VE – rieka.
 
 
===Hydraulický model nádrže===
 
Model pre nádrž má za úlohu stanoviť na základe počiatočnej a koncovej hladiny v nádrži volný objem (pre plnenie) a zásobný objem (pre prázdnenie), ktorý treba započítať k prietokom z hydrologického modelu, čím je ďalej stanovený priemerný denný odtok z daného stupňa vážskej kaskády. Priebežne bude z vypočítaných objemov v jednotlivých hodinách určovať pozíciu hladiny v nádrži. Teoreticky budeme uvažovať s hydrostatickou hladinou v nádrži pri rôznych prietokoch z dôvodu zjednodušenia. Deformáciu nádrže budeme riešiť nepriamo, ale táto metóda bude bližšie rozpracovaná pri úseku v kap. 2.5.5 a 2.5.6. Základné rovnice hydraulického modelu nádrže budú:
 
 
{{vzorec|<math>{{V}_{Z\acute{a}sobn\acute{y}\ \left( krok \right)}}={{V}_{Z\acute{a}sobn\acute{y}\ \left( krok-1 \right)}}+{{Q}_{Celkov\acute{y}\ pr\acute{i}tok\left( krok \right)}}\cdot 3600-{{Q}_{Celkov\acute{y}\ odtok\left( krok \right)}}\cdot 3600</math>|2.3}}
 
{{vzorec|<math>H{{H}_{N\acute{a}dr\check{z}e}}=f\left( {{V}_{Z\acute{a}sobn\acute{y}\left( krok \right)}} \right)</math>|2.4}}
 
 
kde:
 
{|
 
|-
 
|krok
 
|krok výpočtu [hod];
 
|-
 
|<math>V_{Zásobný}</math>
 
|objem nádrže na konci kroku výpočtu [m3];
 
|-
 
|<math>V_{Zásobný(krok-1)}</math>
 
|objem nádrže na začiatku kroku výpočtu [<math>m^3</math>];
 
|-
 
|<math>Q_{Celkový prítok}</math>
 
|celkový prítok z vyššieho stupňa na konci kroku výpočtu [<math>m^3.s^{-1}</math>];
 
|-
 
|<math>Q_{Celkový odtok}</math>
 
|celkový odtok z nádrže na konci kroku výpočtu [<math>m^3.s^{-1}</math>];
 
|-
 
|<math>HH_{Nádrž}</math>
 
|pozícia hladiny vodnej nádrže na konci výpočtového kroku [m.n.m Bpv].
 
 
|}
 
|}
  
K dispozícii máme objemové krivky jednotlivých vodných nádrží, z ktorých budeme pomocou lineárnej interpolácie získavať hodnoty objemov a hladín aj na pozíciách, ktoré nie sú v krivkách uvedené. Jednoduchosť inverzie tejto funkcie nám zabezpečí absolútnu rovnosť hodnôt aj po niekoľkonásobnej transformácii objemu na hladinu a naopak pri následnom posune parametrov a cyklickom opakovaní (Obr. 2.13).
 
  
[[Súbor:rb_dp2.13.png|framed|center|Obr. 2.13  Objemová krivka nádrže]]
+
Časy boli namerané na PC s procesorom Intel® Core™2 CPU 6400 @ 2,13 GHz, 1,99 GB RAM, OS Microsoft Windows XP Professional a Office 2003.
  
===Hydraulický model kanálov===
+
==Hydromodelovanie==
Úlohou hydraulického modelu kanálov je stanoviť hydraulické straty spôsobené drsnosťou kanálov a ďalšími faktormi (ako napr. prietok, nábehová hladina). Tieto straty sa premietnu do spádových pomerov pre jednotlivé objekty VE. Môžeme povedať, že derivačné kanále sa nachádzajú pri každej elektrárni s výnimkou VD KA. Preto tvorba tohto modelu bude z hľadiska výpočtu a vyhodnocovania jednotlivých parametrov jednou s najnáročnejších úloh tejto práce.
+
Po úspešnom zápise vstupov do modelu, priamo alebo pomocou aproximácie výroby, máme model pripravený na spustenie. V tejto časti popíšeme jednotlivé interakcie medzi modulmi hydromodelu. Po spustení nasleduje načítanie vstupných dát optimalizačným modelom, ktorý na základe vstupov z kapitoly 2.9 navrhne nasadenie prietoku cez aktuálne počítanú kaskádu. SVK sú postupne spracovávané v smere toku vody tak, aby bola zabezpečená hydrologická a hydraulická väzba.
  
Pri výpočtoch sa bude uvažovať v kanáloch s ustáleným prúdením vody napriek tomu, že pri zmenách prietoku ustálené prúdenie nastáva po dlhšom čase. Takýto prístup bude zvolený z dôvodu hľadania možného nasadenia vodných elektrární s časovým intervalom deň, týždeň, mesiac a rok. Nebudeme teda skúmať dynamiku vody v kanály nad, ani pod VE pri nábehoch, odstaveniach alebo zmenách v prietoku pri prevádzke.
 
  
Základná rovnica hydraulického modelu kanálov môžeme vyjadriť v tvare:
+
Poradie modelovania stupňov vážskej kaskády:
{{vzorec|<math>\Delta {{z}_{(krok)}}=Q_{(krok)}^{2}\left[ \left( 1\pm \zeta  \right)\cdot \frac{\alpha }{2g}\cdot \left( \frac{1}{S_{2(krok)}^{2}}-\frac{1}{S_{1(krok)}^{2}} \right)+\frac{Q_{(krok)}^{2}}{K_{p(krok)}^{2}}\Delta l \right]</math>|2.5}}
+
#SVK OR-TV
 +
#SVK LM-BE
 +
#SVK KR-SU-LI
 +
#SVK ZI
 +
#SVK HR-MI-PB
 +
#SVK NO
 +
#SVK LA-IL-DU-TN
 +
#SVK KO-NM-HS
 +
#SVK MA
 +
#SVK KA
  
kde:
 
{|
 
|-
 
|krok
 
|krok výpočtu v jednej hodine;
 
|-
 
|Δz
 
|diferencia od nábehovej hladiny [m];
 
|-
 
|Q
 
|prietok prívodným alebo odpadným kanálom [<math>m^3.s^{-1}</math>];
 
|-
 
 
|súčiniteľ miestnej straty;
 
|-
 
 
|Coriolisovo číslo;
 
|-
 
|g
 
|gravitačné zrýchlenie [<math>m.s^{2}</math>];
 
|-
 
|S1
 
|prietoková plocha (horný prietokový prierez) [<math>m^2</math>];
 
|-
 
|S2
 
|prietoková plocha (dolný prietokový prierez) [<math>m^2</math>];
 
|-
 
|Kp
 
|Sp, Cp, Rp sú aritmetické priemery plôch [m2], rýchlostných súčiniteľov [<math>m^{0,5}.s^{-1}</math>] a hydraulických polomerov [m] v hornom a dolnom priereze kde
 
<math>{{K}_{p}}={{C}_{p}}\cdot {{S}_{p}}\sqrt{{{R}_{p}}}</math>
 
|-
 
Δl
 
|dĺžka kanála [m].
 
|}
 
  
Súčasťou rovnice sú aj konštrukčné parametre kanála popísané v kap. 1.2.1. Vstupnými parametrami pre výpočet budú prietok cez kanál a nábehová hladina v kanály (ktorú bude možné pre celý rozsah hodín dňa meniť v rozsahu prevádzkových hladín).  
+
Prekročenie prevádzkových hladín optimalizačný model kontroluje prostredníctvom hydraulického modelu vodných nádrží. Prekročenie koriguje a po poslednej hodine dňa hľadá riešenie dosiahnutia plánovanej koncovej hladiny. Počas optimalizácie verifikuje maximálne nasadenie prietoku cez SVK a jednotlivo cez TG, z ktorých vyberá pre výpočet najmenšiu hodnotu maximálneho obmedzenia. Tieto kontroly prebiehajú aj zo strany minimálnych prietokov cez TG. Vo výpočte optimalizácie SVK prebieha priebežný výpočet pohybu nábehových hladín v derivačných kanáloch. V čase veľkých vôd optimalizačný model porovnáva obmedzenia SVK pred prítokom a po prítoku veľkej vody starým korytom. V prípade rozdielu celú optimalizáciu daného stupňa opakuje s tým, že časť vody nad obmedzením prevedie z kanála na hať. Celý proces sa opakuje do chvíle rovnosti obmedzení pred prítokom a po prítoku vody starým korytom pod poslednú VE danej kaskády. Tento stav vidieť na dialógu Priebeh spracovania výpočtu počas hydromodelovania (Obr. 3.3) v položke Hodiny (SVK bez zmeny a cyklicky sa opakujúce  hodnoty tejto položky).
  
===Úsek nádrž - VE===
+
[[Súbor:Narh_3.3.png|framed|center|Obr. 3.3 Dialóg - priebeh spracovania]]
Úsek nádrž - VE môžeme rozdeliť na úseky s prívodným kanálom a bez prívodného kanála na objekt VE. Nábehová hladina pre kanál alebo VE bude vyhodnocovaná z hladiny v nádrži pre daný výpočtový krok.
 
  
Medzi úseky s kanálom patria:
+
Po ukončení optimalizácie sú postupne údaje importované do hydraulicko-transformačného modulu, kde výpočet prebehne veľmi rýchlo. Výpočet nastaví hydrauliku v kanáloch a na prietokové a spádové pomery transformačný model stanový počet strojov a výkonové zaťaženie v danej hodine. Výsledky sú poslané späť do súboru Hydromodel.xls.
*VD Dolné Kočkovce – VE Ladce;
 
*VD Trenčianske Biskupice – VE Kostolná;
 
*VD Dráhovce (VN Sĺňava) – VE Madunice.
 
Medzi úseky bez prívodného kanála patria (nazývajú sa príhaťové VE):
 
*VN Tvrdošín – VE Tvrdošin;
 
*VN Bešeňová – VE Bešeňová;
 
*VN Krpeľany – VE Krpeľany;
 
*VN Žilina – VE Žilina;
 
*VN Hričov – VE Hričov;
 
*VN Nosice – VE Nosice;
 
*VN Kráľová – VE Kráľová.
 
  
[[Súbor:rb_dp2.14.png|framed|center|Obr. 2.14  Pozdĺžny rez prívodného kanála k VE z vodnej nádrže]]
+
==Analýza výsledkov hydromodelovania==
===Úsek VE - VE===
 
Na úseku VE – VE sa nachádza kanál, ktorý je rozdelený na dva (popísané technickou dokumentáciou):
 
*odpadný od VE;
 
*prívodný k VE.
 
  
Súčasťou týchto kanálov budú aj objemové krivky v rozsahu prevádzkových hladín. Významný rozmer nadobudnú vo chvíli, keď diferencia od nábehovej hladiny poklesne pod minimálnu prevádzkovú hladinu pre daný objekt, alebo prietok cez kanál dosiahne také hodnoty, ktoré znížia HH a zároveň zvýšia DH natoľko, že sa nedodrží minimálny spádový parameter pre VE. V takom prípade hydromodel začne upravovať pozíciu nábehovej hladiny. Pri zmene pozície nábehovej hladiny v kanály sa prepočíta plnenie alebo prázdnenie kanála, čo sa premietne do prietokových pomerov cez celý stupeň.
+
Pri prehliadaní výsledkov a kontrole v súbore Hydromodel.xls prechádzame záložkami prostredníctvom tlačidiel na obrázku 3.4.
  
Medzi úseky tohto typu patria:
+
[[Súbor:Narh_3.4.png|framed|center|Obr. 3.4 Panel ovládania - hydromodel]]
*VE Krpeľany – VE Sučany;
 
*VE Sučany - VE Lipovec;
 
*VE Hričov – VE Mikšová;
 
*VE Mikšová – VE Považská Bystrica;
 
*VE Ladce – VE Ilava;
 
*E Ilava – VE Dubnica;
 
*VE Dubnica – VE Trenčín;
 
*E Kostolná – VE Nové Mesto;
 
*VE Nové Mesto - VE Horná Streda.
 
  
[[Súbor:rb_dp2.15.png|framed|center|Obr. 2.15  Pozdĺžny rez prívodného a odpadného kanála medzi VE]]
+
Postupne budeme prechádzať prostredím pre zadávanie údajov. Vstupy, ktoré môžeme editovať na obrázku 3.5, sú podfarbené modrým pozadím. Ostatné údaje nie je možné meniť. Ovládanie celého programu je jednoduché a intuitívne. Pri prepínaní sa zobrazia grafické výstupy. Je možné prepnúť sa na dátovú časť tlačidlom alebo sa prerolovať myšou. Na každej záložke pre SVK je v ľavom rohu príznak spustenia modelu. Pokiaľ je tento príznak v stave Stop model, na záložke je problém nekorektných údajov, ktoré treba opraviť.
  
===Úsek VE - nádrž===
+
[[Súbor:Narh_3.5.png|framed|center|Obr. 3.5 Panel zadávania vstupov a prehliadanie výstupov - hydromodel]]
Úsek VE – nádrž je kľúčový pri určovaní nepriamej deformácie nádrže v bezprostrednej blízkosti odpadného kanála poslednej VE kaskády. Aby sme mohli hydraulicky prepojiť tento úsek s nádržou bolo potrebné vyhodnotiť vzdutie na konci odpadného kanála. Najlepším riešením tohto problému bolo nájsť najmenší potenciál vzdutia DH pod VE. Vstupný parameter nábehovej hladiny v kanály je hladina vodnej nádrže. Hydraulický model hľadá riešenie pre daný prietok, kým dovtedy klesajúca DH nenaberie stúpajúcu tendenciu. V tomto okamihu je vyhodnotený najmenší potenciál DH a nastáva hydraulické spojenie odpadného kanála VE s nádržou nasledujúceho stupňa vážskej kaskády.
 
  
Medzi úseky tohto typu patria:
+
=Záver=
*VE Orava – VN Tvrdošín;
+
V práci sme rozobrali problematiku hydromodelovania vážskej kaskády, ktorej praktický výstup je nástroj na modelovanie. Postupne, ako sme prechádzali jednotlivé kapitoly, mali sme možnosť zoznámiť sa s danou problematikou navrhovania hydrologického, hydraulického, transformačného a optimalizačného modelu. Spolu jednotlivé časti tvoria systém umožňujúci veľmi náročné matematicko-logické výpočty. Rozobrané nové prístupy modelovania v tejto oblasti, oproti stávajúcemu modelu na HED, umožnili hydraulické spojenie jednotlivých skupín SVK, prácu s hladinami v derivačných kanáloch, lepšiu metódu transformácie dotoku vody, modelovanie veľkej vody od bežných prevádzkových prietokov až po prietoky storočných vôd, autonómnosť optimalizácie, integritu údajov a celého modelu. Veľký rozdiel v  nastavovaní hydromodelu oproti pôvodnému na HED je v zadávaní vstupov. Pôvodný model umožňuje nastaviť vstupy odhadom, využívajúc chybu ľudského faktora multinásobná korekcia vstupných dát s postupným približovaním dosiahnutia požadovanej výroby, kontrola nasadenia turbogenerátorov do prevádzky a v mnohých prípadoch nutnosť manuálne tvarovať DDZ. Ak by sme porovnali obidva hydromodely zistíme, že pôvodný je kompilovaný, ale časová náročnosť výpočtu je porovnateľná s časom výpočtu nástroja vytvoreného v aplikácii Excel pomocou interpretovacieho jazyka VBA. Pokiaľ by sme mali vyhodnotiť celkovú časovú náročnosť modelovania jednej úlohy, stávajúci model vyžaduje plné nasadenie užívateľa v rozsahu jednej hodiny prípadne dlhšie. Naproti tomu aktuálne vytvorený model je schopný poskytnúť výsledok úlohy v priebehu pár minút od začiatku modelovania. Ak by sme chceli riešiť úlohu veľkej vody, pôvodný model na HED nie je schopný riešiť takúto úlohu.
*VE Liptovská Mara - VN Bešeňová;
 
*VE Žilina – VN Hričov;
 
*VE Považská Bystrica – VN Nosice;
 
*VE Nosice – VN Dolné Kočkovce;
 
*VE Trenčín – VN Trenčianske Biskupice;
 
*VE Horná Streda VN Sĺňava.
 
  
[[Súbor:rb_dp2.16.png|framed|center|Obr. 2.16  Pozdĺžny rez odpadného kanála od VE do vodnej nádrže]]
+
=Literatúra=
 +
# Mäsiar, E., Kamenský, J. Hydraulika pre stavebných inžinierov II, Bratislava, ALFA 1989.
 +
# Dušička, P., Šulek, P., Popis algoritmov hydromodelovania navrhnutých pre SW model prípravy prevádzky VE, Bratislava, Technická dokumentácia 2006.
 +
# Walkenbach, J. Excel 2003 – Programování ve VBA, Brno, Computer Press 2006.
 +
# Weber, M., Breden, M. Veľká kniha řešení Excel VBA, Brno, Computer Press 2007.
 +
# Benický J.: Manipulačný poriadok pre vodné stavby Orava a Tvrdošín. Bratislava, Výskumný ústav vodného hospodárstva 2003.
 +
# SVP. š.p.: Manipulačný poriadok pre vodné stavby Liptovská Mara a Bešeňová. Ružomberok, SVP. š.p. Odštepný závod Piešťany, závod Povodie horného Váhu 2004.
 +
# SVP. š.p.: Manipulačný poriadok pre vodnú stavbu Krpeľany - Súčany - Lipovec. Ružomberok, SVP. š.p. Odštepný závod Piešťany, závod Povodie horného Váhu 2004. Aktualizácia vydania z roku 1989 vypracovaného Hydroconsultom Bratislava.
 +
# Cábel J.: Vodné dielo Žilina. Dočasný manipulačný poriadok IV. časť. Bratislava, VVB š.p. 1998.
 +
# SVP. š.p.: Manipulačný poriadok pre vodné dielo Hričov - Mikšová - Považská Bystrica. Púchov, SVP. š.p. Odštepný závod Piešťany, závod Púchov 1999. Aktualizácia vydania z roku 1987 vypracovaného Hydroconsultom Bratislava.
 +
# SVP. š.p.: Manipulačný poriadok pre vodné dielo Nosice. Púchov, SVP. š.p. Odštepný závod Piešťany, závod Púchov 2003.
 +
# SVP. š.p.: Manipulačný poriadok pre kaskádu Dolné Kočkovce - Ladce - Ilava - Dubnica - Trenčín. Púchov, SVP. š.p. Odštepný závod Piešťany, závod Púchov 2002.
 +
# SVP. š.p.: Manipulačný poriadok pre vodnú stavbu II Važská kaskáda Trenčianske Biskupice - Kostolná - Nové Mesto - Horné Streda. Piešťany, SVP. š.p. Odštepný závod Piešťany 2004.
 +
# Benický J. a kol.: Manipulačný poriadok pre vodné stavby Drahovce - Madunice. Bratislava, Hydroconsult Bratislava 1989. Prepracoval: SVP. š.p. Odštepný závod Piešťany a závod Povodia stredného Váhu II, Piešťany 2005.
 +
# SVP. š.p.: Manipulačný poriadok pre vodné dielo Kráľová. Šaľa, SVP. š.p. Odštepný závod Šaľa 1997.
  
2.5.6  Úsek VE - rieka
+
=Príloha A=
Úsek VE – rieka môžeme rozdeliť na úseky s odpadným kanálom a bez odpadného kanála od objektu VE. V prípade s odpadným kanálom bude nábehová hladina vypočítaná funkciou, kde vstupný parameter na určenie nábehovej hladiny bude prietok  z VE. Úsek, kde sa kanál nenachádza bude DH vypočítaná priamo funkciou  .
+
Tabuľka hodnôt rýchlostného súčiniteľa ''C'' k parametrom&nbsp;''R'' a&nbsp;''n'' <nowiki>[</nowiki>1<nowiki>]</nowiki>
  
Medzi úseky s kanálom patria:
+
{|border="2" cellspacing="0" cellpadding="4" width="100%"
*VE Tvrdošín – rieka Orava;
+
|      n <nowiki>[</nowiki>-<nowiki>]</nowiki><br>R <nowiki>[</nowiki>m<nowiki>]</nowiki>
*VE Bešeňová – rieka Váh;
+
|align = "right"|'''0,014'''
*VE Lipovec – rieka Váh;
+
|align = "right"|'''0,017'''
*VE Madunice – rieka Váh.
+
|align = "right"|'''0,020'''
Medzi úseky bez kanála patria:
+
|align = "right"|'''0,025'''
*VE Kráľová – rieka Váh.
+
|align = "right"|'''0,028'''
 +
|align = "right"|'''0,030'''
 +
|align = "right"|'''0,035'''
 +
|align = "right"|'''0,040'''
  
[[Súbor:rb_dp2.17.png|framed|center|Obr. 2.17 Pozdĺžny rez odpadného kanála od VE do koryta rieky]]
+
|-
 +
|align = "right"|'''0,20'''
 +
|align = "right"|55,21
 +
|align = "right"|43,35
 +
|align = "right"|35,65
 +
|align = "right"|26,76
 +
|align = "right"|23,56
 +
|align = "right"|21,23
 +
|align = "right"|17,34
 +
|align = "right"|14,48
  
==Vyhodnotenie parametrov pre kanál==
+
|-
Podľa konštrukčných parametrov derivačných kanálov vyhľadaných v technickej dokumentácii pre jednotlivé stupne vážskej kaskády bolo možné vytvoriť základy hydraulického modelu. Aby sme boli schopný použiť rovnicu na výpočet hydraulických strát na jednotlivých úsekoch kanálov nevyhnutnou podmienkou je:
+
|align = "right"|'''0,40'''
*Vytvoriť náhradu funkcie rýchlostného súčiniteľa (Chézyho súčiniteľ);
+
|align = "right"|61,72
*Vyhodnotiť drsnosť kanálov z nameraných hodnôt HH a DH pri rôznych prietokových pomeroch.
+
|align = "right"|49,41
 +
|align = "right"|41,25
 +
|align = "right"|31,80
 +
|align = "right"|28,40
 +
|align = "right"|25,80
 +
|align = "right"|21,51
 +
|align = "right"|18,30
  
===Rýchlostný súčiniteľ===
+
|-
Rýchlostný súčiniteľ je hodnota vystupujúca v Chézyho vzťahu na výpočet strednej prierezovej rýchlosti v alebo Q. V našom prípade nastáva malý problém a to
+
|align = "right"|'''0,60'''
*nemôžeme použiť Pavlovského empirický vzťah 1.15 z dôvodu platnosti hraníc 0,1m < R < 3,0m, kde v mnohých prípadoch hydraulický polomer kanálov dosahuje R > 5;
+
|align = "right"|65,86
*máme k dispozícii diskrétne hodnoty v rozsahu 0,1m až 5,0m, ale my potrebujeme vyhodnotiť rýchlostný súčiniteľ aj nad uvedený parameter hydraulického polomeru.
+
|align = "right"|53,41
 +
|align = "right"|44,90
 +
|align = "right"|35,20
 +
|align = "right"|31,67
 +
|align = "right"|28,90
 +
|align = "right"|24,40
 +
|align = "right"|21,03
  
Tento problém vieme riešiť preložením bodov vhodne opisujúcou polynómickou funkciou, ktorá z najmenšou odchýlkou bude vyhovovať našim požiadavkám. Pretože v tejto chvíli nevieme z určitosťou povedať, či pri zmene pozície nábehovej hladiny smerom k maximálnej prevádzkovej hladine zväčšime hydraulický polomer kanála natoľko, že funkcia rýchlostného súčiniteľa sa začne vyhodnocovať nekorektne, musíme zabezpečiť jej plynulý prechod lineárnou interpolačnou funkciu cez posledné dve hodnoty. Teraz aj keď nevieme povedať aké skutočné hodnoty rýchlostného súčiniteľa sú nad R > 5 zabezpečíme korektné správanie celého hydraulického modelu (Obr. 2.18). Pre potreby hydromodelu budeme pracovať so stupňom drsnosti n = {0,014; 0,017; 0,020; 0,025; 0,030; 0,035; 0,040}, pre ktoré vytvoríme funkcie polynómu.
+
|-
 +
|align = "right"|'''0,80'''
 +
|align = "right"|68,93
 +
|align = "right"|56,35
 +
|align = "right"|47,70
 +
|align = "right"|37,84
 +
|align = "right"|34,25
 +
|align = "right"|31,30
 +
|align = "right"|26,66
 +
|align = "right"|23,18
  
[[Súbor:rb_dp2.18.png|framed|center|Obr. 2.18  Príklad funkcie opisujúcu rýchlostný súčiniteľ z hydraulického polomeru kanála]]
+
|-
 +
|align = "right"|'''1,00'''
 +
|align = "right"|71,43
 +
|align = "right"|58,82
 +
|align = "right"|50,00
 +
|align = "right"|40,00
 +
|align = "right"|36,36
 +
|align = "right"|33,33
 +
|align = "right"|28,57
 +
|align = "right"|25,00
  
===Drsnosť kanálov===
+
|-
Vyhodnotenie drsnosti jednotlivých kanálov na vážskej kaskáde musí prebehnúť z rôznych vzoriek parametrov hladín a prietokov, ktoré ďalej vieme analyzovať v predpripravenom hydraulickom modely. Princíp vyhľadávania drsností kanálov je v nájdení najvhodnejšej kombinácii drsností prívodného a odpadného kanála. Pomocným nástrojom, ktorý bude na tento účel naprogramovaný v kanáloch hydromodelu nastaví nábehovú hladinu ako HHVE a cez prietok bude hľadať najbližšiu pozíciu DHVE. K dispozícii bude n vzoriek a pre každú urobí kompletné vyhodnotenie. Z najväčšieho počtu súhlasných kombinácii drsností všetkých vzoriek tieto použijeme do technických parametrov pre jednotlivé kanále.
+
|align = "right"|'''1,20'''
 +
|align = "right"|73,29
 +
|align = "right"|60,47
 +
|align = "right"|51,50
 +
|align = "right"|41,40
 +
|align = "right"|37,71
 +
|align = "right"|34,63
 +
|align = "right"|29,79
 +
|align = "right"|26,18
  
==Podporné metódy a funkcie pre hydromodel==
+
|-
Všetky doterajšie riešenia sa týkali úsekov derivačných kanálov na stupňoch vážskej kaskády bez ich priameho hydraulického prepojenia. V tejto chvíli vieme matematicky opísať hydraulický model od nádrže po poslednú VE každého stupňa. Podmienka dobrého návrhu hydromodelu je celkové hydraulické prepojenie stupňov vážskej kaskády. Každý stupeň má poslednú VE (v niektorých prípadoch je zároveň prvá), za ktorou je odpadný kanál ústiaci do vodnej nádrže alebo koryta rieky.  Práve tieto úseky nás budú zaujímať a v návrhu modelu je potrebné rozpracovať nasledujúce:
+
|align = "right"|'''1,40'''
*Navrhnúť hydraulické prepojenie odpadného kanála s vodnou nádržou nasledujúceho stupňa;
+
|align = "right"|74,86
*Vytvoriť funkciu, ktorá bude opisovať pozíciu hladiny v rieke na konci odpadného kanála ústiaceho do rieky;
+
|align = "right"|61,88
*Vytvoriť funkciu, ktorá bude opisovať pozíciu hladiny v rieke pri jalových vodách cez hať pod poslednou VE daného stupňa.
+
|align = "right"|52,75
 +
|align = "right"|42,64
 +
|align = "right"|38,91
 +
|align = "right"|35,76
 +
|align = "right"|30,86
 +
|align = "right"|27,20
  
===Hydraulické prepojenie odpadného kanála s vodnou nádržou nasledujúceho stupňa===
+
|-
Ako bolo spomenuté pri úseku VE – nádrž jedná sa o miesto, ktoré nie je možné spojiť priamo z nádržou spôsobom nábehová hladina = hladina nádrže. Budeme hľadať riešenie ako sa bude teoreticky správať hladina nádrže na konci odpadného kanála pri rôznych prietokoch a hladinách v nádrži. Konštrukčné parametre kanála sú známe a vieme, že rôzne prietoky potrebujú rôzne minimálne prietokové plochy nato, aby sa cez kanál dostali do nádrže. Preto budeme hľadať najmenšiu prietokovú plochu prostredníctvom potenciálu najmenšieho vzdutia DH, ktorý začneme vyhodnocovať cez zmeny pozície nábehovej hladiny kanála. Samozrejme ak hladina nádrže svojou pozíciou zabezpečí väčšiu prietokovú plochu ako je potrebné pre daný prietok, posun nábehovej hladiny sa zastaví na hladine nádrže. Princíp metódy bude nasledovný. Pozícia nábehovej hladiny na začiatku bude v dostatočnej výške a to optimálne nad maximálnou DH pod VE. Metóda hľadania najmenšieho potenciálu DH spočíva v znižovaní nábehovej hladiny dokedy klesajúca DH nenaberie stúpajúcu tendenciu. V tomto bode sa hľadanie ukončí a máme tak vytvorené hydraulické spojenie jednotlivých stupňov vážskej kaskády.
+
|align = "right"|'''1,60'''
 +
|align = "right"|76,29
 +
|align = "right"|63,11
 +
|align = "right"|53,90
 +
|align = "right"|43,72
 +
|align = "right"|39,96
 +
|align = "right"|36,80
 +
|align = "right"|31,83
 +
|align = "right"|28,13
  
===Funkcia opisujúca pozíciu hladiny v rieke na konci odpadného kanála ústiaceho do rieky===
+
|-
Podobne ako je to pre spojenie odpadného kanála s nádržou aj v tomto prípade musíme nepriamo vyhodnotiť hladiny rieky cez predpripravený hydraulický model pre daný úsek Váhu. Opäť máme k dispozícii vzorku hladín v závislosti od prietoku cez VE. V tomto prípade nám pomocný výpočtový nástroj bude hľadať cez potenciál DH hladinu rieky pri danom prietoku. Vyhodnotené budú vzorky hladín cez stanovené drsnosti a konečný výber hladín rieky bude preložený vhodnou polynómickou funkciou.
+
|align = "right"|'''1,80'''
 +
|align = "right"|77,57
 +
|align = "right"|64,23
 +
|align = "right"|54,95
 +
|align = "right"|44,72
 +
|align = "right"|40,91
 +
|align = "right"|37,70
 +
|align = "right"|32,71
 +
|align = "right"|28,95
  
===Funkcia opisujúca pozíciu hladiny v rieke pri jalových vodách cez hať pod poslednou VE daného stupňa===
+
|-
Jalové vody, ktoré sa opakujú raz za 100 rokov budeme z časových dôvodov vyhodnocovať z meraní uvedených v technickej dokumentácii pre jednotlivé stupne vážskej kaskády. Vzhľadom k tomu, že vo väčšine prípadov sa udáva hladina 100 ročných vôd pod VE, závislosť funkcie na objektoch TV, ZI, HS, MA bude podľa známych bodov lineárna. Pri objektoch HR, TN bude závislosť tiež lineárna, ale vznikne nová dynamická lineárna funkcia spojením dvoch lineárnych funkcii, z ktorých jedna opísuje závislosť minimálnej a druhá maximálnej hladiny v nádrži od prietoku jalovej vody.
+
|align = "right"|'''2,00'''
Príklad spojenia týchto funkcii priblížime pre VD NO.
+
|align = "right"|78,72
 +
|align = "right"|65,29
 +
|align = "right"|55,85
 +
|align = "right"|45,64
 +
|align = "right"|41,78
 +
|align = "right"|38,56
 +
|align = "right"|33,51
 +
|align = "right"|29,73
  
Funkcia pre max. hladinu nádrže: <math>y=0,000820513x+279,19</math>
+
|-
Funkcia pre min. hladinu nádrže: <math>y=0,00138462x+277,808</math>
+
|align = "right"|'''2,20'''
 +
|align = "right"|79,79
 +
|align = "right"|66,23
 +
|align = "right"|56,70
 +
|align = "right"|46,48
 +
|align = "right"|42,58
 +
|align = "right"|39,33
 +
|align = "right"|34,26
 +
|align = "right"|30,45
  
kde <math>x={{Q}_{Celkov\acute{y}\ Jalov\acute{y}}}</math>
+
|-
 +
|align = "right"|'''2,40'''
 +
|align = "right"|80,72
 +
|align = "right"|67,05
 +
|align = "right"|57,50
 +
|align = "right"|47,24
 +
|align = "right"|43,30
 +
|align = "right"|40,06
 +
|align = "right"|34,94
 +
|align = "right"|31,13
  
Výpočet koeficienta <math>x^1</math>:
+
|-
 +
|align = "right"|'''2,60'''
 +
|align = "right"|81,64
 +
|align = "right"|67,88
 +
|align = "right"|58,25
 +
|align = "right"|47,96
 +
|align = "right"|44,03
 +
|align = "right"|40,73
 +
|align = "right"|35,60
 +
|align = "right"|31,75
  
<math>\begin{align}
+
|-
  & Pomer\_koe{{f}^{1}}=\frac{0,000820513-0,00138462}{279,6-274,6}=-0,0001128214 \\
+
|align = "right"|'''2,80'''
& {{k}^{1}}=0,000820513-Pomer\_koe{{f}^{1}}\left( 279,6-274,6 \right) \\
+
|align = "right"|82,50
\end{align}</math>
+
|align = "right"|68,64
 +
|align = "right"|58,95
 +
|align = "right"|48,64
 +
|align = "right"|44,69
 +
|align = "right"|41,46
 +
|align = "right"|36,20
 +
|align = "right"|32,35
  
 +
|-
 +
|align = "right"|'''3,00'''
 +
|align = "right"|83,29
 +
|align = "right"|69,35
 +
|align = "right"|59,60
 +
|align = "right"|49,28
 +
|align = "right"|45,30
 +
|align = "right"|42,00
 +
|align = "right"|36,70
 +
|align = "right"|32,80
  
Výpočet koeficienta <math>x^0</math>:
+
|-
 +
|align = "right"|'''3,20'''
 +
|align = "right"|84,07
 +
|align = "right"|70,05
 +
|align = "right"|60,25
 +
|align = "right"|49,88
 +
|align = "right"|45,89
 +
|align = "right"|42,56
 +
|align = "right"|37,00
 +
|align = "right"|33,00
  
<math>\begin{align}
+
|-
  & Pomer\_koe{{f}^{0}}=\frac{279,19-277,808}{279,6-274,6}=0,2764 \\
+
|align = "right"|'''3,40'''
& {{k}^{0}}=279,19-Pomer\_koe{{f}^{0}}\left( 279,6-274,6 \right) \\
+
|align = "right"|84,72
\end{align}</math>
+
|align = "right"|70,64
 +
|align = "right"|60,80
 +
|align = "right"|50,40
 +
|align = "right"|46,10
 +
|align = "right"|42,80
 +
|align = "right"|37,30
 +
|align = "right"|33,20
  
Funkcia popisujúca dynamické závislosti:
+
|-
 +
|align = "right"|'''3,60'''
 +
|align = "right"|85,43
 +
|align = "right"|71,29
 +
|align = "right"|61,35
 +
|align = "right"|50,70
 +
|align = "right"|46,50
 +
|align = "right"|43,15
 +
|align = "right"|37,60
 +
|align = "right"|33,45
  
<math>\begin{align}
+
|-
  & H{{H}_{Koniec\ OK}}={{k}^{1}}x+{{k}^{0}} \\
+
|align = "right"|'''3,80'''
& H{{H}_{Koniec\ OK}}=\left[ 0,000820513+0,0001128214\cdot \left( 279,6-H{{H}_{N\acute{a}dr\check{z}\ Aktu\acute{a}\ln a}} \right) \right]\cdot {{Q}_{Celkov\acute{y}\ Jalov\acute{y}}}+ \\
+
|align = "right"|86,07
& +279,19-0,2764\cdot \left( 279,6-H{{H}_{N\acute{a}dr\check{z}\ Aktu\acute{a}\ln a}} \right) \\
+
|align = "right"|71,88
\end{align}</math>
+
|align = "right"|61,90
 +
|align = "right"|51,00
 +
|align = "right"|46,90
 +
|align = "right"|43,50
 +
|align = "right"|37,90
 +
|align = "right"|33,70
  
[[Súbor:rb_dp2.19.png|framed|center|Obr. 2.19  Grafické spojenie lineárnych funkcii do dynamickej funkcie]]
+
|-
 +
|align = "right"|'''4,00'''
 +
|align = "right"|86,72
 +
|align = "right"|72,41
 +
|align = "right"|62,40
 +
|align = "right"|51,30
 +
|align = "right"|47,20
 +
|align = "right"|43,75
 +
|align = "right"|38,10
 +
|align = "right"|33,90
  
Nakoniec závislosť hladín z jalovej vody cez hať na objektoch BE, KR, NO, KA bude opísaná polynómom n-tého stupňa.
+
|-
 +
|align = "right"|'''4,20'''
 +
|align = "right"|87,29
 +
|align = "right"|72,94
 +
|align = "right"|62,90
 +
|align = "right"|51,55
 +
|align = "right"|47,50
 +
|align = "right"|44,00
 +
|align = "right"|38,30
 +
|align = "right"|34,15
  
Všetky uvedené funkcie opisujú závislosť HH (nádrže, rieky) na konci odpadného kanála poslednej VE od prietoku nad hltnosť VE alebo derivačného kanála. Pokiaľ nastane obmedzenie užívateľom cez stupeň do hltnosti objektov (ostatnú vodu prevedie model cez hať), hladina na konci odpadného kanála bude vyhodnotená metódou alebo funkciou, akoby sa jednalo o prietok cez objekty VE. Hydromodel v tomto prípade vyhodnotí aj časový posun, ktorý je potrebný pri prietoku starým korytom rieky z hate do cieľa pod objekt VE.
+
|-
 +
|align = "right"|'''4,40'''
 +
|align = "right"|87,86
 +
|align = "right"|73,47
 +
|align = "right"|63,40
 +
|align = "right"|51,80
 +
|align = "right"|47,70
 +
|align = "right"|44,20
 +
|align = "right"|38,45
 +
|align = "right"|34,35
  
==Transformačný model==
+
|-
Aby sme sa dopracovali k hľadanému výsledku, ktorým je elektrický výkon jednotlivých objektov VE, výsledky hydrologického a hydraulického modelu musíme transformovať. Na základe prietoku, spádu, účinnosti a dostupnosti technologického zariadenia transformačný model vypočíta elektrický výkon s presným počtom nasadených turbogenerátorov v jednotlivých hodinách dňa. Ďalšou úlohou modelu bude pri nedostupnosti jedného TG vyhodnotiť nasadenie ďalšieho dostupného TG do prevádzky.
+
|align = "right"|'''4,60'''
 +
|align = "right"|88,43
 +
|align = "right"|73,94
 +
|align = "right"|63,85
 +
|align = "right"|52,00
 +
|align = "right"|47,90
 +
|align = "right"|44,50
 +
|align = "right"|38,60
 +
|align = "right"|34,45
  
Základnú rovnicu transformačného modelu môžeme zapísať v tvare:
+
|-
 +
|align = "right"|'''4,80'''
 +
|align = "right"|89,00
 +
|align = "right"|74,41
 +
|align = "right"|64,25
 +
|align = "right"|52,25
 +
|align = "right"|48,10
 +
|align = "right"|44,70
 +
|align = "right"|38,70
 +
|align = "right"|34,55
  
{{vzorec|<math>{{P}_{(krok)}}=\frac{{{Q}_{(krok)}}\cdot g\cdot {{h}_{(krok)}}\cdot {{\eta }_{TG}}}{1000}</math>|2.6}}
+
|-
 +
|align = "right"|'''5,00'''
 +
|align = "right"|89,50
 +
|align = "right"|74,88
 +
|align = "right"|64,70
 +
|align = "right"|52,50
 +
|align = "right"|48,30
 +
|align = "right"|44,90
 +
|align = "right"|38,80
 +
|align = "right"|34,65
  
kde:
 
{|
 
|krok
 
|krok výpočtu v jednej hodine;
 
|-
 
|P
 
|Elektrický výkon [MW];
 
|-
 
|Q
 
|Prietok vody cez TG [<math>m^3.s^{-1}</math>];
 
|-
 
|g
 
|konštanta gravitačného zrýchlenia [<math>m.s^{-2}</math>];
 
|-
 
|h
 
|spád (h = HH – DH) [m];
 
|-
 
|ηTG
 
|celková účinnosť (turbína + generátor) [%].
 
 
|}
 
|}
  
===Inverzná transformácia elektrického výkonu===
 
Inverzná transformácia bude slúžiť na vyjadrenie prietoku z regulačného výkonu (SRV) z prvej VE na derivačnom kanály. Aby sme dosiahli vyrovnanú bilanciu vody pri pôsobení regulačného signálu na objekty VE pri minimálnej a maximálnej požiadavke na SRV je potrebné, aby v rovnakom čase boli navrhnuté aj prietoky cez VE pri zmenách rovnako. Pomocou výkonu pre SRV bude vyhodnotený prietok na SRV prvej VE a pomocou spätnej transformácie môžeme počítať rezervu SRV pri daných spádových a prietokových pomeroch na ostatných VE príslušného derivačného kanála. Transformačný model bude aj na základe týchto výpočtov vyhodnocovať prevádzkové spustenie, prípadne odstavenie TG. Ak kapacita dosiahnuteľného výkonu bude vyčerpaná (niektorý z TG bude z technických príčin mimo prevádzky), vodu ktorú hydrologický a hydraulický model vyhodnotí na danom stupni vážskej kaskády prevedie transformačný model na jalovo cez kanál.
 
Základnú inverznú rovnicu pre vyjadrenie prietoku z výkonu SRV môžeme zapísať v tvare:
 
  
{{vzorec|<math>{{Q}_{SRV(krok)}}=\frac{{{P}_{SRV(krok)}}\cdot 1000}{g\cdot {{h}_{(krok)}}\cdot {{\eta }_{TG}}}</math>|2.7}}
+
=Príloha B=
 +
Súčasťou tejto diplomovej práce je CD médium s&nbsp;nasledujúcim obsahom:
 +
 
 +
1. Dokument tejto práce vo formáte DOC;
  
kde:
+
2. Výber z&nbsp;plnej verzie hydromodelu pre SVK OR-TV, LM-BE, KR-SU-LI.
{|
 
|-
 
|PSRV
 
|Elektrický výkon sekundárnej regulácie výkonu [MW];
 
|-
 
|QSRV
 
|Prietok vody pre sekundárnu reguláciu výkonu cez TG [<math>m^3.s^{-1}</math>];
 
|}
 
Ostatné členy rovnice sú totožné s členmi transformačnej rovnice.
 
  
=Záver=
+
* Hydromodel.xls
=Literatúra=
+
* HMVN_Váh.xls
 +
* HT_modul.xls
 +
* Dotok_vody_modul.xls
 +
* MLP_Qmax_Váh.xls
 +
* MLP_Qmax_kanál_Váh.xls
 +
* MLP_Qmin_kanál_Váh.xls

Aktuálna revízia z 14:58, 9. december 2010

Modelovanie a simulácia je jednou z hlavných metód ako overovať chovanie analyzovaných systémov. Matematický model je abstrakcia pretransformovaná do počítačového programu, ktorý je podobný pôvodnému objektu, na ktorom vykonávame experimenty. Modelový experiment je riešenie vhodne formulovanej matematickej úlohy s rešpektovaním vzťahov opisujúcich model pre konkrétne podmienky. Komunikácia užívateľa s programom prebieha cez príkazy a cez zadanie počiatočných a koncových podmienok konkrétnych parametrov. Výsledky sú z dôvodu prehľadnosti zobrazené v grafickej aj číselnej forme. Optimalizovaným modelovaním hľadáme také riešenie úlohy, ktoré spĺňa naše predpoklady výsledku, prípadne nám umožní sa vyhnúť problémovým situáciám v modelovanom systéme. Našou úlohou, ktorú budeme ďalej rozoberať, je optimalizované modelovanie vážskej kaskády. V druhej kapitole boli uvedené všetky časti a moduly súvisiace so samotným hydromodelovaním. Postupnosť krokov modelovania je nasledovná:

  • Vstup pre zadanie počiatočných a koncových podmienok;
  • Hydromodelovanie;
  • Analýza výsledkov hydromodelovania.

Vstup pre zadanie počiatočných a koncových podmienok

Hydromodel je komplexný nástroj pozostávajúci z jednotlivých modulov, ktoré sú naprogramované v prostredí VBA. Práca s nástrojom si však vyžaduje hlbšie znalosti v danej problematike. Pre zadávanie vstupov, počiatočných a koncových podmienok, nám slúži súbor s programom Hydromodel.xls, prostredníctvom ktorého je možné editovať vstupy pre celý model. Vstupy sú kontrolované proti zadaniu nekorektných údajov a zároveň v prípade, že zadaná hodnota presahuje obmedzujúce parametre, model príslušné vstupy vyznačí (podfarbí pozadie a zmení farbu písma) a zablokuje spustenie hydromodelovania. Takéto zadávanie vstupov môžeme nazvať - zadávanie vstupov naslepo.


Pri spustení (otvorení súboru Hydromodel.xls) sa zobrazí priestor s výsledkami tvorený číselným zobrazením nasadeného výkonu jednotlivých objektov VE, ktorých výkony sumarizuje skupina SVK. Nakoniec súčet všetkých skupín je vyjadrený v položke DDZ Váh a graficky zobrazený spolu s položkou Ocenenie DDZ (Obr. 3.1). Suma ocenenia DDZ vyjadruje prácu a veľkosti hodnôt ocenenia vyjadrujú pomer nasadenia v priebehu celého dňa.

Obr. 3.1 Ovládací panel

V súčasnej dobe nástroj nedisponuje databázou vzhľadom k rozsahu aktuálne riešenej úlohy. Rozšírenie bude vyžadovať analýzu a následne tvorbu vhodne nastaveného manažmentu dát, ktorý zabezpečí rozumné riešenie pre použitie rozšírenia na týždennú, mesačnú a ročnú prípravu prevádzky a väzbu na riadiaci systém. V súčasnej dobe databázu nahrádza interná pamäť, ktorá zabezpečuje potrebný manažment údajov v rámci D-1, D a D+1. Plánovanie nasadenia výkonu počas dňa pomocou modelu vieme realizovať v rozsahu 23 až 25 hodín.


Do jednotlivých medzipovodí sú prítoky vstupom, ktorý kryje časť výroby. Model sám nedokáže na požadovanú výrobu samočinne nastaviť odtoky z vyrovnávacích nádrží TV a BE tak, aby kryl ostatnú časť výroby. Na tento účel je modul, súčasť hydromodelu, ktorý umožňuje aproximovať výrobu s presnosťou, ktorú rozoberieme neskôr. Modul sprehľadňuje zadávanie vybraných vstupov a v okamihu zadania dáva odozvu v podobe odhadu výroby. Modul sa volá zo súboru Hydromodel.xls. Pri tomto volaní si odovzdajú informácie o aktuálnych pomeroch prietoku a výkonu. V prípade veľkých vôd, keď dochádza k odstaveniu skupín SVK a prietokové, výkonové pomery sa rovnajú nule, modul s názvom Aproximácia_výroby_Váh.xls, ktorý aproximuje výrobu, použije prednastavené technické parametre TG jednotlivých objektov VE. Prostredie modulu je na obrázku 3.2. Pri odovzdávaní informácii o pomeroch modul aproximácie výroby preberie z riadiacej časti aj prítoky za D-1 z jednotlivých medzipovodí a v rámci výpočtu prostredníctvom aktuálnych prítokov z D-1 a priemerných hodnôt prietoku cez skupiny SVK počas dňa D stanový najlepší odhad výroby pre hydromodel. Tak ako riadiaca časť programu aj modul aproximácie používa časť hydraulickej časti pre vodné nádrže a modul pre dotok vody. Modul aproximácie integruje prvky prehľadného nastavenia celého riadiaceho súboru Hydromodel.xls.

Obr. 3.2 Panel aproximácie výroby

Jednotlivé položky sú:

  1. Minimálny počet nasadenia TG skupiny SVK do prevádzky;
  2. Maximálny počet nasadenia TG skupiny SVK do prevádzky;
  3. Maximálny prietok cez skupiny SVK;
  4. Hodnotu počiatočnej hladiny VN;
  5. Hodnotu plánovanej koncovej hladiny VN;
  6. Prítoky vody do jednotlivých medzipovodí;
  7. Biologické odtoky z medzipovodí do koryta;
  8. Priemerný denný prietok cez SVK;
  9. Stanovenie špičkovej alebo priebežnej prevádzky skupín SVK;
  10. CheckBox tlačidlá stanovujeme prepúšťanie Qj do hltnosti kanálov cez kanál alebo hať;
  11. Prepínač umožňujúci prenášať obmedzenia oblastí 1, 2, 3;.
  12. CheckBox povolenia priebežného výpočtu po zadaní hodnôt v oblastí 4, 5, 6, 7.


Celý nástroj zabezpečuje integritu dát, nekorektné vstupy výstražne označí a v prípade odoslania (napr. záporný priemerný prietok SVK) blokuje odoslanie. Samotná presnosť aproximácie je závislá od oblasti prietokov, v ktorej sa modeluje. Pri prevádzke do hltnosti kanálov a mierne nad túto hltnosť je odhad uspokojivý. V oblasti prietokov veľkých vôd, kde nastáva aj odstavenie SVK z prevádzky a pomerové hodnoty prietoku a výkonu sú načítané z technických parametrov TG objektov VE, je odhad slabý. V tomto prípade sa využíva prednosť v jednoduchosti zadávania údajov, nie presnosť odhadu.


Príklad presností odhadov: Zaujímavosť zistená pri analýze odpúšťania vody z VN TV na krytie výkonu Váhu je v tom, že plné výkonové nasýtenie Váh dosiahne približne za štyri dni. Hladiny sú pritom modelované tak, aby počiatočná sa rovnala koncovej hladine. Dotoky z dňa D model po ukončení modelovania SVK pripraví pre D+1. V modelovom príklade nebudeme uvažovať jalovú vodu starým korytom. Všetky medzipovodia majú nastavený bočný prítok rovnajúci sa biologickému odtoku a odtok z vyrovnávacích nádrží TV, BE je nastavený na minimálny odtok.

Tab. 3.1 Výkonová nasýtenosť VE Váh

P.č. Aproximácia
[MWh]
Hydromodel
[MWh]
Rozdiel
[MWh]
Čas spracovania úlohy
[sek.]
1. 1517,6 1693,7 176,1 16,55
2. 3280,4 3383,7 103,3 25,14
3. 4699,2 4827,9 128,7 17,94
4. 4928,4 5066,7 138,3 18,88
5. 4946,4 5110,8 164,4 20,59


V tabuľke 3.1 prvý bod je počiatočný stav Váhu s minimálnymi prítokmi do medzipovodí a odtokmi z vyrovnávacích nádrží. V druhom bode je zmenený iba odtok z VN TV z 4 m3.s-1 skokovo na 65 m3.s-1. Je vidieť postupné výkonové nasycovanie na Váhu a plné nasýtenie nastáva na štvrtý deň. Pri aproximácii výroby je v tomto príklade rozdiel v rozsahu od 100 MWh do 180 MWh, čo môžeme považovať za uspokojivý výsledok. Chyba je spôsobená spriemerovaním odtokov z SVK počas dňa, pričom hydromodel nasadzuje prietok do špičiek a väčšina objemu ostáva v dni D a nie je presunutá do dňa D+1. Pokiaľ by modul aproximácie nevyužíval model dotokov vody, hodnota rozdielu odhadu voči hodnote z hydromodelu by mohla presiahnuť 1500 MWh v čase väčších zmien prítokov, čo je neprípustné!


Odľahčenie výroby namodelujeme analogicky. VN TV nastavíme pôvodný odtok 4 m3.s-1. Tabuľka 3.2 opisuje tak, ako v prípade nasýtenia výkonu v prvom bode, počiatočný stav Váhu a v nasledujúcom bode zníženie odtoku z VN TV. Rozdiel odhadu je pri odľahčovaní v rozsahu od 160 MWh do 220 MWh, čo môžeme opäť považovať za uspokojivé. Opäť dostávame štyri dni do úplného odľahčenia výroby VE.

Tab. 3.2 Odľahčenie výroby VE Váh

P.č. Aproximácia
[MWh]
Hydromodel
[MWh]
Rozdiel
[MWh]
Čas spracovania úlohy
[sek.]
1. 4946,4 5110,8 164,4 20,59
2. 3195,2 3412,9 217,7 29,38
3. 1828,5 2031,7 203,2 20,41
4. 1526,1 1709,3 183,2 17,41
5. 1514,6 1693,7 179,1 16,98


Časy boli namerané na PC s procesorom Intel® Core™2 CPU 6400 @ 2,13 GHz, 1,99 GB RAM, OS Microsoft Windows XP Professional a Office 2003.

Hydromodelovanie

Po úspešnom zápise vstupov do modelu, priamo alebo pomocou aproximácie výroby, máme model pripravený na spustenie. V tejto časti popíšeme jednotlivé interakcie medzi modulmi hydromodelu. Po spustení nasleduje načítanie vstupných dát optimalizačným modelom, ktorý na základe vstupov z kapitoly 2.9 navrhne nasadenie prietoku cez aktuálne počítanú kaskádu. SVK sú postupne spracovávané v smere toku vody tak, aby bola zabezpečená hydrologická a hydraulická väzba.


Poradie modelovania stupňov vážskej kaskády:

  1. SVK OR-TV
  2. SVK LM-BE
  3. SVK KR-SU-LI
  4. SVK ZI
  5. SVK HR-MI-PB
  6. SVK NO
  7. SVK LA-IL-DU-TN
  8. SVK KO-NM-HS
  9. SVK MA
  10. SVK KA


Prekročenie prevádzkových hladín optimalizačný model kontroluje prostredníctvom hydraulického modelu vodných nádrží. Prekročenie koriguje a po poslednej hodine dňa hľadá riešenie dosiahnutia plánovanej koncovej hladiny. Počas optimalizácie verifikuje maximálne nasadenie prietoku cez SVK a jednotlivo cez TG, z ktorých vyberá pre výpočet najmenšiu hodnotu maximálneho obmedzenia. Tieto kontroly prebiehajú aj zo strany minimálnych prietokov cez TG. Vo výpočte optimalizácie SVK prebieha priebežný výpočet pohybu nábehových hladín v derivačných kanáloch. V čase veľkých vôd optimalizačný model porovnáva obmedzenia SVK pred prítokom a po prítoku veľkej vody starým korytom. V prípade rozdielu celú optimalizáciu daného stupňa opakuje s tým, že časť vody nad obmedzením prevedie z kanála na hať. Celý proces sa opakuje do chvíle rovnosti obmedzení pred prítokom a po prítoku vody starým korytom pod poslednú VE danej kaskády. Tento stav vidieť na dialógu Priebeh spracovania výpočtu počas hydromodelovania (Obr. 3.3) v položke Hodiny (SVK bez zmeny a cyklicky sa opakujúce hodnoty tejto položky).

Obr. 3.3 Dialóg - priebeh spracovania

Po ukončení optimalizácie sú postupne údaje importované do hydraulicko-transformačného modulu, kde výpočet prebehne veľmi rýchlo. Výpočet nastaví hydrauliku v kanáloch a na prietokové a spádové pomery transformačný model stanový počet strojov a výkonové zaťaženie v danej hodine. Výsledky sú poslané späť do súboru Hydromodel.xls.

Analýza výsledkov hydromodelovania

Pri prehliadaní výsledkov a kontrole v súbore Hydromodel.xls prechádzame záložkami prostredníctvom tlačidiel na obrázku 3.4.

Obr. 3.4 Panel ovládania - hydromodel

Postupne budeme prechádzať prostredím pre zadávanie údajov. Vstupy, ktoré môžeme editovať na obrázku 3.5, sú podfarbené modrým pozadím. Ostatné údaje nie je možné meniť. Ovládanie celého programu je jednoduché a intuitívne. Pri prepínaní sa zobrazia grafické výstupy. Je možné prepnúť sa na dátovú časť tlačidlom alebo sa prerolovať myšou. Na každej záložke pre SVK je v ľavom rohu príznak spustenia modelu. Pokiaľ je tento príznak v stave Stop model, na záložke je problém nekorektných údajov, ktoré treba opraviť.

Obr. 3.5 Panel zadávania vstupov a prehliadanie výstupov - hydromodel

Záver

V práci sme rozobrali problematiku hydromodelovania vážskej kaskády, ktorej praktický výstup je nástroj na modelovanie. Postupne, ako sme prechádzali jednotlivé kapitoly, mali sme možnosť zoznámiť sa s danou problematikou navrhovania hydrologického, hydraulického, transformačného a optimalizačného modelu. Spolu jednotlivé časti tvoria systém umožňujúci veľmi náročné matematicko-logické výpočty. Rozobrané nové prístupy modelovania v tejto oblasti, oproti stávajúcemu modelu na HED, umožnili hydraulické spojenie jednotlivých skupín SVK, prácu s hladinami v derivačných kanáloch, lepšiu metódu transformácie dotoku vody, modelovanie veľkej vody od bežných prevádzkových prietokov až po prietoky storočných vôd, autonómnosť optimalizácie, integritu údajov a celého modelu. Veľký rozdiel v nastavovaní hydromodelu oproti pôvodnému na HED je v zadávaní vstupov. Pôvodný model umožňuje nastaviť vstupy odhadom, využívajúc chybu ľudského faktora – multinásobná korekcia vstupných dát s postupným približovaním dosiahnutia požadovanej výroby, kontrola nasadenia turbogenerátorov do prevádzky a v mnohých prípadoch nutnosť manuálne tvarovať DDZ. Ak by sme porovnali obidva hydromodely zistíme, že pôvodný je kompilovaný, ale časová náročnosť výpočtu je porovnateľná s časom výpočtu nástroja vytvoreného v aplikácii Excel pomocou interpretovacieho jazyka VBA. Pokiaľ by sme mali vyhodnotiť celkovú časovú náročnosť modelovania jednej úlohy, stávajúci model vyžaduje plné nasadenie užívateľa v rozsahu jednej hodiny prípadne dlhšie. Naproti tomu aktuálne vytvorený model je schopný poskytnúť výsledok úlohy v priebehu pár minút od začiatku modelovania. Ak by sme chceli riešiť úlohu veľkej vody, pôvodný model na HED nie je schopný riešiť takúto úlohu.

Literatúra

  1. Mäsiar, E., Kamenský, J. Hydraulika pre stavebných inžinierov II, Bratislava, ALFA 1989.
  2. Dušička, P., Šulek, P., Popis algoritmov hydromodelovania navrhnutých pre SW model prípravy prevádzky VE, Bratislava, Technická dokumentácia 2006.
  3. Walkenbach, J. Excel 2003 – Programování ve VBA, Brno, Computer Press 2006.
  4. Weber, M., Breden, M. Veľká kniha řešení Excel VBA, Brno, Computer Press 2007.
  5. Benický J.: Manipulačný poriadok pre vodné stavby Orava a Tvrdošín. Bratislava, Výskumný ústav vodného hospodárstva 2003.
  6. SVP. š.p.: Manipulačný poriadok pre vodné stavby Liptovská Mara a Bešeňová. Ružomberok, SVP. š.p. Odštepný závod Piešťany, závod Povodie horného Váhu 2004.
  7. SVP. š.p.: Manipulačný poriadok pre vodnú stavbu Krpeľany - Súčany - Lipovec. Ružomberok, SVP. š.p. Odštepný závod Piešťany, závod Povodie horného Váhu 2004. Aktualizácia vydania z roku 1989 vypracovaného Hydroconsultom Bratislava.
  8. Cábel J.: Vodné dielo Žilina. Dočasný manipulačný poriadok IV. časť. Bratislava, VVB š.p. 1998.
  9. SVP. š.p.: Manipulačný poriadok pre vodné dielo Hričov - Mikšová - Považská Bystrica. Púchov, SVP. š.p. Odštepný závod Piešťany, závod Púchov 1999. Aktualizácia vydania z roku 1987 vypracovaného Hydroconsultom Bratislava.
  10. SVP. š.p.: Manipulačný poriadok pre vodné dielo Nosice. Púchov, SVP. š.p. Odštepný závod Piešťany, závod Púchov 2003.
  11. SVP. š.p.: Manipulačný poriadok pre kaskádu Dolné Kočkovce - Ladce - Ilava - Dubnica - Trenčín. Púchov, SVP. š.p. Odštepný závod Piešťany, závod Púchov 2002.
  12. SVP. š.p.: Manipulačný poriadok pre vodnú stavbu II Važská kaskáda Trenčianske Biskupice - Kostolná - Nové Mesto - Horné Streda. Piešťany, SVP. š.p. Odštepný závod Piešťany 2004.
  13. Benický J. a kol.: Manipulačný poriadok pre vodné stavby Drahovce - Madunice. Bratislava, Hydroconsult Bratislava 1989. Prepracoval: SVP. š.p. Odštepný závod Piešťany a závod Povodia stredného Váhu II, Piešťany 2005.
  14. SVP. š.p.: Manipulačný poriadok pre vodné dielo Kráľová. Šaľa, SVP. š.p. Odštepný závod Šaľa 1997.

Príloha A

Tabuľka hodnôt rýchlostného súčiniteľa C k parametrom Rn [1]

n [-]
R [m]
0,014 0,017 0,020 0,025 0,028 0,030 0,035 0,040
0,20 55,21 43,35 35,65 26,76 23,56 21,23 17,34 14,48
0,40 61,72 49,41 41,25 31,80 28,40 25,80 21,51 18,30
0,60 65,86 53,41 44,90 35,20 31,67 28,90 24,40 21,03
0,80 68,93 56,35 47,70 37,84 34,25 31,30 26,66 23,18
1,00 71,43 58,82 50,00 40,00 36,36 33,33 28,57 25,00
1,20 73,29 60,47 51,50 41,40 37,71 34,63 29,79 26,18
1,40 74,86 61,88 52,75 42,64 38,91 35,76 30,86 27,20
1,60 76,29 63,11 53,90 43,72 39,96 36,80 31,83 28,13
1,80 77,57 64,23 54,95 44,72 40,91 37,70 32,71 28,95
2,00 78,72 65,29 55,85 45,64 41,78 38,56 33,51 29,73
2,20 79,79 66,23 56,70 46,48 42,58 39,33 34,26 30,45
2,40 80,72 67,05 57,50 47,24 43,30 40,06 34,94 31,13
2,60 81,64 67,88 58,25 47,96 44,03 40,73 35,60 31,75
2,80 82,50 68,64 58,95 48,64 44,69 41,46 36,20 32,35
3,00 83,29 69,35 59,60 49,28 45,30 42,00 36,70 32,80
3,20 84,07 70,05 60,25 49,88 45,89 42,56 37,00 33,00
3,40 84,72 70,64 60,80 50,40 46,10 42,80 37,30 33,20
3,60 85,43 71,29 61,35 50,70 46,50 43,15 37,60 33,45
3,80 86,07 71,88 61,90 51,00 46,90 43,50 37,90 33,70
4,00 86,72 72,41 62,40 51,30 47,20 43,75 38,10 33,90
4,20 87,29 72,94 62,90 51,55 47,50 44,00 38,30 34,15
4,40 87,86 73,47 63,40 51,80 47,70 44,20 38,45 34,35
4,60 88,43 73,94 63,85 52,00 47,90 44,50 38,60 34,45
4,80 89,00 74,41 64,25 52,25 48,10 44,70 38,70 34,55
5,00 89,50 74,88 64,70 52,50 48,30 44,90 38,80 34,65


Príloha B

Súčasťou tejto diplomovej práce je CD médium s nasledujúcim obsahom:

1. Dokument tejto práce vo formáte DOC;

2. Výber z plnej verzie hydromodelu pre SVK OR-TV, LM-BE, KR-SU-LI.

  • Hydromodel.xls
  • HMVN_Váh.xls
  • HT_modul.xls
  • Dotok_vody_modul.xls
  • MLP_Qmax_Váh.xls
  • MLP_Qmax_kanál_Váh.xls
  • MLP_Qmin_kanál_Váh.xls