Masterproef Optimalisatie van de verwerkingsprocedure voor het ruimen van het Groot Schijn

Masterproef Optimalisatie van de verwerkingsprocedure voor het ruimen van het Groot Schijn Studiegebied Industriële wetenschappen en technologie Opleiding Master in de industriële wetenschappen: milieukunde Promotoren Kathelijne Velghe Ward De Cooman, VMM Academiejaar 2010-2011

Birgit Van Boxstael

Howest – departement Academische bachelor- en masteropleidingen, Graaf Karel de Goedelaan 5, 8500 Kortrijk

Woord vooraf Deze masterproef werd geschreven met het oog op het behalen van de graad van Master in de Industriële Wetenschappen: Milieukunde. Als ik terugblik op dit jaar kan ik concluderen dat het geheel aan acties (veldwerk, onderzoek, …) om tot dit werk te komen, een absoluut lichtpunt waren in dit hectische masterjaar. Ik ben dan ook enige dank verschuldigd aan het Team Waterbodemmeetnet van de Vlaamse Milieumaatschappij voor het opfleuren van elke masterproefdag die in het teken van het praktisch onderzoek stond. Daarbij wil ik in het bijzonder de heer W. De Cooman, mijn externe promotor, bedanken voor de geweldige begeleiding. Maar ook de heren L. De Temmerman, N. Meyskens, K. De Wit, K. Dehaemers en T. Beckx zou ik graag bedanken voor hun geduld, hulp en steun tijdens het voorbije jaar. Vervolgens wil ik mevrouw K. Velghe bedanken voor de begeleiding en het veelvuldig nalezen van deze masterproef. Tenslotte wil ik nog mijn ouders bedanken voor het mogelijk maken van deze studie.

Birgit Van Boxstael

I

Abstract Flanders gets quantity problems in the watercourses caused by human needs such as shipping, housing or agriculture that require greater flow capacity and water depth. Also the various floods are a quantity problem. On the other hand there’s a need to clear the watercourses. All these problems can be solved by dredging. There is a need to optimize the processing procedure for measuring the quantity of sediment in the Flemish waterways. In the past no studies were available that did research on how to survey bends in watercourses in an efficient and accurate way. There is also need for a good method for averaging (smoothing) digital terrain models (DTM). This study was made to understand the processing with Terramodel software and understand the mathematical background. With this knowledge it was possible to optimize smoothing and surveying. The most efficient and accurate method for surveying bends and smoothing a DTM was made. Three types of curves were characterized and classified for practical surveying. Investigated curves can be classified in 3 classes. Class 1 contains the curves from 100°135° and needed to be measured in 8 profiles. Curves of 80° needed 5 profiles and curves of 30°-50° needed 4 profiles. The best result for smoothing the DTM was obtained by using a grid size of 0.9-0.9 m. These results gives the VMM a good support for surveying on an accurate and efficient way.

II

Inhoudsopgave WOORD VOORAF ............................................................................................................ I ABSTRACT ...................................................................................................................... II LIJST VAN AFKORTINGEN ............................................................................................ V LIJST VAN TABELLEN .................................................................................................. VI LIJST VAN FIGUREN .................................................................................................... VII 1

INLEIDING ................................................................................................................ 1

2

BESCHRIJVING VMM............................................................................................... 3

3

2.1

Algemeen ........................................................................................................... 3

2.2

Geschiedenis ...................................................................................................... 4

2.3

Activiteiten .......................................................................................................... 4

2.4

Vestigingen ......................................................................................................... 4

2.5

Organogram ....................................................................................................... 5

LITERATUURSTUDIE ............................................................................................... 6 3.1

Waterbodem ....................................................................................................... 6

3.1.1

Transport van contaminanten in de waterbodem ......................................... 7

3.1.2

Kwaliteit van de waterbodem ....................................................................... 8

3.1.3

Baggeren en saneren .................................................................................. 9

3.1.4

Ruimingsmeetnet ........................................................................................11

3.2

Kwantiteitsbepaling ............................................................................................13

3.2.1

Handmatige technieken ..............................................................................13

3.2.2

Akoestische technieken ..............................................................................16

3.2.3

Elektromagnetische techniek ......................................................................19

3.2.4

Positiebepaling ...........................................................................................20

3.2.5

Keuze meetmethoden.................................................................................27

3.2.6

Spreiding meetmethoden ............................................................................27

3.3

Dataverwerking ..................................................................................................28

III

3.3.1

Algemeen ...................................................................................................28

3.3.2

Terramodel .................................................................................................29

3.4 4

MATERIAAL EN METHODEN .................................................................................38 4.1

Inventarisatie materiaal...............................................................................38

4.1.2

Methode .....................................................................................................39

Smoothing .........................................................................................................40

4.2.1

Inventarisatie materiaal...............................................................................40

4.2.2

Methode .....................................................................................................41

RESULTATEN EN CONCLUSIES............................................................................44 5.1

Onderzoek Vlaamse bochten .............................................................................44

5.1.1

Bocht van 40°.............................................................................................44

5.1.2

Bocht van 80°.............................................................................................46

5.1.3

Bocht van 105°...........................................................................................47

5.1.4

Bocht van 120°...........................................................................................49

5.1.5

Bocht van 135°...........................................................................................51

5.1.6

Besluit: procedure inmetingen ....................................................................53

5.2

6

Onderzoek van Vlaamse bochten ......................................................................38

4.1.1

4.2

5

Groot Schijn .......................................................................................................37

Smoothing .........................................................................................................54

5.2.1

Volume vergelijken .....................................................................................54

5.2.2

Visuele beoordeling ....................................................................................54

5.2.3

Besluit smoothing .......................................................................................55

BESLUIT ..................................................................................................................56

BIBLIOGRAFIE BIJLAGEN

IV

Lijst van afkortingen AOW

Afdeling Operationeel Waterbeheer

ARW

Afdeling Rapportering Water

DGPS

Differentieel Global Positioning System

DEM

Digital Elevation Model

DTM

Digital Terrain Model

GPS

Global Positioning System

IVA

Intern Verzelfstandigd Agentschap

LIDAR

Laser Imaging Detection And Ranging

MIRA

Milieurapport

NM

Navigation message

RTK

Real Time Kinematic

RWZI

Rioolwaterzuiveringsinstallatie

TIN

Triangular Irregular Network

VLAREA

Vlaams Reglement voor Afvalvoorkoming en -beheer

VMM

Vlaamse Milieumaatschappij

V

Lijst van tabellen Tabel 3-1 Overzicht plaatsbepalingsmethodes i.f.v. meettechnieken ...............................26 Tabel 3-2 Overzicht gebruik van meetmethoden ............................................................27 Tabel 5-1 Resultaten bocht van 40°.................................................................................44 Tabel 5-2 Resultaten bocht van 80°.................................................................................46 Tabel 5-3 Resultaten bocht van 105°...............................................................................47 Tabel 5-4 Resultaten bocht van 120°...............................................................................49 Tabel 5-5 Resultaten bocht van 135°...............................................................................51 Tabel 5-6 Resultaten smoothing ......................................................................................54

VI

Lijst van figuren Figuur 2-1 Beleidsdomeinen departement leefmilieu, natuur en energie .......................... 3 Figuur 2-2 Organogram VMM .......................................................................................... 5 Figuur 2-3 Organogram Team Meetnet Waterbodem....................................................... 5 Figuur 3-1 Waterbodem .................................................................................................. 6 Figuur 3-2 Factoren die het transport beïnvloeden .......................................................... 7 Figuur 3-3 Resultaten Triade 2006-2009 ......................................................................... 9 Figuur 3-4 Baggerwerken ...............................................................................................10 Figuur 3-5 Saneringscyclus van een waterbodem...........................................................12 Figuur 3-6 Uitschuifbare stok ..........................................................................................13 Figuur 3-7 Niet inschuifbare peilstok ...............................................................................14 Figuur 3-8 Peilhengel ......................................................................................................15 Figuur 3-9 Singlebeam en Multibeam ............................................................................17 Figuur 3-10 Werking Side Scan sonar ............................................................................18 Figuur 3-11 Subbottom profiler .......................................................................................18 Figuur 3-12 Georadar ....................................................................................................19 Figuur 3-13 Waterdruktechniek .......................................................................................20 Figuur 3-14 Uitzetwerk ....................................................................................................21 Figuur 3-15 Totaalstation ................................................................................................24 Figuur 3-16 Laser altimetrie ............................................................................................25 Figuur 3-17 Van links naar rechts: contour, grid, TIN .....................................................29 Figuur 3-18 Verwerking via Terramodel ..........................................................................30 Figuur 3-19 Omgeschreven cirkel ...................................................................................31 Figuur 3-20 Edge flipping ................................................................................................31 Figuur 3-21 Divide and Conquer .....................................................................................32 Figuur 3-22 Toevoeging van een punt in de driehoek. ....................................................33 Figuur 3-23 Interpolatie ...................................................................................................34 Figuur 3-24 TIN Interpolatie ...........................................................................................35

VII

Figuur 3-25 Situering Groot Schijn ..................................................................................37 Figuur 4-1 Inmeting .........................................................................................................39 Figuur 4-2 Onderzochte traject .......................................................................................41 Figuur 4-3 Inmeting .........................................................................................................42 Figuur 4-4 Raster grootte ................................................................................................42 Figuur 5-1 Visuele voorstelling foutenpercentage van een bocht van 40°.......................44 Figuur 5-2 Verschilkaart bocht van 40°, a (5 profielen), b (4 profielen), c (3 profielen) ....45 Figuur 5-3 Visuele voorstelling foutenpercentage van een bocht van 80°.......................46 Figuur 5-4 Verschilkaart bocht van 80°, a (8 profielen), b (7 profielen), c (6 profielen), d (5 profielen), e (4 profielen) ..................................................................................................47 Figuur 5-5 Visuele voorstelling foutenpercentage van een bocht van 105°.....................48 Figuur 5-6 Verschilkaart bocht van 105°, a (10 profielen), b (9 profielen), c (8 profielen), d (7 profielen), e (6 profielen) ..............................................................................................48 Figuur 5-7 Visuele voorstelling foutenpercentage van een bocht van 120°.....................49 Figuur 5-8 Verschilkaart bocht van 120°, a (9 profielen), b (8 profielen), c (7 profielen), d (6 profielen), e (5 profielen), f (4 profielen) .......................................................................50 Figuur 5-9 Visuele voorstelling foutenpercentage van een bocht van 135°.....................51 Figuur 5-10 Verschilkaart bocht van 135°, a (12 profielen), b (11 profielen), c (10 profielen), d (9 profielen), e (8 profielen), f (7 profielen), g (6 profielen), h (5 profielen) ....52 Figuur 5-11 Indeling in klassen .......................................................................................53 Figuur 5-12 Visuele voorstelling foutenpercentage smoothing ........................................54 Figuur 5-13 Visuele voorstelling smoothing 0.9 x 0.9 (m) ................................................55

VIII

1 Inleiding Een natuurlijke waterloop ontwikkelt een toestand waarbij waterafvoer, erosie, sedimentatie en sedimenttransport in evenwicht zijn. Menselijke ingrepen brengen dit evenwicht uit balans met een toename van de toevoer van slib en sedimenten in de waterlopen tot gevolg. De toevoer van slib is afkomstig van de lozing van ongezuiverd afvalwater terwijl de toevoer van sediment afkomstig is van landerosie en oevererosie. Deze verhoogde sedimenttoevoer veroorzaakt niet alleen een kwantiteitsprobleem maar ook een kwaliteitsprobleem. Het kwaliteitsprobleem komt tot stand door de continue wisselwerking tussen water en bodem. Een vervuilde bodem zal door nalevering het water blijven vervuilen maar ook door diffuse lozingen van het omliggende land. Denk vb. aan afspoeling van land en oever, dat verontreinigd is met pesticiden. Ook puntlozingen van zwevende stof (slib uit RWZI of overstorten) beladen met contaminanten zorgt nog steeds voor een impact op de waterloop/waterbodem. Hoewel de meeste verontreiniging van waterbodems wordt toegeschreven aan historische verontreiniging, worden toch nog steeds concentraties van contaminanten teruggevonden die niet te verklaren zijn. Sommige contaminanten zijn reeds lang verboden maar worden nog steeds teruggevonden. Daarentegen krijgt men te maken met een kwantiteitsprobleem door menselijke behoeften zoals scheepvaart, huisvesting of landbouw die een grotere afvoercapaciteit en waterdiepte vereisen. Ook de diverse overstromingen vallen onder het kwantiteitsprobleem. Deze problemen kunnen gedeeltelijk opgelost worden door voorzorgsmaatregelen die een aangroeiende sedimentatie voorkomen. Volgens verschillende studies is gebleken dat erosiebestrijdende middelen zeer efficiënt zijn. Maar dat vergt een grote aanpassing van verschillende gebruiken van land en water (zie sensibiliseringsprogramma van landbouwers). Enkele voorbeelden hiervan zijn: de beperking van overstorting van de rioolstelsels, goed ontwikkelde oeverstroken die oevererosie voorkomen, … . Soms zijn ruimingswerken van een waterloop, gekoppeld aan afvalwaterzuivering en erosiebestrijding, noodzakelijk. Vooraleer ze uitgevoerd kunnen worden dient de kwaliteit en kwantiteit in kaart gebracht worden. M.b.v. deze gegevens kan het kostenplaatje van de werken bepaald worden en is verdere evaluatie van de waterloop mogelijk. De kwantiteitsbepaling werd in het verleden uitbesteed aan externe bureaus. Na een kosten-batenanalyse kwam de Vlaamse Milieumaatschappij tot de conclusie dat een kwantiteitsbepaling in eigen beheer kostenefficiënter is. De grote uitdaging van deze masterproef was het optimaliseren van de verwerkingsprocedure die uiteindelijk leidt tot een kwantiteitsbepaling. De vijf doelstellingen van deze masterproef kunnen als volgt omschreven worden: 1. Inzicht krijgen in de verwerkingswijze in Terramodel ®. 2. De theoretische achtergrond van Terramodel koppelen aan de praktische inmeting met een veldstrategie tot gevolg. 3. Onderzoek naar de meest accurate en efficiëntste methode voor het inmeten van enkele typische ‘Vlaamse’ bochten in het Groot Schijn.

1

4. Opstellen van een gestandaardiseerde procedure voor het praktisch inmeten van 3 types bochten. 5. Onderzoek naar de mogelijkheid om de modellen uit te middelen (smoothen).

Mijn belangrijkste rol in de VMM was inzicht te verwerven in de verwerkingsprocedure van de kwantiteitsbepaling en deze (software matige) te optimaliseren. In de eerste weken van mijn masterproef periode maakte ik kennis met Terramodel. Door op zoek te gaan naar de wiskundige beginselen van dit programma, kon ik mijn doelstellingen omtrent de optimalisatie volmaken.

2

2 Beschrijving VMM 2.1

Algemeen

Figuur 2-1 Beleidsdomeinen departement leefmilieu, natuur en energie (Departement Leefmilieu, natuur en energie van de Vlaamse overheid, 2010)

De Vlaamse overheid werd ingedeeld in 13 beleidsdomeinen. Elk beleidsdomein bestaat uit een departement en Verzelfstandigde Agentschappen. (Vlaanderen.be, 2010) Enerzijds staan de departementen in voor de beleidsplanning. Anderzijds houden de verschillende agentschappen zich bezig met de beleidsuitvoering (Van Boxstael, 2009). 2009) Door de hervorming werden 13 ministeries gevormd. Elk ministerie bestaat uit het departement en de IVA (Intern Verzelfstandigde Agentschappen) zonder rechtspersoonlijkheid. De Intern Verzelfstandigde Agentschappen staan o.l.v. een bevoegde minister. Wanneer een IVA rechtspersoonlijkheid heeft mag het in eigen naam contracten sluiten zonder tussenkomst van de minister. De Vlaamse Milieumaatschappij is een IVA met rechtspersoonlijkheid. Concreet wil dit zeggen dat de VMM werd opgericht bij decreet, door het Vlaams Parlement. Mevr. J. Schauvliege is momenteel el Vlaams minister bevoegd voor Openbare werken, werken Energie, Leefmilieu en Natuur (Van Boxstael, 2009). 2009)

3

2.2

Geschiedenis

De VMM werd opgericht bij decreet van het bestuurlijk beleid van 12 december 1990 als Vlaamse openbare instelling. Door het decreet van 7 mei 2004, dat in werking trad op 1 april 2006, werd de VMM omgevormd tot een IVA met rechtspersoonlijkheid (Van Boxstael, 2009). De taken werden in het decreet van 7 mei 2004 vastgelegd (VMM, 2010).

2.3

Activiteiten

De taken van de VMM zijn, zoals eerder vermeld, decretaal vastgelegd en onderverdeeld in drie thema’s: •

lucht

•

milieu

•

water

De VMM heeft de volgende taken (VMM, 2010): •

meten en controleren van de kwaliteit en kwantiteit van water

•

watersystemen beheren

•

heffing op watervervuiling innen

•

heffing op grondwaterwinning innen

•

adviseren over milieuvergunningen

•

planning en toezicht op de zuiveringsinfrastructuur

•

luchtkwaliteit bewaken

•

doet beleidsvoorstellen

•

inventarisatie van de lozingen

•

opstellen MIRA

2.4

Vestigingen

De hoofdzetel van de VMM bevindt zich in Erembodegem. Verder bezit de VMM vestigingen in de 5 provincies (Van Boxstael, 2009).

4

2.5

Organogram

Figuur 2-2 Organogram VMM (VMM, 2010)

Figuur 2-3 Organogram Team Meetnet Waterbodem (VMM, 2010)

5

3 Literatuurstudie 3.1

Waterbodem

Het begrip waterbodem werd gedefinieerd, in het decreet van 18 juli 2003 betreffende het integraal waterbeleid, als de bodem van een oppervlaktewaterlichaam die altijd of een groot gedeelte van het jaar onder water staat , zie Figuur 3-1 (Belgisch Staatsblad, 2010). 2010)

Figuur 3-1 Waterbodem (VMM, 2003)

Onder waterbodem verstaan we de historische bodem en de sedimentlaag. Deze sedimentlaag wordt gevormd door afzetting van zwevende stoffen. Zwevende stoffen kunnen van industriële of huishoudelijke oorsprong zijn, maar ook landbouwactiviteit of afspoeling kunnen aan het ontstaan van zwevend stof liggen. De waterbodem waterbodem doorstaat een continue cyclus van erosie en sedimentatie. Het is vanzelfsprekend dat de Vlaamse waterbodems, gezien de oorsprong van de zwevende stoffen, in het merendeel van de gevallen verontreinigd zijn.. De afzetting van zwevend stof zal het grootst grootst zijn in gebieden met een lage stroomsnelheid. Zo zal verontreinigd zwevend stof eerder terug te vinden zijn in de gebieden den met een lage stroomsnelheid (Vignati, Pardos, & Diserens, 2003). 2003) Dit verklaart waarom er in de Vlaamse laamse waterlopen een dikkere (vervuilde) sedimentlaag te vinden is dan in de Waalse waterlopen. De waterbodem is van cruciaal belang in het waterecosysteem. Het is de standplaats, rustplaats en voedselbron voor diverse fauna en flora.

6

3.1.1

Transport van contaminanten conta in de waterbodem

Er bestaat een continue wisselwerking tussen water en bodem. Als gevolg van een steeds zoeken naar een chemisch evenwicht tussen water en bodem zal een contaminatie met microverontreinigingen van het water een contaminatie van de waterbodem veroorzaken. Een vervuilde bodem zal door nalevering het al dan niet propere water vervuilen. Contaminanten kunnen in verschillende vormen worden getransporteerd: • • •

Opgelost in het bodemwater; onder gasvormige toestand in de gasfase; geadsorbeerd op bodemcolloïden die zelf door het bodemwater getransporteerd worden.

Verder gaat men van de ene transportvector over naar de andere via via evenwichtsreacties. Zo bepaalt de constante van Henry het evenwicht tussen de gasfase en waterfase en de verschillende distributiecoëfficiënten het evenwicht tussen de waterige fase en bodemcolloïden, zie Figuur 3-2 (Kluwer, 1998).

Figuur 3-2 Factoren die het transport beïnvloeden (Kluwer, 1998)

In hoeverre een contaminant in staat is te migreren in de bodem,, en op welke manier (via waterfase, vastgehecht aan bodemmatrix of via gasfase), hangt af van verschillende polluenteigenschappen. igenschappen. De belangrijkste worden kort besproken.

3.1.1.1

Wateroplosbaarheid

De wateroplosbaarheid geeft aan in welke mate een product zich in de waterige fase zal verdelen. Hoe groter de wateroplosbaarheid, hoe sneller de contaminanten door een bodem getransporteerd orteerd worden via de verzadigde zone1.

1

Poriën die volledig gevuld zijn met water.

7

3.1.1.2

Dampspanning

De dampspanning is een maat voor de druk die door een product wordt uitgeoefend op de wanden van een gesloten ruimte. Hoe hoger de dampspanning, hoe groter de vervluchtigbaarheid. Meestal zijn wateroplosbare contaminanten niet geneigd te vervluchtigen.

3.1.1.3

Constante van Henry

Deze eigenschap geeft weer in welke mate een product vervluchtigt van de waterige fase naar de gasfase. Hieruit kan men bepalen in welke fase een bepaald product zich concentreert (Kluwer, 1998).

3.1.1.4

Adsorptie/desorptie

Anorganische en organische contaminanten hebben de neiging te adsorberen aan de bodemmatrix. Deze adsorptie gebeurt voornamelijk aan de kleifractie en organisch materiaal. Klei is verantwoordelijk voor de adsorptie van polaire en hydrofiele producten. Organisch materiaal is in staat om zowel hydrofiele als hydrofobe producten te adsorberen. Adsorptie en desorptie worden weergegeven via distributiecoëfficiënten. Een belangrijke coëfficiënt is de octanol/waterverdelingcoëfficiënt (Kow). Deze geeft de verdeling weer van een stof tussen n-octanol en water. De Kow- waarde geeft een idee over het hydrofobisch gedrag van een bepaalde contaminant. Hoe hoger deze waarde, hoe groter de neiging om geadsorbeerd te worden aan de bodemmatrix. Microverontreinigingen met een grote Kow- waarde, worden gemakkelijk geadsorbeerd, waardoor ze lang in de waterbodem aanwezig blijven en een diepe penetratie in de bodem mogelijk is (Kluwer, 1998). Tijdens baggerwerken moet rekening gehouden worden met de mogelijkheid dat sediment gebonden contaminanten terug worden vrijgegeven (Daniels, House, Rae, & Parker, 2000).

3.1.2

Kwaliteit van de waterbodem

De kwaliteit van de Vlaamse waterbodems is één van de verantwoordelijkheden van de VMM. In maart 2000 werd het 'routinematig meetnet waterbodemkwaliteit' opgestart. Dit meetnet heeft volgende doelen: •

De inventarisatie en monitoring van de kwaliteit van de waterbodem binnen het Vlaamse gewest aan een frequentie van één maal per 4 jaar, om prioritaire saneringsplaatsen voor de waterbodems te inventariseren.

•

Het bepalen van de waterbodemkwaliteit door middel van indexen en waterbodemkwaliteitsklassen aan de hand van de triade beoordeling.

•

Het toetsen van de meetgegevens aan de huidige wettelijke normen, zoals voorwaarden voor afvalvoorkoming (VLAREA).

•

Het nagaan van de impact van specifieke acties op de waterbodemkwaliteit.

•

Het kennen van de kwaliteit als basisinformatie voor het uitwerken van BBP's, advisering vergunningen, enz. 8

Het meetnet bevat 600 meetplaatsen. Per jaar worden 150 meetplaatsen onderzocht volgens de Triade. De Triade is een beoordelingsmethode waarbij elke meetplaats onderworpen wordt aan drie verschillende onderzoeken: onderzoeken: fysicochemisch, biologisch biolog en ecotoxicologisch. De e combinatie van deze onderzoeken stelt de VMM in staat om waterbodems te rangschikken in functie van toenemende prioriteit voor saneringsonderzoek in het kader van het ecologisch herstel van rivieren/beken. rivieren/bek Uiteraard zijn deze saneringen slechts zinvol wanneer men het lozen van verontreinigingen stop zet. 2006 2009 weergegeven. In Figuur 3-3 worden de resultaten van de Triade van periode 2006-2009 Dit taartdiagram geeft weer weer in hoever het saneren van waterbodems prioritair is in Vlaanderen.

Figuur 3-3 Resultaten Triade 2006-2009 2006 (VMM, 2010)

3.1.3

Baggeren en saneren

Baggeren wordt algemeen beschouwd als al het weghalen van overtollig sediment uit waterlopen/zeeën. Bij verontreiniging van dit sediment dient men over te gaan op het saneren/zuiveren van dit sediment. Baggeren is van cruciaal belang bij volgende situaties: •

Verbeteren van de kwaliteit van de waterloop. wa

•

Kans op overstromingen door een te hoge sedimentafzetting.

•

Scheepvaartproblemen door een te dikke sedimentlaag.

Indien men baggert om de kwaliteit van een waterloop te verbeteren kan men spreken over het saneren van de waterloop. Men zal de verontreinigde verontreinigde specie weghalen, waardoor de continue wisselwerking tussen de waterbodem en het water geen verontreinigd water meer kan teweeg brengen. Dit zal uiteindelijk leiden tot een verbetering van de kwaliteit van de waterloop.

9

Figuur 3-4 Baggerwerken

Het uitgegraven sediment wordt baggerspecie genoemd. De toepassingsmogelijkheden van baggerspecie hangen af van de milieuhygiënische en bouwtechnische kwaliteit van de specie. Of baggerspecie al dan niet mag gebruikt worden als bodem of bouwstof is via wettelijke normen vastgelegd in het VLAREA (Vlaams Reglement voor Afvalvoorkoming en beheer). Deze zijn van toepassing op het sediment vanaf het moment dat het wordt gebaggerd. Vanaf dit ogenblik wordt de specie als afvalstof beschouwd (zelfs al is i de specie niet-verontreinigd) (Vito, 2007). 2007) Van de jaarlijks te ruimen baggerspecie voldoet 82% aan de VLAREA-normen VLAREA voor hergebruik. In dit getal zit de specie van zowel vaarbare als onbevaarba onbevaarbare waterlopen. Hierbij dient opgemerkt gemerkt te worden dat 23% voldoet voor hergebruik als bodem en 59% voldoett voor hergebruik als bouwstof. Baggerspecie wordt in Vlaanderen voor volgende toepassingsmogelijkheden gebruikt (Vito, 2007): • • • •

verspreiding op de oever (oeverdeponie); terugklepping van baggerspecie in de waterloop of de zee; behandeling van de specie, gevolgd door hergebruik als bodem of bouwstof; storting van specie.

10

3.1.4

Ruimingsmeetnet

Aangezien de VMM verantwoordelijk is voor het onderhoud van de 1ste categorie waterlopen, werd een ruimingsmeetnet opgestart. Dit ruimingsmeetnet staat, binnen de VMM, onder de bevoegdheid van de Afdeling Operationeel Waterbeheer (AOW). Deze afdeling houdt zich ook bezig met de voorbereiding van het dossier waarin de kwaliteit en kwantiteit van de betrokken waterloop wordt verduidelijkt. Op basis van dit dossier kan een aannemer de omvang in kostprijs van de gevraagde werken inschatten en een offerte maken. Het team meetnet waterbodem van de Afdeling Rapportering Water (ARW) heeft de taak op zich genomen om de kwantiteitsbepaling van de sedimenten in eigen beheer te doen. In het verleden werden deze inmetingen uitbesteed aan externe bedrijven. De bepaling in eigen beheer brengt verschillende voordelen met zich mee: •

ARW kan kwaliteit en kwantiteit gemakkelijk aan elkaar koppelen omdat deze afdeling over een monitoringsmeetnet beschikt die de toestand en de trend van de Vlaamse waterlopen inventariseert en rapporteert;

•

alle gegevens (meetplaatsen, werkplannen, parameterprojecten, ...) worden gecentraliseerd in een waterbodemkwaliteitsdatabank;

•

in de Waterbodemkwaliteitsdatabank werd een kwantiteitsmodule toegevoegd.

In actie 8.2 van het sectoraal Uitvoeringsplan Bagger- en Ruimingspecie werd de opdracht gegeven voor het opstarten van een centralisatie van kwantitatieve gegevens zoals deze uitgevoerd werd in de kwantiteitsdatabank. De koppeling van kwantiteit en kwaliteit vereenvoudigt, verbetert, versnelt en standaardiseert het rapporteren. In de kwantiteitsmodule van de waterbodemkwaliteitsdatabank worden volgende gegevens weergegeven: •

hoeveelheid geruimde specie;

•

kwaliteitstoestand voor en na een sanering (is de toestand verbeterd na de ruiming?);

•

hoeveelheid geruimde specie op jaarbasis per provincie, per bekken … ;

•

verzameling van gegevens van onder andere aannemer, offerte … .

Het ruimingsmeetnet gaat als volgt te werk:

11

Figuur 3-5 Saneringscyclus van een waterbodem

12

3.2

Kwantiteitsbepaling

Aan de hand van de locatie en het doel van de kwantiteitsbepaling, wordt een bepaalde methode aangewend. In dit deel wordt een opsomming van de verschillende kwantiteitsbepaling methoden gegeven.

3.2.1

Handmatige technieken

3.2.1.1 Peilstok De peilstok of prikstok is de meest gebruikte methode voor het handmatig peilen. Er worden verschillende uitvoeringsvormen van de peilstok gebruikt. De meest gebruikte peilstok is een niet inschuifbare stok (zie Figuur 3-7). Bij de bepaling van de grenslaag renslaag tussen water en slib wordt een rooster bevestigd op deze stok. De bepaling van de grens tussen vaste bodem en slib gebeurt met een smallere voet of de achterzijde van een prikstok. Naargelang de locatie, varieert men de grootte van het rooster. (STOWA, 2006). De VMM maakt gebruik van een uitschuifbare stok (zie Figuur 3-6),, omdat deze minder arbeidsintensief is. Het is tijdsrovend om tijdens de inmeting steeds het rooster te vervangen. ver De mate van uitschuiving is afhankelijk van de locatie/diepte. Voor zowel de bepaling van de sliblaag als de vaste bodem, maakt men gebruik van een smalle pin.

Figuur 3-6 Uitschuifbare stok

De werkwijze wijze van de peilstok is onafhankelijk van de uitvoeringsvorm. In beide gevallen wordt de stok loodrecht in het water gelaten tot wanneer men weerstand waarneemt. Op dit moment oment wordt de positie bepaald. Naargelang de locatie, varieert men de grootte van de voet. Indien dien de overgang tussen slib en vaste bodem moeilijk voelbaar is, maakt men soms gebruik van een aantal boringen (STOWA, 2006). De prikstokken worden zo licht mogelijk gemaakt zodat het gewicht van de stok geen invloed invloe heeft op de indrukking in de bodem.

13

Figuur 3-7 Niet inschuifbare peilstok

14

3.2.1.2 Peilhengel Het principe van de peilhengel is vergelijkbaar met dat van de peilstok, al heeft de peilhengel het grote voordeel dat at men vanaf de oever of vanuit een boot kan inmeten. De peilhengel kan allerlei grootten ootten aannemen met als maximaal bereik 8 m. Inmeten vanaf een boot is noodzakelijk wanneer de watergang te breed of slecht begaanbaar is (STOWA, 2006). Dit type van inpeilen wordt niet gebruikt door de VMM. De peilhengel is ideaal voor het inmeten van de bovenzijde van de sedimentlaag. De positiebepaling van de vaste bodem daarentegen is niet zo aan de orde, aangezien men heel goed de weerstandsverschillen erschillen in de bodem moet voelen om nauwkeurig te meten. Verder is het meten met een peilhengel veel arbeidsintensiever naarmate men verder van de kant of boot meet (Kennis Centrum Waterbodem, 2006). 2006) Volgende aspecten zijn heel eel belangrijk bij het gebruik van deze methode (STOWA, 2006): 2006) • • • •

loodrechte oodrechte stand van de peilhengel; gevoel van de surveyor voor de weerstand die overeenkomt met het niveau van bovenzijde en onderzijde van de sedimentlaag; onnauwkeurigheid keurigheid van het plaatsbepalingsysteem; …

Figuur 3-8 Peilhengel (STOWA, 2006)

15

3.2.2

Akoestische technieken

Onder de akoestische technieken verstaan we het singlebeam -en multibeam echolood, side scan sonar en het subbottom profiler systeem. Enkel singlebeam wordt door de VMM gebruikt.

3.2.2.1

Single - en multibeam echolood

Deze technieken worden bathymetrische technieken genoemd. Bij deze technieken wordt vanuit een boot, met behulp van geluidspulsen, de afstand tot de bodem bepaald (zie Figuur 3-9). De geluidsgolven worden uitgezonden, weerkaatsen op de bodem, keren terug en worden gedetecteerd. De looptijd van de echo is recht evenredig met de afstand. Uiteraard hangt de looptijd van de echo af van de geluidssnelheid in water. De geluidssnelheid is afhankelijk van de temperatuur, saliniteit en waterdiepte (Magelas, 2006). d=

vt 2

d= diepte v= geluidssnelheid t= aankomsttijd De akoestische technieken worden gebruikt op plaatsen waar de waterloop zich hiertoe leent. Vanaf een waterdiepte van 150 cm en een breedte van 3 m is het gebruik van bathymetrie aan te raden. Hierdoor is het logisch dat deze techniek enkel gebruikt wordt voor de inmeting van bevaarbare waterlopen of recreatievijvers. Men maakt gebruik van verschillende frequenties naargelang het te meten onderdeel van de waterloop. Voor het inmeten van de bovenzijde van de sliblaag maakt men gebruik van 210 of 700 kHz. Het bepalen van de onderzijde van de sliblaag gebeurt bij 33 kHz (Magelas, 2006). Het cruciale verschil tussen singlebeam en multibeam is het aantal uitgezonden bundels geluidsgolven. In tegenstelling tot multibeam wordt met singlebeam slechts één bundel uitgezonden loodrecht onder de boot. Multibeam kan niet gebruikt worden in waterlopen, aangezien de taluds voor interferentie zorgen. Deze methode is ideaal voor het inmeten van zeeën. Bij beide methoden is het gebruik van een geluidsmeter noodzakelijk (Magelas, 2006). De VMM maakt, indien de locatie zich ertoe leent, gebruik van de singlebeam methode.

16

Figuur 3-9 Singlebeam en Multib beam (STOWA, 2006)

17

3.2.2.2

Side Scan Sonar

Het Side Scan Sonar systeem wordt veel minder gebruikt dan single en multibeam echolood. Dit systeem bestaat uit een processor, een elektronische transmissiekabel en sleepkabel en uit een “towfish”. owfish”. Deze towfish is vergelijkbaar met de transducer van singlebeam en zendt (en ontvangt) akoestische signalen uit. Het systeem zweeft boven de waterbodem en zendt zijwaarts een akoestisch signaal uit. Dit signaal wordt door de bodem gereflecteerd en terug opgevangen door de towfish (zie Figuur 3-10). De towfish genereert per keer één energiepuls en wacht op de gereflecteerde golf. Hoe meer pulsen uitgezonden worden, hoe nauwkeuriger het beeld. Aangezien harde objecten sterker reflecteren dan zachte objecten wordt de Side Scan Sonar vooral gebruikt voor het detecteren van grof vuil (fietsen, autowrakken, …). …)

Figuur 3-10 Werking Side Scan sonar (Metaldec Survey BV, 2010)

3.2.2.3 Subbottom profiler Een laatste akoestische techniek is de subbottom profiler (Figuur 3-11),, deze wordt gebruikt voor de detectie van leidingen en objecten en het positioneren van van het grensvlak tussen de sliblaag en de harde bodem. Deze techniek is vergelijkbaar met singlebeam inglebeam met als enige verschil dat men een geluidspuls van lagere frequentie uitzendt (STOWA, 2006). 2006)

Figuur 3-11 Subbottom profiler (GSE Rentals, 2010)

18

3.2.2.4 Foutenbronnen De grootste onnauwkeurigheden bij het gebruik van akoestische methoden zijn (STOWA, 2006): • • • • •

Plaatsbepaling Diepgang Bootgeometrie Squat (diepgang van de boot) Interpolatie

Aangezien een sterkere reflectie zichtbaar is tussen vuil en sediment, ondervinden de akoestische technieken geen hinder van grof vuil op de bodem. Bij deze technieken is het zeer belangrijk om steeds de geluidssnelheid te meten omdat deze verschilt van waterlaag tot waterlaag. Deze is afhankelijk van de temperatuur en saliniteit van het water (Magelas, 2006).

3.2.3

Elektromagnetische techniek

Deze methoden worden niet gebruikt door d de VMM.

3.2.3.1 Georadar Georadar, ook wel grondradar genoemd, wordt gebruikt voor het peilen van watergangen. Dit systeem zendt, via een zendantenne, kort magnetische pulsen uit. Deze pulsen reflecteren op inhomogeniteiten in de bodem. Deze methode wordt gebruikt gebruikt in ondiepe (tot 3.5 m) zoete watergangen. Brakke of zoute waterlopen zijn ongeschikt (STOWA, 2006). 2006) Georadar wordt vooral gebruikt voor het onderzoeken van waterbodems op de aanwezigheid van leidingen, kabels e.d. Deze techniek techniek is niet geschikt voor alle waterlopen, maar is volop in ontwikkeling voor het meten van waterdieptes.

Figuur 3-12 Georadar (Prospección Geofísica y Geofísica Pura, 2010)

19

3.2.3.2 Waterdruk techniek Een bekende waterdruk techniek is het Medusa systeem. Dit systeem maakt gebruik van de waterdruk voor het bepalen van de waterdiepte. Een dieptesensor sleept achter de boot over de bodem van de waterloop. De diepte wordt bepaald door rekening te houden met het feit dat de dichtheid bij de overgang van water naar waterbodem tot op een bepaalde waarde gaat. De hoogteligging van de overgang wordt verder afgeleid uit een waterdrukmeter die de hoogte van de waterkolom bepaalt. bepaalt Het Medusa systeem steem kan op alle bodemsoorten bodemsoorten toegepast worden. In tegenstelling tot de akoestische technieken moet de geluidssnelheid niet gemeten worden.

Figuur 3-13 Waterdruktechniek (Kennis Centrum Waterbodem, 2006)

3.2.4

Positiebepaling

Zonder nauwkeurige positiebepaling zijn de inmetingen waardeloos. Hieronder volgt een opsomming van de verschillende mogelijke methoden. Het is zeer belangrijk dat tijdens de kwantiteitsmetingen gebruik gemaakt wordt w van een nauwkeurige plaatsbepaling. De gebruikte techniek is volledig afhankelijk van locatie tot locatie. Volgende methoden zijn mogelijk (STOWA, 2006): •

Uitzetwerk (meetlint en/of piketten)

•

(D)GPS

•

RTK-GPS

•

Totaalstation

•

Laser

20

3.2.4.1

Uitzetwerk (meetlint en/of piketten)

Hierbij wordt de locatie van een punt bepaald d.m.v. het totaalstation.. De D exacte locatie wordt opgeslagen en met behulp van een nagel of piket wordt de plaats aangeduid. Vanaf dit punt kan men verdere locaties uitzetten. uitz

Figuur 3-14 Uitzetwerk

21

3.2.4.2 (D)GPS GPS staat voor Global Positioning System en werd ontwikkeld door het Departement of Defense (DoD) van de USA. Het systeem bestaat uit minstens 24 satellieten op een hoogte van 20 000 km. Elke satelliet draait tweemaal per dag om de aarde in 6 verschillende banen (Magelas). Iedere satelliet zendt een uniek GPS signaal uit. Dit signaal bestaat uit vier soorten informatie: (Magelas, 2006) •

C/A code en een P code. bevatten de tijdinformatie die wordt gegenereerd door de atoom klokken in de satelliet. De Coarse Aqcuisition code (C/A code) bevat de identificatie code van de satelliet. Deze code wordt elke milliseconde uitgezonden. De P-code is enkel decodeerbaar door het Amerikaanse leger.

•

Almanak data: levert de GPS ontvanger informatie over de plaats van de satelliet op een bepaald tijdstip (deze informatie is op voorhand bekend)

•

Ephemeris data: bevat actuele informatie over de status van satelliet, alsook de datum en de tijd.

De almanak en de ephemeris vormen samen de 'navigation message', de NM. Het principe is gebaseerd op de afstandsmeting tussen de satelliet en de ontvanger. De GPS ontvanger berekent de afstand tot de satellieten door de tijd te meten waarin de GPS signalen de ontvanger bereiken. M.a.w. het tijdsverschil tussen het verzonden en ontvangen signaal wordt gemeten (Van Boxstael, 2009). De ontvanger synchroniseert de interne klok met de kloksignalen van de ontvangen GPS signalen. De ontvanger keert als het ware terug in de tijd, naar de tijd van uitzending door de satelliet. Tijdens de synchronisatie wordt aan de hand van de C/A code de exacte positie van de satelliet bepaald door deze op te zoeken in de ephemeris gegevens uit de NM (Van Boxstael, 2009). Hoeveel tijd de ontvanger moet worden teruggeschoven, is de reistijd van het GPS signaal. Uit dit tijdsignaal kan de afstand tot de satelliet bepaald worden. Het tijdsverschil wordt vermenigvuldigd met de snelheid van het signaal, de lichtsnelheid (Van Boxstael, 2009). Eén satelliet geeft één afstand die overeenkomt met een sfeer. Indien men drie sferen laat kruisen, verkrijgt men twee mogelijke plaatsbepalingen. Soms is één punt verwerpbaar. Pas vanaf vier satellieten, kan exact de plaatsbepaling gebeuren door de vier sferen te laten snijden. Waar de sferen snijden kan ook de pseudorange genoemd worden (Magelas, 2006). De nauwkeurigheid van de GPS is onderhevig aan een aantal foutenbronnen: •

Vertragingen in ionosfeer en troposfeer;

•

multipath: vertraging van signaal op grote objecten;

•

onnauwkeurigheid van de klok in de ontvanger;

•

orbitale fouten: onnauwkeurigheden in de posities van de satellieten;

22

•

aantal satellieten (Magelas, 2006).

Deze foutenbronnen worden opgelost door DGPS. DGPS staat voor Differentieel Global Positioning System. Het principe van DGPS bestaat eruit dat een (referentie) GPS ontvanger op een bekende plaats wordt gepositioneerd. Met het ontvangen signaal wordt een omgekeerde berekening uitgevoerd waardoor de ontvangsttijden bepaald worden aan de hand van de gekende positie. Vervolgens vergelijkt het systeem de normale ontvangen reistijden met de berekende reistijden, waarna een correctiefactor wordt bepaald. De referentieontvanger stuurt zijn correctiefactor door via radiolink, satellietsignaal of telefoonverbinding naar de gewone ontvangers. Hierdoor worden de metingen van de gewone ontvangers gecorrigeerd (Van Boxstael, 2009).

3.2.4.3

RTK-GPS

De werking van RTK ligt in het verlengde van DGPS. Toch is er een groot verschil tussen beide: DGPS werkt met een gecodeerd signaal en RTK met de faseverschillen in de satellietsignalen. Bij RTK is het satellietsignaal zelf (C/A code, …) niet interessant, wel de fase van de draaggolf van het signaal. I.t.t. DGPS is RTK tot op enkele cm nauwkeurig (Van Boxstael, 2009). RTK maakt gebruik van een referentieontvanger en een mobiele ontvanger (rover) waarbij het faseverschil gemeten wordt tussen het uitgestuurde signaal van de satelliet en het door de ontvangers gegenereerde identieke signaal. De referentieontvanger bevindt zich op een plaats met een nauwkeurige positiebepaling. Met een geheel aantal golflengten en de faseverschillen wordt d.m.v. de afstanden naar de satellieten een basislijn berekend tussen de twee ontvangers (Van Boxstael, 2009). Aan de hand van de basislijn wordt de positie van de mobiele ontvanger precies gemeten. De afstanden worden berekend door volgende vergelijking: ∆X=Xsatሺtሻ- Xrec

∆X =

Afstand tussen satelliet en ontvanger.

Xsat (t ) =

X-coördinaat van de satelliet op het tijdstip van uitzenden.

Xrec =

X-coördinaat van de ontvanger.

Deze vergelijking wordt eveneens gebruikt voor het bepalen van de Y, Z coördinaten. De afstand van verschillende satellieten tot de ontvangers wordt gelijktijdig gemeten, waardoor een verschilvector kan opgesteld worden tussen beide ontvangers. Door de vector te koppelen aan de referentieontvanger, wordt een nauwkeurige plaatsbepaling van de rover verkregen (Van Boxstael, 2009).

23

3.2.4.4 Totaalstation Een totaalstation of tachymeter bestaat uit een theodoliet en een afstandsmeter. Een theodoliet meet horizontale en verticale verticale hoeken doordat een kijker om een verticale en een horizontale as draait. Door een meetsysteem op beide assen wordt de hoekgrootte afgelezen. Het is zeer belangrijk dat het toestel waterpas staat. Dit is mogelijk door een ingebouwde ‘waterpas’ en verstelbare rstelbare schroeven. De afstandsmeter bepaalt de schuine afstand tot een punt. Deze schuine afstand wordt omgerekend naar de horizontale afstand. D.m.v. de horizontale hoek en de horizontale afstand worden de (x, (x y)- coördinaten van een bepaald punt bepaald. d. Bij betrekking van de verticale hoek bekomt men de (x, (x y, z) coördinaten, waarbij z de hoogte weergeeft (Van Boxstael, 2009).

Figuur 3-15 Totaalstation

24

3.2.4.5

Laser

Laser is een vrij rij nieuwe plaatsbepalingsmethode. Via laser altimeters of laserscanners is het mogelijk de plaats nauwkeurig te bepalen. Deze techniek wordt ook wel laser altimetry of LIDAR genoemd. De laser altimeter wordt in een vliegtuig of satelliet geplaatst en stuurt stu een infrarood lasersignaal uit op de in te meten locatie. Het signaal wordt door de aarde gereflecteerd en terug opgevangen door de laser altimeter. Met behulp van de looptijd van het h signaal en rekening houdend met de geluidssnelheid, wordt de afstand tot de aarde bepaald. Uiteraard is het de bedoeling dat de (x, (x y, z)- coördinaten berekend erekend worden. Naast de afstand tot het oppervlak moet men de locatie van het vliegtuig van waaruit de laser vertrok kennen, alsook de richting van de laseraltimeter. Deze Deze gegevens worden verkregen via GPS (GeoLas Consulting, 2010). In de tachymeter van de VMM zit een laserfunctie die plaatsbepaling mogelijk maakt. Deze functie is vrij langzaam in gebruik en wordt niet toegepast bij de VMM.

Figuur 3-16 Laser altimetrie (AGIV, 2011)

25

3.2.4.6

Conclusie

Voor een hoge accuraatheid is het gebruik van RTK-GPS (zie Tabel 3-1) aan te raden. Deze techniek is toepasbaar op zowel handmatige als akoestische technieken. Jammer genoeg Is het gebruik van GPS niet altijd mogelijk door te veel begroeiing, geen ontvangst, bruggen, … Als dit het geval is, gaat men over tot het gebruik van het totaalstation. Bij deze methode wordt een spiegel bevestigd op de peilstok, peilhengel, boot …. Het nadeel van het totaalstation is dat deze methode vrij arbeidsintensief is. Tabel 3-1 Overzicht plaatsbepalingsmethodes i.f.v. meettechnieken (STOWA, 2006)

Meettechniek

Peilstok

Peilhengel

Singlebeam

Multibeam

Georadar

Medusa

Uitzetwerk

5

5

n.v.t.

n.v.t.

n.v.t.

n.v.t.

GPS

5m

n.v.t.

5m

5m

5m

5m

RTK-GPS

1 cm

n.v.t.

1 cm

1 cm

1 cm

1 cm

(D)GPS

2m

n.v.t.

2m

2m

2m

2m

Totaalstation

1 cm

1 cm

n.v.t.

n.v.t.

n.v.t.

n.v.t.

Laser

1 cm

1 cm

n.v.t.

n.v.t.

n.v.t.

n.v.t.

26

3.2.5

Keuze meetmethoden

Zoals eerder aangehaald hangt het gebruik van een bepaalde methode af van de locatie. Indien men gebruik maakt van een boot, zoals bij de akoestische methoden, moet men er rekening mee houden dat een waterdiepte van 1.5 m noodzakelijk is. Verder kan men handmatig inpeilen wanneer ook hier de waterdiepte dit toelaat. Tabel 3-2 Overzicht gebruik van meetmethoden (STOWA, 2006)

Meettechniek

Peilstok

Peilhengel

Singlebeam

Multibeam

Georadar

Medusa

Waterdiepte

<1.5 m

<2 m

>1.5 m

>1.5 m

0.3 – 0.8 m

>0.2 m

Waterbreedte

onbeperkt

<8 m

>3 m

>3 m

>3 m

onbeperkt

Bovenkant slib

+

+

+

+

+/-

+

Onderkant slib

+/-

+/-

+/-

-

?

-

Zout/brak water

++

++

++

++

n.v.t.

++

Reactie bodemvreemd materiaal

+

+

+

+

++

+

3.2.6

Spreiding meetmethoden

De metingen moeten zo worden uitgevoerd dat een interpolatie tussen de verschillende punten een accuraat beeld geeft van de werkelijkheid. Zo zal men bij het gebruik van handmatig inpeilen, singlebeam echolood en georadar gebruik maken van dwarsraaien. De afstand tussen 2 dwarsraaien is van volgende aspecten afhankelijk: •

Type waterloop: in kunstmatig rechte aangelegde waterlopen en meren varieert het bodemprofiel minder sterk dan in bochtige natuurlijke beken.

•

Type gebied: in stedelijk gebied bevatten de waterlopen meer grof vuil en verhardingen dan in landelijk gebied.

•

Doel peiling: de inpeilingen kunnen nodig zijn voor onderhoudsbaggerwerken of voor saneringsbaggerwerken. In het geval van saneringsbaggerwerken is de nauwkeurigheid van de metingen van primordiaal belang, gezien men veelal te maken heeft met verontreinigd slib en de baggerkosten veel hoger liggen.

De afstand tussen 2 dwarsraaien varieert tussen 25-50m bij saneringsbaggerwerken en om de 200 m bij onderhoudsbaggerwerken in landelijk gebied. In beken werkt men sowieso met een afstand van 25-50 m, gezien de onvoorspelbaarheid van deze waterlopen.

27

3.3 3.3.1

Dataverwerking Algemeen

Na het inmeten van de gekozen waterloop, worden data verkregen in de vorm van (x, y, z)punten. Deze punten ontstaan door inmeten met behulp van één van de besproken methodes op pagina 13. Bij verder onderzoek van de waterloop is omzetting van deze punten in een 3D- model wenselijk. Een 3D-model, of een Digital Terrain Model, is een digitale weergave van het ingemeten oppervlak. Soms wordt een DTM ook wel DEM (Digital Elevation Model) genoemd (Droj, 2008). DTM’s worden gebruikt voor volgende toepassingen: • • • • • • •

Visualisatie van het terrein; reductie van zwaartekrachtmetingen; cartografische en morfologische terrein analyse; verbeteren van satellietfoto’s; extractie van terreinparameters onderzoek naar erosie; etc.

Verder worden DTM’s veel gebruikt om het reliëf van terreinen weer te geven in GIS systemen (Dalyot, Doytsher, & Katzil, 2009). De nauwkeurigheid van het DTM hangt af van interpolatiefouten, gekozen resolutie en de nauwkeurigheid van de bron data. Uiteraard wordt de nauwkeurigheid van de databron beïnvloed door de gekozen positiebepaling. (RTK-GPS en totaalstation zijn accurater dan uitzetwerk) (Keim et al., 1999) (Dadson, 1999). Een DTM kan uit verschillende soorten databronnen gegenereerd worden. De voornaamste vormen zijn fotometrisch (LIDAR), terreinsurveys en digitale cartografische data (Rahman, 1994). Digitale cartografische data worden veel gebruikt gezien de grote beschikbaarheid aan topografische kaarten. Deze vorm van dataverwerking is niet mogelijk bij onderzoek van subaquatische terreinen. In deze cases maakt men veelal gebruik van terreinsurveys. DTM kaarten kunnen in verschillende ‘structuren’ weergegeven worden, de voornaamste zijn: • • •

TIN (Triangular Irregular Network) Contouren Grid

Het TIN werd reeds in de jaren 70 ontwikkeld met de bedoeling een netwerk te vormen van onregelmatig geplaatste punten (Wadembere, 2000).Het is een netwerk van vectoren waarbij de gevormde driehoeken elkaar niet overlappen. Deze driehoeken variëren door de terreinruwheid (EGIS foundation, 1994). Een volgende weergave voorstelling van een DTM zijn contouren. Deze worden vooral gebruikt bij topografische weergaven. Contouren verbinden de punten met gelijke hoogte met elkaar.

28

De laatste mogelijke weergave is een grid, dit is gebaseerd op een regelmatig re verdeeld rooster. Een grid kan rechtstreeks vanuit een puntenbestand gegenereerd worden, maar conversie van TIN of contouren naar grid is ook mogelijk.

Figuur 3-17 Van links naar rechts: contour, grid, g TIN (Dalyot, Doytsher, & Katzil, 2009)

De echte berekening van deze modellen gebeurt door spatial interpolation algorithms (Droj, 2008). Elk van de eerder opgesomde weergaven werkt met een specifiek sp interpolatie algoritme (zie Figuur 3-17).

3.3.2

Terramodel

In deze masterproef wordt gebruik gemaakt van Terramodel. Dit is een softwarepakket geproduceerd door Trimble, een Amerikaans bedrijf gespecialiseerd in positiebepalingsystemen. Terramodel wordt omschreven als een flexibel software pakket, geschikt voor de snel veranderende wereld van survey, ontwerp en constructie (Trimble, 2010).. Gezien VMM gebruik maakt van Trimble materiaall (GPS, totaalstation, totaalst …) is het logisch dat de keuze op Terramodel viel. Het verwerkingsschema van Terramodel ziet er als volgt uit (zie Figuur 3-18 18):

29

Figuur 3-18 Verwerking via Terramodel (Rijkswaterstaat directie Limburg, 2003)

Hierbij worden de (x, y, z)-punten punten omgezet in een continu vorm/model. Bij deze conversie wordt gebruik gemaakt van interpolatie. In het geval van Terramodel gaat men over tot een DTM in de vorm van een TIN. Een eerste belangrijke stap is het vormen van het TIN (Wadembere, 2000). Hierbij is het belangrijk dat bij de inmeting, punten zo werden gekozen dat elke (belangrijke) verandering verandering in het oppervlak vertegenwoordigd wordt door een punt. Verder is de vorming van de driehoeken zeer belangrijk. Deze kan op verschillende manieren gebeuren: • •

The Delaunay triangulatie (Leica, 2010) Ordenen volgens afstand (Klinkenberg, 1997)

In Terramodel wordt gebruikt gemaakt van de Delaunay triangulatie.

3.3.2.1 The delaunay triangulatie De vorming van driehoeken is gebaseerd op de Delaunay triangulatie. Deze methode werd ontdekt door Boris Delaunay in i 1934 (McEachron, 2009).. Hij gaf volgende definitie aan zijn methode: de omgeschreven cirkel van een driehoek wordt gevormd door drie punten en kan geen andere punten bevatten dan de drie die de cirkel bepalen. De omgeschreven cirkel wordt verkregen door de middelloodlijnen middelloodlijnen van twee zijden te tekenen en met het verkregen snijpunt de cirkel te construeren (zie Figuur 3-19).

30

Figuur 3-19 Omgeschreven mgeschreven cirkel

Het vormen van Delaunay Triangles berust op algoritmische berekeningen. De vijf meest gebruikte algoritmes zijn Incremental, Gift Wrap, Divide and Conquer, Sweepline en QuickHull (Su & Drysdale, 1996). 1996) Van deze ze vijf algoritmen zijn Divide and Conquer en Incremental het meest aangewezen voor de constructie con van een TIN (De Wulf, Hennau, & Constales, 2006). Een veel gebruikte methode voor het vormen van Delaunay Triangles is edge flipping. fli Edge flipping, kortweg flipping genoemd, wordt gebruikt wanneer de driehoeken niet voldoen aan het Delaunay principe. Bij het flippen van zijden is het de bedoeling dat twee driehoeken die niet voldoen aan de Delaunay driehoeken na de ingreep wel voldoen. vo Deze werkwijze wordt weergegeven in Figuur 3-20.. De driehoek abc en acd voldoen niet aan Delaunay Triangles aangezien de omschrijvende cirkels respectievelijk d en b bevatten. Na ‘flipping’ zijde ac tot bd, worden de driehoeken abd en bcd gecreëerd die wel voldoen aan Delaunay (De Wulf, Hennau, & Constales, 2006).

Figuur 3-20 Edge flipping (De Wulf, Hennau, nnau, & Constales, 2006)

Divide and Conquer algoritme Het divide and Conquer onquer algoritme gebruikt de hele dataset en verdeelt deze in subsets zodat er één of twee driehoeken naast elkaar overblijven. Hierna worden de deeloplossingen aan elkaar verbonden (zie Figuur 3-21). 3 Deze methode is t.o.v. het Incremental algoritme het meest aangewezen om grote datasets te verwerken erken gezien de hogere snelheid (De Wulf, Hennau, & Constales, 2006). 2006)

31

Figuur 3-21 Divide and Conquer (De Wulf, Hennau, & Constales, 2006)

Incremental algoritme Dit algoritme start met de constructie van een driehoek die alle punten van de dataset bevat. Eén voor één worden de punten toegevoegd aan de driehoek en komen de Delaunay driehoeken tot stand. Dit proces wordt weergegeven in Figuur 3-22. Eerst wordt het nieuwe punt P toegevoegd (a). Vervolgens wordt het punt verbonden met de zijden van de bestaande driehoek waarin het punt zich bevindt. Hierna wordt gecheckt of de nieuwe verbinding nog steeds voldoet aan het Delaunay principe. Er mag zich, met andere woorden, geen nieuw punt in de omgeschreven cirkel bevinden. Als de gevormde driehoek hieraan voldoet, verandert er niets en werd een deel van de Delaunay driehoek gevormd (c, d). Indien, zoals in (d), de gevormde driehoek niet voldoet, wordt de zijde ‘geflipt’ en worden 2 andere diagonalen kandidaat om een zijde te worden. Dit proces gaat verder tot er geen kandidaat zijden meer te vinden zijn (streepjeslijn). (Lischinski, 1993) Doordat de snelheid daalt naarmate de dataset stijgt, is dit algoritme minder geschikt voor de constructie van een TIN uit grote datasets. Dit algoritme is eerder geschikt voor het updaten van een TIN en niet voor het vormen van een TIN (De Wulf, Hennau, & Constales, 2006).

32

Figuur 3-22 Toevoeging van een punt in de driehoek. (Lischinski, 1993)

3.3.2.2 Interpolatie Na het construeren van het TIN gaat Terramodel over op interpolatie. Aangezien tijdens het inmeten een dataset ontstond die niet de gehele oppervlakte oppervlakte bedekt, is interpolatie nodig om alle tussenliggende punten te ‘voorspellen’ (zie Figuur 3-23).. Interpolatie wordt onderverdeeld in globale en lokale interpolatie. In het eerste geval is er sprake van een interpolatie nterpolatie die van alle punten in de dataset gebruik maakt. Logischerwijs zal het veranderen van één punt, de hele berekening veranderen (Caverly).. Lokale interpolatie maakt enkel gebruik van de omringende punten, waardoor bij verandering van één punt enkel de omringende punten beïnvloed worden (Fürst) (Caverly).. De interpolatiekeuze wordt bepaald door het type variabele en de afstand tussen de punten.

33

z z2 z3

z0 ?

z1

y (x 3,y3) (x0 ,y0 )

(x 1,y1)

(x 2,y2)

x Figuur 3-23 Interpolatie (Fürst)

Volgende interpolatietechnieken worden veelvuldig gebruikt: • • • • •

Trend IDW (Inverse Distanced Weighted method) Kriging Lineaire TIN interpolatie …

Terramodel maakt gebruikk van een lineaire TIN interpolatie. interpolatie

34

Lineaire TIN interpolatie

Figuur 3-24 TIN Interpolatie (Hemphill, 2009)

Eens een TIN gevormd werd, erd, zijn alle gekende punten lineair verbonden. verbonden. Hieruit wordt een driehoek verkregen die een vergelijking vergelijk vormt van volgende vorm: Ax + By + Cz + D = 0

35

Via deze vergelijking is het mogelijk om de onbekende punten in elke driehoek te bepalen. A, B, C en D zijn gekende punten die bepaald worden door de (x, y, z)- coördinaten bekomen uit de opmeting. Veronderstel dat de 3 gekende punten er als volgt uit zien: (x1, y1, z1), (x2, y2, z2), en (x3, y3, z3) (EMS-i, 2011). De punten A, B, C en D worden als volgt bepaald: A = y1 (z2 - z3) + y2 (z3 - z1) + y3 (z1 - z2) B = z1(x2 - x3) + z2(x3 - x1) + z3(x1 - x2) C = x1(y2 - y3) + x2(y3 - y1) + x3(y1 - y2) D = -Ax1 - By1 - Cz1 Hierdoor kan de vergelijking herschreven worden als volgt: z=fሺx,yሻ=-

A B D x- yC C C

Met deze vergelijking kan men voor gelijk welk punt in het TIN de hoogte bepalen (EMS-i, 2011).

36

3.4

Groot Schijn

Het Groot Schijn is een waterloop van eerste categorie gelegen in de Benedenschelde (zie Figuur 3-25). Het oorsprongsgebied bevindt zich in Westmalle en stroomt doorheen Zoersel, Oelegem, Schilde, Wijnegem, Wommelgem, Deurne en Borgerhout. De Schijnmonding werd in de loop der jaren meermaals verlegd voor de Antwerpse Haven. De rivier was tot 2009 overwelfd in het et havengebied nabij stortplaats de Hoge Maey. Maey. Het Groot Schijn wordt met een vijzelgemaal naar het Lobroekdok gepompt om waterellende in Merksem te voorkomen en vismigratie mogelijk te maken. Het onderzoek naar de efficiëntste inmetingsmethode om typische Vlaamse bochten in te meten, werd uitgevoerd aan de hand van het Groot Schijn. Deze waterloop is enorm bochtig, waardoor de keuze snel gemaakt was.

Figuur 3-25 Situering Groot Schijn

37

4 Materiaal en methoden 4.1 4.1.1

Onderzoek van Vlaamse bochten Inventarisatie materiaal

Vervoer • • • • • •

Boot o Roeispanen o Meertouw Trailer Auto Buitenboordmotor Benzinebak Benzine

Inpeiling • •

• • •

Uitschuifbare peilstok GPS o Antenne o COAX-kabel (10m) o octopuskabel TSC2 o Antenne o Markeerstokje GSM Waadpak

Verwerking •

Computer met Terramodel

38

4.1.2

Methode

Het onderzoek naar de meest efficiënte methode voor het inmeten van verschillende types van meandering (bochten) van Vlaamse waterlopen werd uitgevoerd aan de hand van het Groot Schijn. Dit is een waterloop van eerste categorie die behoorlijk kronkelt en bochten bevat die veel voorkomen in de Vlaamse waterlopen. waterlopen Vijf typische bochten werden uitgekozen en ingemeten. De karakterisatie van deze de bochten gebeurt door het aanduiden duiden van de graadsverandering die de waterloop water maakt. Als plaatsbepalingsmethode werd voor het gebruik van GPS gekozen. In dit geval werd gebruik gemaakt van een Trimble SPS751 GPS en een antenne met modelnaam Zephyr. De antenne werd bevestigd op de uitschuifbare prikstok. pri Via een COAX- kabel werd de antenne met de GPS verbonden. Deze werd op zijn beurt verbonden met de TSC2. Dit is een innovatieve zakcomputer die alle ingemeten data opslaat. Aangezien het waterpeil tijdens de inmeting hoog stond was waden niet mogelijk en werd gebruik gemaakt van de boot. Met behulp van de peilstok werd telkens de positie van het topslib en de vaste bodem bepaald en opgeslagen in de TSC2. Deze posities werden bepaald door door het verschil in weerstand die de surveyor waarneemt. Dit is puur intuïtief en persoonsgebonden. persoonsgebonden

Figuur 4-1 Inmeting

Elke bocht werd met een maximaal aantal profielen ingemeten. Dit beeld is het meest accuraat en wordt als waarheidsgetrouw beschouwd. Uit deze data werd een DTM gemaakt, alsook het volume lume sediment werd berekend. Vervolgens werd op basis van het resultaat steeds het minst relevante profiel verwijderd. Hieruit werd terug een DTM en volumeberekening gemaakt. Dit proces werd herhaald tot een minimaal aantal profielen werd bekomen, die nog steeds beantwoorden aan de oorspronkelijke oorspronkelijk vorm van de bocht. In bijlage 1 wordt een gedetailleerd overzicht gegeven van de gebruikte procedure.

39

4.2 4.2.1

Smoothing Inventarisatie materiaal

Vervoer • • • • • •

Boot o Roeispanen o Meertouw Trailer Auto Buitenboordmotor Benzinebak Benzine

Inpeiling • •

• • •

Uitschuifbare peilstok GPS o Antenne o COAX-kabel (10m) o octopuskabel TSC2 o Antenne o Markeerstokje GSM Waadpak

Verwerking •

Computer met Terramodel

40

4.2.2

Methode

Het onderzoek naar een passende ‘smoothing’ is noodzakelijk met het oog op de uit te voeren baggerwerken. Na het inpeilen en e de verwerking van de data in Terramodel, erramodel, ontstaat een DTM. Uit dit model worden profielen gegenereerd waaruit de aannemer, die de baggerwerken zal uitvoeren, akte krijgt van de toestand in de waterloop. Hieruit kan de aannemer nnemer afleiden op welke plaatsen welke hoeveelheden sediment liggen en waar hij logischer wijs zal moeten baggeren. bagger In Figuur 4-2 wordt het onderzochte traject weergegeven in een isopach of verschilkaart. Deze kaart wordt gelezen als een reliëfkaart. Enerzijds geven de rode gedeelten gedeelte de plaatsen met de dikste sedimentlaag weer en anderzijds duiden de donkerblauwe gedeelten gedeelte de plaatsen met de dunste sedimentlaag aan. Uiteraard worden de tussenliggende sedimentdikten voorgesteld door door de tussenliggende kleuren (zie legende). De VMM ondervond dat de aannemers problemen ondervonden bij het lezen van de ruwe profielen. Het is onmogelijk om met graafmachines de werkelijke ‘grillen’ van de waterloop uit te graven waardoor er vraag ag kwam naar naa meer gelijkmatiger// uitgemiddelde/gesmoothe profielen. Het onderzoek naar een passende ‘smoothing’ werd behandeld met behulp van een deeltraject van het Groot Schijn. De methode van inpeilen staat beschreven in 4.1.2.

Figuur 4-2 Onderzochte traject

Als plaatsbepalingsmethode werd voor het gebruik van GPS gekozen. In dit geval werd gebruik gemaakt van een Trimble SPS751 GPS en een antenne met modelnaam Zephyr. De antenne werd bevestigd op de uitschuifbare prikstok. Via een COAXCOAX kabel werd de antenne met de GPS verbonden. Deze werd op zijn beurt verbonden met de TSC2. Dit is een innovatieve zakcomputer die alle ingemeten data opslaat. Aangezien het waterpeil tijdens tij de inmeting hoog stond was waden niet mogelijk en werd gebruik gemaakt van de boot. Met behulp van de peilstok werden en telkens de positie van het topslib en de vaste bodem bepaald en opgeslagen in de TSC2. Deze posities werden bepaald door door het verschil in weerstand die de surveyor waarneemt. Dit is puur intuïtief en persoonsgebonden.

41

Figuur 4-3 Inmeting

Na de inmeting werd een DTM gegenereerd, gegenereerd alsook werd het totale sedimentvolume van het traject berekend. Het smoothen van de data gebeurt door op het TIN een raster (grid) te leggen. Hoe groter het raster, hoe smoother de profielen. Dit valt te verklaren doordat men bij een grotere grid over een kleiner kleiner aantal punten uitmiddelt, waardoor de extremen worden worde afgevlakt.

Figuur 4-4 Raster grootte

42

Op het traject worden volgende rastergroottes vergeleken: • • • • • • • • •

0.25 x0.25 m 0.5 x 0.5 m 0.75 x 0.75 m 0.8 x 0.8 m 0.85 x 0.85 m 0.9 x 0.9 m 1x1m 1.25 x 1.25 m 1.5 x 1.5 m

De verschillende modellen worden visueel beoordeeld op het vlak van smoothing, maar ook het verschil in sedimentvolume t.o.v. de ruwe data wordt bekeken. Hieruit wordt het meest efficiënte raster bepaald.

43

5 Resultaten en conlusies 5.1

Onderzoek Vlaamse bochten

Tijdens de bespreking van de resultaten werd een foutenpercentage tot 5% getolereerd.

5.1.1

Bocht van 40°

In Tabel 5-1 worden de resultaten van een bocht van 40° weergegeven. In Figuur 5-2 wordt het foutenpercentage visueel voorgesteld. Tabel 5-1 Resultaten bocht van 40°

Aantal profielen

Slibvolume (m³)

5 4 3

123 123,6 110

Verschil t.o.v. Verschil t.o.v. max. max. aantal aantal profielen (m³) profielen (%) 0 0 -0,6 0,5 13,6 11

procentuele fout (%)

40° 12 11 10 9 8 7 6 5 4 3 2 1 0 3

4

5

Aantal profielen

Figuur 5-1 Visuele voorstelling foutenpercentage van een bocht van 40°

44

Figuur 5-2 Verschilkaart bocht van 40°, a (5 profielen), b (4 profielen), len), c (3 profielen)

Figuur 5-2 geeft de evolutie weer van de sedimentligging ntligging tijdens de variatie van het aantal profielen. Aangezien we hier te maken hebben met een verschilkaart wordt deze, deze zoals aangehaald in 4.2.2, gelezen zen als een reliëfkaart. Uit de resultaten kan besloten worden dat bij dit type bocht vier profielen de meest efficiënte manier van inmeten is. s. Een foutenpercentage van 0.5 % is verwaarloosbaar.

45

5.1.2

Bocht van 80°

In Tabel 5-2 worden de resultaten van een bocht van 80° weergegeven. In Figuur 5-3 wordt het foutenpercentage visueel voorgesteld. Tabel 5-2 Resultaten bocht van 80°

Aantal profielen

Slibvolume (m³)

8 7 6 5 4

137,8 137,3 138,2 144 145,9

Verschil t.o.v. max. aantal profielen (m³) 0 0,5 0,4 6,2 8,1

Verschil t.o.v. max. aantal profielen (%) 0 0,4 0,3 4,5 5,9

80° Procentuele fout (%)

7 6 5 4 3 2 1 0 4

5

6

7

8

Aantal profielen

Figuur 5-3 Visuele voorstelling foutenpercentage van een bocht van 80°

46

len), c (6 profielen), d (5 Figuur 5-4 Verschilkaart bocht van 80°, a (8 profielen), b (7 profielen), profielen), e (4 profielen) sedimentligging tijdens de variatie in het aantal Figuur 5-4 geeft de evolutie weer van de sedimentligging profielen. In dit geval zijn vijf profielen aangewezen. Het foutenpercentage van 4.5 % ligt binnen de toegelaten waarde en het inmeten is het minst arbeidsintensief.

5.1.3

Bocht van 105°

In Tabel 5-3 worden de resultaten van een bocht van 80° weergegeven. In Figuur 5-5 wordt het foutenpercentage visueel voorgesteld. Tabel 5-3 Resultaten bocht van 105°

Aantal profielen

Slibvolume (m³)

10 9 8 7 6

100.3 101.7 102.7 108.8 112

Verschil t.o.v. max. aantal profielen (m³) 0 1.4 2.4 8.5 11.7

Verschil t.o.v. max. aantal profielen (%) 0 1.3958 2.3928 8.4745 11.665

47

Procentuele fout (%)

105° 13 12 11 10 9 8 7 6 5 4 3 2 1 0 6

7

8

9

10

Aantal profielen

Figuur 5-5 Visuele voorstelling oorstelling foutenpercentage fouten van een bocht van 105°

Figuur 5-6 Verschilkaart bocht van 105°, a (10 profielen), b (9 profielen), c (8 profielen), d (7 profielen), e (6 profielen)

Figuur 5-6 geeft de evolutie weer van de sedimentligging tijdens de variatie in het aantal profielen. Indien we rekening houden met het maximaal toegelaten foutenpercentage van 5% is het noodzakelijk om van een bocht van 105° 8 profielen in te meten, ten einde e een zo goed mogelijke benadering van de werkelijkheid weer te geven.

48

5.1.4

Bocht van 120°

In Tabel 5-4 worden de resultaten van een bocht van 80° weergegeven. In Figuur 5-7 wordt het foutenpercentage visueel voorgesteld. Tabel 5-4 Resultaten bocht van 120°

Aantal profielen

Slibvolume (m³)

9 8 7 6 5 4

72.8 70.4 72.3 67.7 72.4 63

Verschil t.o.v. max. aantal profielen (m³) 0 2.4 0.5 5.1 0.4 9.8

Verschil t.o.v. max. aantal profielen (%) 0 3.3 0.7 7 0.5 13.5

Procentuele fout (%)

120° 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 4

5

6

7

8

9

Aantal profielen

Figuur 5-7 Visuele voorstelling foutenpercentage van een bocht van 120°

49

Figuur 5-8 Verschilkaart bochtt van 120°, a (9 profielen), b (8 profielen), len), c (7 profielen), d (6 ( profielen), e (5 profielen), f (4 profielen)

Figuur 5-8 geeft de evolutie weer van de sedimentligging tijdens de variatie in het aantal profielen. Bij deze bocht zien we een bijzonder fenomeen. Naarmate het aantal profielen vermindert zien we dat ookk de vorm van de bocht verandert. verandert Hierdoor wordt het volume verschil (m³) opgeheven door de verandering van de bocht door het weghalen van profielen (en aldus verandering ng van de waterlijn). In dit geval kunnen vijf profielen reeds volstaan voor het inmeten.

50

5.1.5

Bocht van 135°

In Tabel 5-5 worden de resultaten van een bocht van 80° weergegeven. In Figuur 5-9 wordt het foutenpercentage visueel voorgesteld. Tabel 5-5 Resultaten bocht van 135°

Aantal profielen

Slibvolume (m³)

12 11 10 9 8 7 6 5

97 95.7 94.3 93.5 97.8 99.1 91.1 86.7

Verschil t.o.v. max. aantal profielen (m³) 0 1.3 2.7 3.5 0.8 2.1 5.9 10.3

Verschil t.o.v. max. aantal profielen (%) 0 1.3 2.8 3.6 0.8 2.2 6.1 10.6

9

12

PRocentuele fout (%)

135° 12 11 10 9 8 7 6 5 4 3 2 1 0 5

6

7

8

10

11

Aantal profielen

Figuur 5-9 Visuele voorstelling foutenpercentage van een bocht van 135°

51

Figuur 5-10 Verschilkaart bocht van 135°, a (12 profielen), b (11 profielen), len), c (10 profielen), d (9 profielen), e (8 profielen), f (7 profielen), g (6 profielen), h (5 profielen)

Figuur 5-10 geeft de evolutie tie weer van de sedimentligging tijdens de variatie in het aantal profielen. Ook hier zien we het zelfde fenomeen als bij een bocht van 120°. In principe kunnen voor een bocht van 135° 7 profielen volstaan voor het inmeten. Het spreekt voor zich dat indien indie de omstandigheden qua arbeidsintensiteit het toelaten, 8 profielen een nog meer accuraat resultaat zullen opleveren. Hier wordt een foutenpercentage van 0.8% genoteerd.

52

5.1.6

Besluit: procedure inmetingen

In Figuur 5-11 wordt een overzicht gegeven van het aantal profielen dat wordt aanbevolen in functie van de aard van de bocht. Na het bekijken van de verschillende bochten kwamen drie hoofdtypes naar voor. Deze types maken het tijdens het veldwerk veel gemakkelijker om de bochten te beoordelen en accuraat in te meten. Dit is als het ware een soort van gestandaardiseerde procedure voor het praktisch inmeten. Figuur 5-11 Indeling in klassen

Klasse

Karakterisering

Aantal profielen

1

135-100°

8

2

80°

5

3

30-50°

4

Grafische weergave

53

5.2

Smoothing

5.2.1

Volume vergelijken

In Tabel 5-6 worden de resultaten van de smoothing weergegeven. In Figuur 5-12 wordt het foutenpercentage visueel voorgesteld. Tabel 5-6 Resultaten smoothing

Rastergrootte

Slibvolume (m³)

TIN 0.25 x 0.25 0.5 x 0.5 0.75 x 0.75 0.8 x 0.8 0.85 x 0.85 0.9 x 0.9 1x1 1.25 x 1.25 1.5 x 1.5

1062.4 1060.3 1054.4 1043.5 1042.4 1038.6 1037.8 1030.1 1011.7 992.9

Verschil t.o.v. max. aantal profielen (m³) 0 2.1 8 18.9 20 23.8 24.6 32.3 50.7 69.5

Verschil t.o.v. max. aantal profielen (%) 0 0.2 0.7 1.8 1.9 2.2 2.3 3.0 4.8 6.5

Smoothing 7 6

Procentuele fout (%)

5 4 3 2 1 0 0

0,2

0,4

0,6

0,8

1

1,2

1,4

1,6

Grootte raster

Figuur 5-12 Visuele voorstelling foutenpercentage smoothing

5.2.2

Visuele beoordeling

In Figuur 5-13 wordt de visuele voorstelling van een smoothing met raster 0.9 x 0.9 m weergegeven. De oorspronkelijke vaste bodem wordt weergegeven in donkerblauw, de gesmoothe vaste bodem wordt weergegeven in roos. Het oorspronkelijk waargenomen beeld van de bovenste laag van het slib (topslib) wordt in rood weergegeven, deze laag in gesmoothe toestand kan herkend worden aan de groene kleur. Tenslotte wordt in licht blauw de waterlijn gevisualiseerd.

54

Figuur 5-13 Visuele voorstelling smoothing 0.9 x 0.9 (m)

In bijlage 2 zijn de resultaten weergegeven van de visuele voorstelling van de verschillende rasters toegepast bij de smoothing evaluatie.

5.2.3

Besluit smoothing

Op basis van een visuele beoordeling van het smoothing resultaat en de procentuele fout die werd bevonden, blijktt het raster van 0.9 x 0.9 (m) het meest efficiënt. De procentuele fout is bij dit raster klein (2.3%) %) en men bekomt een mooie smoothing zonder uitstulpingen.

55

6 Besluit Het belangrijkste doel van dit onderzoek was het verwerven van wetenschappelijke inzichten in de verwerkingswijze en de technieken die gebruikt worden door het software pakket Terramodel ® van Trimble. Via literatuur werd de werking van het programma, in de context van de huidige verwerkingsprocedure, onderzocht. Door het opbouwen van deze kennis van de software werd het mogelijk om bij de inmeting op het veld veel efficiënter te werk te gaan. Het bepalen en het vastleggen van dwarsprofielen, nodig voor het verzamelen van sedimentgegevens, kon op deze manier effectiever en efficiënter gebeuren. Eveneens werden met deze kennis en kunde overbodige profielen vermeden. Op deze manier werd heel wat tijd gespaard zonder verlies aan informatie. Verder bestond deze studie uit het onderzoek naar de efficiëntste methode om verschillende types van meandering (bochten) van een typische Vlaamse laaglandbeek te karakteriseren. Via Terramodel werd per type bocht onderzocht hoeveel dwarsraaien nodig zijn om een accuraat beeld te krijgen maar tegelijkertijd ook efficiënt te werken. Door het volume aan slib bij een maximaal aantal dwarsraaien te vergelijken met een steeds afnemend aantal dwarsraaien werd per bocht een ideaal scenario opgesteld. Na de verwerking van de gegevens werd besloten dat bochten van Vlaamse waterlopen kunnen geclassificeerd en gekarakteriseerd worden in drie types, nodig voor het inpeilen van waterlopen/waterbodems. Deze types hebben een ideaal aantal dwarsraaien en maken het de surveyors van de Vlaamse Milieumaatschappij makkelijker om accuraat in te meten in een kortere tijdspanne, zonder verlies aan informatie. Op basis van deze informatie werd een ontwerp gestandaardiseerde procedure voor het praktisch inmeten van drie types bochten opgemaakt. Bij het opvolgen van ruimingsprojecten worden de verzamelde gegevens verwerkt tot een digitaal terrein model wordt bekomen (DTM). Aan de hand van dit model wordt ingeschat hoeveel volume sediment dient geruimd te worden. Deze informatie wordt aan de hand van profielen voorgelegd aan de aannemer die instaat voor de ruimingswerken. Na verloop van tijd ondervond de Vlaamse Milieumaatschappij dat aannemers problemen hadden met het lezen en toepassen van deze profielen. Aangezien ruimingswerken, in de meeste gevallen, uitgevoerd worden met graafmachines zijn meer gestroomlijnde, uitgemiddelde profielen gewenst. Gelet op de voorgaande kennis van Terramodel werd het mogelijk gemaakt om een ‘smoothing’ uit te voeren en dus zonder verlies aan informatie aannemers een bruikbaar en op het terrein toepasbaar DTM af te geven.

56

Bibliografie Artikels Dalyot, S., Doytsher, Y., & Katzil, Y. (2009). Digital Terrain Models: A Tool for Establishing Reliable and Qualitative Environmental Control Processes. NATO Science for Peace and Security Series C: Environmental Security , 215-234. Daniels, W., House, W., Rae, J., & Parker, A. (2000). The distribution of micro-organic contaminants in river bed- sediment cores. The Science of The Total Environment , 81-91. De Wulf, A., Hennau, M., & Constales, D. (2006). TIN-Based Digital Terrain Modelling Using Multibeam Data. Hydro , 186-189. Droj, G. (2008). Improving the accuracy of Digital Terrain Model. Informatica , 65-72. Keim, R., Skaugset, A., & Bateman, D. (1999). Digital terrain modeling of small stream channels with a total-station theodolite. Advances in Water Resources , 41-48. Lischinski, D. (1993). Incremental Delaunay Triangulation. New York: Academic Press. Su, P., & Drysdale, R. L. (1996). A Comparison of Sequential Delaunay Triangulation Algorithms. Vignati, D., Pardos, M., & Diserens, J. (2003). Characterisation of bed sediments and suspension of the river Po (Italy) during normal and high flow conditions. Water Research , 2847-2864.

Boeken Kluwer. (1998). Bodem en Bodemkunde.

Rapporten Rijkswaterstaat directie Limburg. (2003). Testplan Bodemligging Zomerbed, Inwinnen. STOWA. (2006). Inventarisatiemethodes voor het bepalen van baggervolumes. Utrecht: STOWA. Van Boxstael, B. (2009). Opmeten van de sedimentdikte m.b.v. Single Beam - Opstellen en implementatie van een gestandaardiseerde methodiek. Gent: Hogeschool Gent. Vito. (2007). Inzet van bagger-en ruimingsspecie ter vervanging van primaire grondstoffen in Vlaanderen.

Varia Belgisch Staatsblad. (2010, oktober 14). DECREET van 18 juli 2003 betreffende het integraal waterbeleid. Magelas. (2006). Bathymetrie, overzicht dieptemetingen. VMM. (2003). De verrekijker , 3.

Websites AGIV. (2011). Laseraltimetrie. Opgeroepen op April 1, 2011, van Agentschap voor Geografische Informatie Vlaanderen: http://www.agiv.be/gis/knowhow/?artid=548 Caverly, K. (sd). Spatial Interpolation, Overview. Opgeroepen op november 12, 2011, van btcsure1.ferris.edu/NGA/Interpolation_Krag.ppt Dadson, S. (1999). Digital Terrain Modelling. Opgeroepen op november 13, 2010, van http://www.geog.ubc.ca/courses/geog570/notes/dtm.html Departement Leefmilieu, natuur en energie van de Vlaamse overheid. (2010). Over onze organisatie - Organogrammen. Retrieved Oktober 23, 2010, from http://www.lne.be/organisatie/organogrammen EGIS foundation. (1994). Design and evaluation of TIN interpolation algorithms. Retrieved november 18, 2010, from http://libraries.maine.edu/Spatial/gisweb/spatdb/egis/eg94038.html EMS-i. (2011). Linear interpolation. Retrieved april 1, 2011, from Environmental Modeling Systems, Inc.: http://www.emsi.com/gmshelp/Interpolation/Interpolation_Schemes/Linear_Interpolation.htm Fürst, J. (sd). Interpolation of Hydrological Variables. Opgeroepen op december 2, 2010, van http://iwhw.boku.ac.at/gisww/bokudownloads/9_interpolation_en.ppt GeoLas Consulting. (2010). Imaging Laser Altimetry. Opgeroepen op februari 28, 2010, van http://www.geolas.com/Pages/laser.html GSE Rentals. (2010). Product information. Opgeroepen op december 28, 2010, van http://www.gserentals.co.uk/viewproduct2.cfm?prod_ref=150 Hemphill, J. J. (2009). TIN interpolation. Opgeroepen op december 22, 2010, van http://www.geog.ucsb.edu/~jeff/projects/la_conchita/tin_interpolation.html Kennis Centrum Waterbodem. (2006). Handmatig peilen. Opgeroepen op december 28, 2010, van http://www.kcwaterbodem.nl/page.php?144 Klinkenberg, B. (1997). The TIN model. Opgeroepen op november 26, 2010, van http://www.geog.ubc.ca/courses/klink/gis.notes/ncgia/u39.html Leica. (2010). Leica Photogrammetry Suite Automatic Terrain Exctraction. Opgeroepen op november 27, 2010, van http://gis.ess.washington.edu/data/erdas_pdfs/LPS_ATE.pdf McEachron, S. (2009). Paving the way. Opgeroepen op november 29, 2010, van http://blog.121pcs.com/terrain-modeling-delaunay-triangulation.html Metaldec Survey BV. (2010). Sidescan sonar survey. Opgeroepen op december 28, 2010, van http://www.metaldec.nl/sidescan.html Prospección Geofísica y Geofísica Pura. (2010, November 12). Georadar. Opgeroepen op februari 18, 2011, van http://www.ingenieriageofisica.com/georadar Trimble. (2010). Terramodel software for survey. Opgeroepen op november 26, 2010, van http://www.agrometius.nl/uploads/leaflets/TerraModelsoftware.pdf

VMM. (2010). Over VMM. Opgeroepen op oktober 23, 2010, van http://www.vmm.be/overvmm VMM. (2010). Situering VMM - organogram. Opgeroepen op oktober 23, 2010, van www.VMM.be/over-vmm/situering-vmm/organogram Wadembere, M. (2000). Demographics Interpolation and Extrapolation. Opgeroepen op november 17, 2010, van http://www.hbp.usm.my/thesis/heritageGIS/research/interpolation.htm

Bijlagen

Inhoudstafel BIJLAGE 1

Procedure voor het onderzoek naar Vlaamse bochten

BIJLAGE 2

Resultaten visuele voorstelling verschillende rasters voor smoothing evaluatie

Bijlage 1: procedure naar het onderzoek van Vlaamse bochten 40° Maximaal aantal profielen

Figuur 14 Vijf profielen

4 profielen

Figuur 15 vier profielen

3 profielen

Figuur 16 drie profielen

80° Maximaal aantal profielen

Figuur 17 Maximaal aantal profielen

7 profielen

Figuur 18 zeven profielen

6 Profielen

Figuur 19 zes profielen

5 profielen

Figuur 20 vijf profielen

4 profielen

Figuur 21 vier profielen

105° Maximaal aantal profielen

Figuur 22 Maximaal aantal profielen

9 profielen

Figuur 23 negen profielen

8 profielen

Figuur 24 acht profielen

7 profielen

Figuur 25 zeven profielen

6 profielen

Figuur 26 zes profielen

120° Alle profielen

Figuur 27 Maximaal aantal profielen

8 profielen

Figuur 28 acht profielen

7 profielen

Figuur 29 Zeven profielen

6 profielen

Figuur 30 Zes profielen

5 profielen

Figuur 31 Vijf profielen

4 profielen

Figuur 32 Vier profielen

135° 12 profielen

Figuur 33 Maximaal aantal profielen

11 profielen

Figuur 34 elf profielen

10 profielen

Figuur 35 Tien profielen

9 profielen

Figuur 36 Negen profielen

8 profielen

Figuur 37 Acht profielen

7 profielen

Figuur 38 Zeven profielen

6 profielen

Figuur 39 Zes profielen

5 profielen

Figuur 40 Vijf profielen

Bijlage 2:: Resultaten visuele voorstelling verschillende rasters voor smoothing evaluatie 0.25 x 0.25 m

0.5 x 0.5 m

0.75 x 0.75 m

0.8 x 0.8 m

0.85 x 0.85 m

0.9 x 0.9 m

1x1m

1.25 x 1.25 m

1.5 x 1.5 m