Souhrn

Monitoring a údržba vodních děl vyžadují znalost jejich referenčního stavu charakterizujícího návrhové parametry. Referenční stav je odvozován primárně z dostupných projektových dokumentací, které jsou vhodně zvolenými metodami převedeny na digitální model terénu dna. Aktuální úroveň dna je stanovována pomocí sonarových prostředků. Rozdílovou analýzou jednotlivých datových sad, které budou výhledově ukládány do specializovaného informačního systému, bude možné provádět efektivní řízení monitoringu údržby.

Předložená práce na příkladech vybraných pilotních lokalit shrnuje problematické aspekty identifikované při georeferencování digitalizovaných historických situačních výkresů, které mohou negativně ovlivňovat korektní lokalizaci jednotlivých částí vodních děl a následnou analýzu morfologických změn. V pilotní lokalitě vodní nádrže je rovněž hodnocena přesnost měřicí technologie použité při zaměřování dna.

Úvod

Problematika sedimentů představuje jedno z významných témat státní politiky v oblasti vodního hospodářství. Správci vodních děl musí podle § 59 zákona č. 254/2001 Sb., o vodách [1], plnit povinnosti a dodržovat podmínky, za kterých bylo vodní dílo povoleno a uvedeno do provozu. Hlavní cíle státní politiky v oblasti vodního hospodářství vycházejí z Vodní rámcové směrnice 2000/60/ES [2] a z obnovené Strategie EU pro udržitelný rozvoj z roku 2006.

Na evropské úrovni je patrná aktivita směrem k problematice sedimentů především v souvislosti s vnitrozemskou plavbou, potažmo rozvojem a údržbou Transevropské dopravní sítě (TEN-T). Ve snaze o implementaci Integrovaného evropského akčního plánu pro vnitrozemskou vodní plavbu z roku 2006 vznikl dokument Good Practice Manual on Inland Waterway Maintenance [3]. Uvedená práce má charakter komplexního manuálu údržby vodních cest vycházející z příkladů dobré praxe. Pozornost je rovněž soustředěna na vliv sedimentů na kvalitu vod, a to zejména zásluhou tzv. SedNet (European Sediment Network) [4].

Obr. 1. Ukázka digitalizace vrstevnic na podkladu situačního výkresu části zátopy vodní nádrže Letovice
Fig. 1. An example of the contour vectorization on the basis of the situational plan of the Letovice Dam

Během dlouhého období existence vodních děl dochází k přirozeným i antropogenně podmíněným a akcelerovaným morfologickým změnám jejich dna, resp. koryta. V exponovaných částech může docházet k zahlubování, naopak v oblastech s nižší dynamikou proudění jsou unášené plaveniny deponovány. Tyto změny mají v konečném důsledku dopad na průtočnou kapacitu a schopnost vodního díla (VD) plnit funkce, pro které bylo vybudováno. Proto je potřebné průběžně monitorovat stav vodních děl a v případě potřeby přistoupit k potřebným opatřením.

Při monitorování stavu vodních děl či jezer z hlediska erozně-akumulačních procesů se v současné době stále více uplatňují moderní metody bezkontaktního průzkumu pomocí sonarů, např. [5–10], či zařízení AUV (z angl. Autonomous Underwater Vehicle) [11]. Aby bylo možné posoudit míru změn morfologie VD, je třeba vycházet z porovnání s informacemi o návrhovém stavu uvedeném v patřičné projektové dokumentaci VD. Ty zahrnují technické zprávy a situační plány v tištěné podobě, mnohdy v různých formátech a v různé kvalitě „čitelnosti“. Tyto podklady je nezbytné převést do digitální podoby a provést jejich lokalizaci (georeferencování) do souřadnicového systému. Problematika digitalizace georeferencování historických map je široce popsána v odborné literatuře, např. [7, 12–14]. Praxe ukazuje, že jednotlivé podklady při zpracovávání vyžadují často individuální přístup vzhledem k době, kdy vznikly a s jakou přesností a pečlivostí byly vyhotoveny.

Předložený článek popisuje na příkladech pilotních lokalit postupy zvolené pro digitalizaci projektových dokumentací vybraných vodních děl různého stáří a velikosti. Výsledky těchto procesů jsou pak využity při prvotním srovnání získaného digitálního návrhového stavu VD se stavem současným, který reprezentují výsledky sonarového měření.

Metodika

Návrhový stav

Návrhový stav vodních děl shrnuje technická dokumentace. Většina podkladů je staršího data a je k dispozici pouze v tištěné podobě v různých formátech. Některé z nich jsou často poškozené, popř. hůře čitelné. Digitalizace situačních výkresů proto probíhala po částech s dostatečnými překryvy, které umožnily navázat jednotlivé díly na sebe.

Obr. 2. Příčné profi ly defi nované koncovými body a průsečíky s osou toku podle projektové dokumentace
Fig. 2. Position of cross profi les defi ned by their intersection with the channel‘s axis according to the situational plan

Přesná lokalizace jednotlivých částí výkresové dokumentace, tzv. georeferencování, byla prováděna v prostředí ESRI ArcGIS. Zpracováním po jednotlivých dílech bylo eliminováno riziko přenesení možných chyb georeferencování, které vyplývá z rozlišení skenovaných snímků, z jedné části výkresu do dalších. Jako vlícovací body byly využity objekty či jejich části zakreslené v situačních výkresech, které bylo možné identifikovat na současných mapách. Obvykle se jedná o rohy budov a význačné tvary projektovaných mostů, hrází či jezů. Pokud to bylo možné, bylo využito rovněž hranic parcel. V některých částech výkresů se však potřebné záchytné body nevyskytovaly. V takových případech byla lokalizace části výkresu provedena pomocí překryvu se sousedními výkresy a s využitím výrazných terénních tvarů, jako jsou charakteristické ohyby vrstevnic ve stržích apod. Následně byly z takto upravených podkladů do vektorových formátů převáděny základní informace o morfologii dna pilotních lokalit podle návrhového stavu.

Vodní nádrže

V případě vodních nádrží byly z georeferencovaných situačních výkresů do mapové vrstvy vykresleny linie vrstevnic, jimž byla přiřazena odpovídající nadmořská výška v souladu s použitým výškovým systémem (obr. 1).

Obr. 3. Schéma stanovení nadmořské výšky a vzdálenosti lomových bodů příčného profilu od průsečíku s osou
Fig. 3. Principle of the fracture points‘ positioning in a cross profi le based on the distance from the intersection of the channel axis with the cross profi le; calculation of the fracture points‘ altitude

Úseky vodních toků

Úprava koryt vodních toků je v situačních výkresech obvykle znázorněna osou koryta a polohou příčných profilů (obr. 2). Tyto profily jsou zaznamenány samostatně jako linie s lomovými body, z nichž některé mají uvedenu hodnotu nadmořské výšky Z podle použitého výškového systému. Každý výkres má znázorněnu referenční hladinu v nadmořské výšce Z0 (obr. 3). Podle této referenční hladiny byly digitalizované výkresy vyrovnány do vodorovné polohy, jednotlivé lomové body uloženy do bodové vrstvy a stanoveny jejich relativní kartézské souřadnice X, Y. Pomocí známých nadmořských výšek některých bodů, referenční hladiny a jejich vzdáleností vyjádřených vzájemnými rozdíly souřadnic XY bylo definováno měřítko zakreslení daného profilu. To pak umožňuje dopočítat nadmořské výšky Z bodů, u nichž tento údaj nebyl uveden, a také vzdálenosti d jednotlivých bodů od osy příčného profilu.

Následná část výpočtu spočívá v určení polohy (kartografických souřadnic) jednotlivých lomových bodů příčných profilů v mapě s využitím mapové vrstvy příčných profilů vykreslených ze situačního výkresu a také mapové vrstvy průsečíků příčných profilů s osou toku (blíže viz obr. 2). Dvojice koncových bodů každého příčného profilu určují směrový vektor u pomyslné přímky, na níž daný příčný profil leží. Snadným přepočtem je možné odvodit normálový vektor n této přímky, jejíž obecná rovnice má tvar:

kde n1, n2 je souřadnice normálového vektoru n,
X, Y  kartografické souřadnice libovolného bodu ležícího na přímce,
c konstanta dané přímky.

Při znalosti zeměpisných souřadnic X0, Y0 průsečíku příčného profilu s jeho osou a na základě vzdálenosti d jednotlivých lomových bodů příčného profilu od průsečíku je možné sestavit soustavu dvou rovnic o dvou neznámých Xi, Y(tj. hledané kartografické souřadnice konkrétního lomového bodu) ve tvaru:

Úpravou této soustavy je získána kvadratická rovnice, jejímž řešením jsou dva body ležící na příčném profilu ve stejné vzdálenosti od průsečíku s osou. Z této dvojice bodů je vybrán ten, který skutečně leží na příčném profilu. Využito bylo atributu vyjadřujícího příslušnost ke straně (P nebo L) příčného profilu a dvou polygonů pokrývajících dostatečně široký prostor v okolí koryta. Oba polygony se dotýkají svými hranicemi a kopírují přitom osu koryta. Hledané body leží v polygonu s totožným označením strany (obr. 4).

Obr. 4. Rozložení vypočítaných bodů vzhledem k jejich příslušnosti k pravé či levé straně vůči ose toku (O – průsečík s osou toku)
Fig. 4. Spatial distribution of the calculated fracture points, considering their affiliation to the right/left side from a channel‘s axis

Digitalizované příčné profily byly následně použity k vytvoření digitálního povrchu dna, resp. koryta, podle návrhového stavu. Body byly importovány do prostředí HEC-RAS, kde byla z příčných profilů vytvořena spojitá geometrie koryta toku, a to tak, že byly propojeny charakteristické lomové body, a na takto vzniklé spojnice byly umístěny meziprofily s rozestupem 5 m (obr. 5). V ojedinělých případech bylo nutné automatické propojení poupravit, aby tvar koryta netvořil překážku ve směru proudění.

Obr. 5. Vytváření 3D koryta v prostředí HEC-RAS
Fig. 5. Creation of the 3D model of a channel using HEC-RAS

Z prostředí HEC-RAS byla vytvořená geometrie přenesena do prostředí ESRI ArcGIS, kde byla s využitím doplňku HEC-Geo RAS dále zpracována do podoby nepravidelné trojúhelníkové sítě TIN (triangulated irregular network) a následně převedena na rastr s velikostí buňky 0,5 m. Stanovení tohoto rozlišení vychází ze sítě bodů zaměřených při echolokaci dna.

Současný stav

Měření batygrafie vodních děl bylo provedeno plavidlem CAPEREA (obr. 6) vycházejícím ze sériově vyráběného typu Sunchser 7518, poháněného dvojicí spřažených elektromotorů. Plavidlo bylo vybaveno skupinou zařízení pro určování zeměpisné polohy, senzorem zajišťujícím korekci náklonu podél hlavní a vedlejší osy a zařízením pro měření rychlosti zvuku a jeho profilu ve vodním sloupci. Samotné měření batygrafie bylo prováděno s pomocí dvou sonarů.

Obr. 6. Multifunkční plavidlo CAPEREA®
Fig. 6. Multipurpose watercraft CAPEREA® used for the sonar measurements of bottom

Sonarová jednotka GeoSwatch Plus Compact 500 kHz je systém určený pro měření batygrafie mělkých vod do hloubky 50 m. Umožňuje vytvářet digitální model dna se statistickou přesností 1,5 mm. Sonar SES-2000 Compact je schopen pronikat do mělkých geologických vrstev a umožňuje zaznamenávání podélného profilu sedimentovaných vrstev a měření jejich mocnosti. Je založen na vysílání dvou různých frekvencí, které spolu interagují a pomocí jejich sumy či rozdílu je možné získat informaci o stratigrafii dna.

Výstupem měření batygrafie je mračno bodů obsahujících údaje o úrovni dna v dané pozici (tj. kartografické souřadnice a nadmořskou výšku). Z této bodové vrstvy byla vytvořena nepravidelná trojúhelníková síť (TIN), která byla následně převedena na rastr s rozlišením 0,5 m. Uvedený postup je znázorněn na obr. 7.

Obr. 7. Postup převodu bodové vrstvy z echolokačního měření části úseku vodního toku (a) na nepravidelnou trojúhelníkovou síť (b) a rastr s rozlišením 0,5 m (c)
Fig. 7. Transformation of the point layer from the sonar measurements of a channel‘s bottom (a) into triangulated irregular network (TIN) (b) and into raster with grid size of 0.5 m (c)

Pilotní lokality

K otestování postupu digitalizace byl primárně vybrán úsek řeky Moravy ve vzdutí nad jezem v Kroměříži po jez v Bezměrově, resp. úsek mezi ř. km 179,230 a 182,515 (obr. 8).

Obr. 8. Zájmový úsek řeky Moravy ve vzdutí jezu mezi Kroměříží a Bezměrovem s vyznačením polohy jednotlivých částí situačního výkresu
Fig. 8. Study river reach of the Morava River between Kroměříž and Bezměrov with the indication of individual parts of the situational plan

Na vodní nádrži Letovice bylo za účelem verifikace sonarového měření provedeno geodetické měření části dna. Měření probíhalo nedaleko obce Vranová při pravém břehu zátopy v pásu přibližně 200 m dlouhém a 30 m širokém (obr. 9). Měrná lať byla z plavidla spouštěna na úroveň dna, jehož nadmořská výška byla poté zaměřena. Podobné kontrolní zaměření bylo provedeno také v pracích Kubinského a kol. [7, 8], ovšem s tím rozdílem, že měrná tyč byla opatřena plochou patkou.

Výsledky

Návrhový a současný stav dna

Popis změny morfologie dna vodních děl za účelem následného vyhodnocení dopadu na jejich funkce a posouzení vhodnosti a efektivity potenciálních zásahů vyžaduje znalost návrhového stavu, k němuž je vztahována současná úroveň dna. Vzhledem ke stáří a povaze situačních výkresů může být jejich digitalizace a správná lokalizace obtížná. Postup, který lze považovat za obecně správný a přesný, je popsán v metodice. Při praktickém zpracování konkrétních podkladů však místy vyvstaly určité problémy.

Obr. 9. Kontrolní geodetické měření úrovně dna části vodní nádrže Letovice
Fig. 9. Verifi cation of the sonar measurement of a bottom in a near-shore zone at the Letovice Dam – geodetic measurement

U historických situačních výkresů je často obtížně zjistitelné, v jakém souřadném systému byly vyhotoveny, což potvrzují i další autoři [12, 15]. Především u delších koryt toků v důsledku toho může docházet k deformaci georeferencovaného podkladu a k odchylování koryt, byť mírnému, od jejich skutečné polohy. Zejména na okrajích koryta proto mohou být při rozdílové analýze zaznamenávány zdánlivé prohlubně či elevace (obr. 10). K tomuto přistupuje i skutečnost, že původně projektované úpravy nemusely být realizovány v celém rozsahu (tzn. po celé délce a na obou stranách koryta). V případech, kdy průběh původního koryta byl velmi podobný průběhu koryta navrhovaného avšak nerealizovaného, může být obtížné rozhodnout, kterou z variant považovat za skutečný návrhový stav.

Obr. 10. Zdánlivé elevace a prohlubně na okrajích měřeného koryta způsobené vychýlením koryta při georeferencování situačního výkresu
Fig. 10. Apparent elevations and depressions at the edges of the river channel caused by the channel deviation when georeferencing the situational plan

Přesnou lokalizaci historických situačních plánů lze provést s pomocí důkladné rekognoskace konkrétního úseku. To představuje geodetické zaměření břehových linií, sklonů břehů, v ideálním případě i zaměření topografie dna. Podle nich je následně možné korektně umístit příčný profil a jedině tak lze určit, zda byla úprava koryta podle projektu provedena, případně v jakém rozsahu. Ani tento postup však nemusí vést k jednoznačnému výsledku. Břehy totiž mohu být za dlouhé období silně ovlivněny prorůstajícími kořeny stromů, erozí nebo dodatečnými změnami v opevnění, což se může projevit ve změně tvaru průtočných profilů. Přesto zůstává posouzení kombinace morfologických charakteristik koryta a následná manuální korekce umístnění konkrétního příčného profilu jediným způsobem, jak provést jeho lokalizaci co nejpřesněji.

Přesnost echolokátoru

Využití sonaru jako spolehlivého zdroje dat bylo prezentováno v řadě prací věnujících se monitoringu dna a erozně-akumulačních procesů probíhajících v přírodních jezerech (např. [6, 10, 12]), vodních nádržích (např. [7, 8]) i vodních tocích (např. [11, 16]). V rámci této práce byla přesnost výstupů měření echolokátoru testována v pilotní lokalitě na vodní nádrži Letovice. Odchylky zaměřených úrovní dna jsou uvedeny na obr. 11, z něhož je patrné, že nepřevažují kladné ani záporné hodnoty, z nichž by bylo možné usuzovat na systematické pod- či nadhodnocovaní úrovně dna echolokačním měřením. Tuto nejednoznačnost je možné odůvodnit velkou členitostí dna, v důsledku které může nevýznamný posun měrné latě z důvodu obtížné dlouhodobé stabilizace plavidla vést k naměření odchýlených hodnot [8]. Obecným problémem při hodnocení úrovně dna touto metodou je rovněž zrnitostní struktura zvodnělých sedimentů, projevující se jejich rozličnou tuhostí vůči spouštěné měrné lati. V případě měkkých sedimentů může docházet k jejímu zaboření do dna a zaměření větší hloubky oproti skutečné. Tento problém by bylo možné jistě alespoň částečně řešit opatřením měrné tyče zmíněnou plochou patkou.

Obecné hodnocení morfologických změn

Porovnáním současné úrovně dna s návrhovým stavem umožnilo identifikovat erozně-akumulační procesy probíhající ve vzdutí jezu v Kroměříži na řece Moravě. Z provedených měření jsou patrná především výrazná zahloubení dna podél konkávních břehů, kde maximální hloubky dosahují hodnot přes 7 m. Nápadné prohlubně se vyskytují také při vyústěních přítoků do hlavní trati a za mostními pilíři. V uvedených místech dochází ke vzniku turbulentního proudění, v němž je sedimentace unášených částic znemožněna. Obecně je však zřetelné, že převažujícím jevem je akumulace plavenin, což je způsobeno zpomalením proudění, a tedy sníženou transportní kapacitou toku v tomto úseku.

Obr. 11. Výsledek porovnání nadmořských výšek v kontrolních bodech vůči úrovni dna zaměřeného echolokátorem v pilotní lokalitě VN Letovice poblíž obce Vranová
Fig. 11. Differential comparison of the altitudes from the verification measurement and the altitudes measured using the sonar in the neer-shore zone of the Letovice Dam

Uvedeným postupem je možné stanovit i konkrétní kubatury deponovaného materiálu. Vzhledem k výše uvedeným problematickým aspektům je však nutno nejprve všechny nejednoznačnosti odstranit. Posouzení potřebnosti případných zásahů však musí vycházet z podrobného vyhodnocení jednotlivých částí koryta, v nichž mohou probíhat odlišné, specifické, erozně-akumulační procesy. Z hlediska logistického je také nutno uvažovat, že skutečná kubatura sedimentů může být nižší z důvodu jejich zvodnění.

Poděkování

Článek vznikl v rámci zpracování projektu KUS QJ1520267 Systém řízení a údržby VH infrastruktury financovaného z prostředků Ministerstva zemědělství ČR.