15.07.2013 Views

Rapporten - Søren Højmark Rasmussen

Rapporten - Søren Højmark Rasmussen

Rapporten - Søren Højmark Rasmussen

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Vandindvinding i Volsted Plantage<br />

– Numerisk modellering af vandføringen i Hølbækken<br />

Naturgeografi, 8. semester, Forår 2006


Titel: Vandindvinding i Volsted Plantage<br />

B-studienævnet<br />

Aalborg Universitet<br />

Sohngaardsholmsvej 57<br />

9000 Aalborg<br />

www.bsn.aau.dk<br />

- Numerisk modellering af vandføringen i Hølbækken<br />

Tema: Grundvandet som ressource<br />

Projektperiode: 1/2 - 2/6 2006<br />

Projektgruppe: D019 Synopsis<br />

Deltagere:<br />

Bastian E. Jensen<br />

Marianne Bismo<br />

Martin Thorsøe<br />

<strong>Søren</strong> H. <strong>Rasmussen</strong><br />

Vejledere:<br />

Jacob Birk Jensen<br />

Niels Arne Wahl<br />

Oplagstal: 7<br />

Sideantal: 105<br />

Bilagsantal: 7 og en Cd-rom<br />

Aalborg kommune planlægger at oprette<br />

en vandindvinding i Volsted Plantage,<br />

hvor der skal indvindes ca. 1 mio. m 3 /år.<br />

Formålet med denne rapport er at<br />

undersøge, hvilken indvirkning dette vil<br />

have på vandføringen i det nærliggende<br />

vandløb Hølbækken, der er klassificeret<br />

som B1 (laksefiskevand) og derfor kun<br />

må påvirkes minimalt. Dette er<br />

undersøgt ved at opsætte en konceptuel<br />

model for området omkring Hølbækken,<br />

der så er brugt til opstillingen af en<br />

numerisk grundvandsmodel. Beregninger<br />

udført med grundvandsmodellen<br />

viser, at en indvinding på 1 mio. m 3 /år<br />

giver en reduktion af Hølbækkens gennemsnitsvandføring<br />

på ca. 20 % og at<br />

denne påvirkes 10 %, hvis der i stedet<br />

indvindes 1/2 mio. m 3 /år. Det konkluderes<br />

derfor, at den planlagte indvinding<br />

påvirker Hølbækkens vand-føring med<br />

mere end det tilladte.


English summary<br />

The municipality of Aalborg plans a groundwater extraction from Volsted Plantage. The<br />

purpose is to supply the citizens of the municipality with 1.000.000 m 3 drinking water<br />

pr. yr. The object of this paper is to estimate the influence this extraction will have on<br />

Hølbækken, a nearby brook. Hølbækken is classified as a salmon habitat. This means<br />

that the median of the annual minimum flow in the brook can only be affected<br />

minimally.<br />

The project area was interpreted by making a conceptual model, in which the geology,<br />

hydrology and hydrogeology of the project area was investigated. On the basis of the<br />

conceptual model a numerical groundwater model was developed. The groundwater<br />

model was generated using GMS 5.1. Two different geological interpretations of the<br />

model area and two different placements of the extraction wells in Volsted Plantage<br />

were simulated using the model. The different calculations and placements produced<br />

similar results. The calculations showed a 20 % decline of the water flow in Hølbækken<br />

when 1.000.000 m 3 /yr was extracted. The water flow declined by 10 % when 500.000<br />

m 3 /yr was extracted. This paper concludes that the planned extraction of 1.000.000<br />

m 3 /yr influences the water flow in Hølbækken by more than the accepted according to<br />

the salmon habitat classification.


Forord<br />

Denne rapport er udarbejdet på 8. semester af Cand. Scient. uddannelsen Geografi med<br />

speciale i Naturgeografi. Temaet for dette semester er Grundvandet som ressource. Det<br />

er valgt at arbejdet med dette tema ved at undersøge den mulige påvirkning ved<br />

etableringen af en ny drikkevandsindvinding i Volsted Plantage. Målet med rapporten er<br />

derfor, at lave en ressourcemæssig vurdering af grundvandsressourcen i oplandet til<br />

plantagen på baggrund af en grundvandsmodel. Til opsætningen af grundvandsmodellen<br />

er der blevet indsamlet data på to feltture den 27.02 - 03.03 og 20.04 - 21.04.<br />

Kildehenvisninger er angivet på formen (forfatters efternavn, årstal: sidetal). Tabeller og<br />

figurer er nummereret fortløbende. På den vedlagte Cd-rom findes de beregninger, som<br />

ligger til grund for resultaterne fra de undersøgelser, der præsenteres i rapporten.<br />

Bastian E. Jensen Marianne Bismo<br />

Martin H. Thorsøe <strong>Søren</strong> H. <strong>Rasmussen</strong>


Indhold<br />

Side<br />

1 Indledning 7<br />

1.1 Problemformulering 10<br />

2 Grundvand- og overfladevandsinteraktion 15<br />

2.1 Ådalen 16<br />

2.2 Ådales geomorfologi 17<br />

2.3 Vandbalanceligningen 19<br />

2.4 Opsummering 23<br />

3 Konceptuel model af projektområdet 25<br />

3.1 Præsentation af fokusområdet 26<br />

3.2 Projektområdets geologi 30<br />

3.3 Projektområdets hydrologi 46<br />

3.4 Projektområdets hydrogeologi 56<br />

3.5 Opsummering 66<br />

4 Opsætning af grundvandsmodellen 69<br />

4.1 Grundvandsmodellering generelt 69<br />

4.2 Grundvandsmodellen 71<br />

4.3 Kalibrering af grundvandsmodellen 76<br />

4.4 Placering af indvindingsboringerne 81<br />

4.5 Anvendelse af grundvandsmodellen 84<br />

4.6 Effekt af grundvandsindvindingen på vandføringen i Hølbækken 85<br />

4.7 Usikkerheder ved grundvandsmodellering 87<br />

5 Konklusion 97<br />

6 Perspektivering 99<br />

Litteraturliste 101<br />

Bilag<br />

A Flowlog<br />

B Vandføring<br />

C Fasebestemmelse<br />

D Hydraulisk ledningsevne<br />

E Slugtest<br />

F Åbundens hydrauliske ledningsevne<br />

G Diffus udstrømning


1<br />

Indledning<br />

Vand er en ressource, der kan anvendes til en lang række forskellige formål. Industrien<br />

har ofte brug for store mængder vand til køling eller procesvand. Landbruget har brug<br />

for vanding af deres marker. Vi mennesker har brug rent drikkevand. Dette vand kan<br />

hentes fra forskellige kilder afhængig af hvad, der er til rådighed, men i Danmark er den<br />

største kilde grundvandet, hvorfra 99 % af vandet indvindes (Danva, 2006). Til industri<br />

og landbrug anvendes stadigvæk i nogen udstrækning vand fra f.eks. vandløb, men<br />

drikkevand indvindes kun fra grundvandet. Rent drikkevand er en livsnødvendighed,<br />

hvilket betyder, at alene af denne grund, er beskyttelsen af grundvandet og sikring af<br />

adgangen til tilstrækkelige mængder en høj prioritet.<br />

Grundvandet kan forurenes af forskellige kilder og dermed blive uegnet til anvendelse.<br />

Forureningen kan f.eks. komme fra industri eller udslip af kemikalier fra utætte tanke.<br />

En væsentlig kilde til forurening af grundvandet er landbruget, hvis anvendelse af<br />

gødning og pesticider kan lede til påvirkning af store arealer. Landbrugets udnyttelsesgrad<br />

af gødning, som ikke er fuldstændig, betyder, at der sker en nedsivning af nitrat og<br />

pesticider til grundvandet. Samtidig er mange af de pesticidtyper, som anvendes,<br />

langsomt nedbrydelige og derfor kan der ske en ophobning af pesticider i grundvandet,<br />

som kan tage mange år om at nedbrydes. Den politiske målsætning i Danmark er, at<br />

uforurenet grundvand skal være basis for drikkevand og at vandrensning kun<br />

undtagelsesvist og i en begrænset periode kan tillades (Regionplan, 2005: 238).<br />

En række problemer kan opstå i forbindelse med selve indvindingen. Når der pumpes<br />

vand fra en boring, sænkes grundvandspotentialet i det omkring denne. Indvindingen<br />

kan også føre til saltindtrængning i kystnære områder eller tørlægning af okkerholdige


Kapitel 1 Indledning<br />

lag i undergrunden. Begge disse problemer betyder, at vandet fra denne boring ikke<br />

længere kan anvendes til drikkevand. Sænkning af grundvandspotentialet kan også<br />

påvirke nærliggende hydrologiske systemer, såsom vådområder eller åer, hvor det kan<br />

føre til sænkning af vandstanden og ændring af de økologiske forhold. Dette er i strid<br />

med den gældende lovgivning, hvor påvirkningen af lokale økosystemer helst skal<br />

undgås eller minimeres. I praksis er der to EU direktiver som omhandler dette problem:<br />

Vandrammedirektivet og Habitatdirektivet. Begge disse direktiver er implementeret i<br />

dansk lovgivning. Habitatdirektivet omhandler beskyttelse af områder med et truet eller<br />

sjældent dyre- eller planteliv og områder med truede, arealmæssigt begrænsede eller<br />

karakteristiske naturtyper. Vandrammedirektivet er en sammentænkning af det hydrologiske<br />

kredsløb, hvor det gælder om at sikre den bedst mulige kvalitet af både<br />

grundvand og overfladevand. Kvaliteten går både på kemiske parametre og økologiske<br />

parametre. Sidstnævnte indebærer bl.a. en tilbageføring af vandløbene til deres<br />

”naturlige” tilstand, dvs. den menneskelige indflydelse skal fjernes. Begge direktiver<br />

indebærer, at der er opstillet en række målsætninger for samtlige danske vandløb, som<br />

skal overholdes (SNS, 2006).<br />

I Nordjyllands Amt har der gennem de seneste fem år været en stigning i antallet af<br />

lukkede boringer, samtidig med, at der ikke er sket et fald i vandforbruget, figur 1.1.<br />

Derfor er der behov for at nye boringer oprettes for at imødekomme behovet for vand.<br />

Derudover er det nødvendigt for vandværkerne at operere med en vis overkapacitet i<br />

deres system, hvilket vil sige, at de har flere boringer med en større kapacitet end det er<br />

nødvendigt for at opfylde behovet for rent drikkevand. Derved kan de stadigvæk levere<br />

den nødvendige mængde vand, selvom en af deres boringer må lukke. Det betyder også,<br />

at de kan sprede deres vandindvinding ud på flere boringer og mindske påvirkningen på<br />

enkelte lokaliteter.<br />

<strong>Rapporten</strong>s projektområde<br />

For at sikre vandforsyningen har Aalborg Kommune i øjeblikket gang i en undersøgelse<br />

af Volsted Plantage, ca. 20 km syd for Aalborg centrum, med henblik på etablering af<br />

en ny kildeplads, figur 1.2. Der er planer om at indvinde 1 mio. m 3 om året fra denne<br />

kildeplads, hvis det viser sig, at området er egnet og de miljømæssige påvirkninger ikke<br />

er for store. Den potentielle kildeplads er beskyttet af skoven og undersøgelser viser, at<br />

det er muligt at indvinde rent drikkevand fra denne lokalitet. Det er imidlertid ikke klart,<br />

hvilken indvirkning en boring vil have på det omkringliggende område.<br />

8


Antal lukkede boringer<br />

100<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

20<br />

10<br />

0<br />

2000 2001 2002 2003 2004 2005<br />

År<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

20<br />

10<br />

0<br />

Forbrug i mio. m 3<br />

Ikke angivet<br />

Andet<br />

Forurening<br />

Forbrugt<br />

drikkevand<br />

Figur 1.1. Forbrug af vand i Nordjyllands Amt og antallet af lukkede boringer pr. år i perioden år 2000 –<br />

2005, fordelt på årsag. Kategorien ”andet” dækker bl.a. over boringer, der er blevet lukkede af<br />

administrative årsager, fordi de ikke længere anvendes, eller fordi de ikke længere fungerer. Data fra<br />

GEUS 1 (2006) og Statistikbanken (2006).<br />

Udover kommunen er der andre interesser i området, som kunne blive påvirket af<br />

vandindvindingen. Der er således et dambrug, lodsejere og beskyttede naturområder.<br />

Dambruget indvinder grundvand fra kilder, som benyttes til produktionen af<br />

opdrætsfisk. Hvis kommunen også begynder at indvinde vand i området, vil det kunne<br />

påvirke mængden og kvaliteten af det vand, som dambruget kan indvinde. Lodsejerne<br />

har gennem deres ejerskab en interesse i området, placeringen af indvindings boringer i<br />

Volsted Plantage vil lægge restriktioner på arealanvendelsen. Endelig er Lindenborg<br />

Ådal et EU Habitatområde, hvor vandføringen i åen ikke må påvirkes, dette område<br />

omfatter også den nedre del af Hølbækken. Hølbækken er desuden målsat som B1,<br />

yngelopvækstvand for laksefisk (Bennedsen, 2002: 81). B1 målsætningen betyder, at<br />

vandløbets vandføring kun må påvirkes med 5-10 % af medianminimumsvandføringen<br />

(Regionplan, 2005: 237). Medianminimumsvandføringen er median af hvert års min.<br />

vandføring normalt over tyve år (DMU 1, 2006).<br />

9


Kapitel 1 Indledning<br />

Figur 1.2. Kortet viser placeringen af Volsted Plantage syd for Aalborg. Data fra TOP10DK (2000).<br />

Signaturen fra dette kort vil gå igen på flere af rapportens kort.<br />

1.1 Problemformulering<br />

I denne rapport er det valgt at fokusere på grundvandindvindingens påvirkning af<br />

Hølbækken. Da Hølbækken ligger tæt på den potentielle indvindingsplads, er der en stor<br />

mulighed for at vandføringen i bækken kan blive påvirket af indvindingen. Samtidig er<br />

undergrunden ved Hølbækken forholdsvist ukendt, så det er vanskeligt at vurdere<br />

præcist, hvilken effekt vandindvinding i Voldsted Plantage kan få. Derfor kunne det<br />

være interessant at undersøge dette område nærmere, med henblik på at belyse effekten<br />

af vandindvinding i forhold til Hølbækkens B1 målsætning. Dette leder frem til<br />

følgende problemstilling for denne rapport:<br />

10<br />

Hvilken indflydelse vil oprettelsen af indvindingsboringer i<br />

Volsted Plantage få for Hølbækkens vandføring?<br />

Spørgsmålet ønskes undersøgt kvantitativt, idet der således kan undersøges, hvorvidt<br />

grundvandsindvindingens indflydelse er af en sådan grad, at Hølbækken ikke kan leve<br />

op til B1 målsætningen. Der vil blive opsat en numerisk grundvandsmodel for området.<br />

Ifølge Refsgaard (2005), består en grundvandsmodellering af fem faser der skal


gennemarbejdes, figur 1.3. Den videre rapport vil blive bygget op med baggrund i<br />

denne opdeling af modelleringsarbejdet.<br />

Modelstudieplan<br />

Konceptualiseringen og<br />

dataindsamling<br />

Modelopsætning<br />

Kalibrering og validering<br />

Simulering og<br />

evaluering<br />

Figur 1.3. Arbejdsgangen i forbindelse med opsætningen af en grundvandsmodel (Refsgaard, 2005).<br />

Afgrænsning af modelområdet<br />

Første skridt i arbejdet med en grundvandsmodel er at lave en modelstudieplan, der er<br />

en udarbejdelse af formålet med modelleringen, en analyse af opgavens omfang og en<br />

plan for udarbejdelsen af modellen. Formålet med modelleringen er præsenteret i det<br />

foregående. For at danne et overblik over modelleringsopgaven, er modellens afgrænsning<br />

blevet undersøgt og vil i det følgende blive præsenteret.<br />

Ved afgrænsning af oplandet for en grundsvandsmodel er det vigtigt, at der er styr på<br />

vandbalancen i området, da der skal opsættes en numerisk grundvandsmodel over<br />

området. Der er anvendt to metoder til at sikre dette. Den første er nul fluks grænser,<br />

vinkelret på potentialekurverne. Da der dermed er det samme tryk på begge sider af<br />

afgrænsningen, vil der ikke ske nogen vandtransport over denne grænse. Den anden<br />

metode er trykgrænser, hvor potentialet langs grænsen er kendt. Den udstrømning, der<br />

sker fra disse områder er kendt, idet trykniveaut er kendt. Projektområdet er blevet<br />

afgrænset ved brug af et potentialekort, hvor vandløb og dræn er indtegnet, figur 1.4.<br />

Mod syd er projektområdet afgrænset af Lindenborg Å og mod vest af åen Skriveren,<br />

11


Kapitel 1 Indledning<br />

der er en tilløbså til Lindenborg Å. Disse to åer fungerer som trykgrænser. Den øvrige<br />

del af området er afgrænset med nul fluks grænser. Denne modelafgrænsning vil blive<br />

anvendt i resten af rapporten og der vil blive henvist til den som området. Hølbækken<br />

og dens ådal defineres som fokusområde for rapporten, hvilket skyldes, at det er i dette<br />

område påvirkningen ønskes bestemt og der vil derfor blive gået i dybden med<br />

undersøgelser i netop dette område.<br />

Figur 1.4. Projektområdets afgrænsning og fokusområdet. Tema for grundvandspotentiale er interpoleret<br />

ud fra punktobservationer af potentialet fra boringer (NJA, 2005). Boredata er af varierende kvalitet og<br />

alder, hvorfor temaet kun er vejledende. Å-tema fra TOP10DK (2000).<br />

Plan for modelleringen<br />

Den resterende rapport har følgende opbygning:<br />

12<br />

Grundvand- overfladevandsinteraktion – I dette kapitel præsenteres<br />

baggrundsviden, der bliver brugt i rapporten i forbindelse med grundvand-<br />

overfladevandsinteraktion. Der lægges vægt på at fastlægge relevante termer og<br />

principper for de undersøgelser, der skal belyse projektområdet.<br />

Konceptuel model af projektområdet – Her vil konceptualiseringen og<br />

dataindsamlingen til grundvandsmodellen foregå. Dette vil foregå i fire dele,<br />

først gennem en beskrivelse af området med en kortanalyse, samt observationer.


Dernæst gennem opstillingen af en geologisk model for området, der opbygges<br />

på baggrund af geologisk viden om projektområdet, GEUS boringer samt<br />

gruppens undersøgelser af undergrunden i Hølbækkens ådal. Udover dette vil<br />

kapitlet indeholde en gennemgang af de hydrologiske processer, herunder<br />

vandføring, infiltration og de forskellige geologiske lags hydrauliske ledningsevne.<br />

Grundvandsmodel – Valg af rammer for opsætningen af grundvandsmodellen<br />

fremlægges og modelopsætningen beskrives. Derudover vil modellen i kapitlet<br />

blive kalibreret og valideret. Herefter præsenteres og diskuteres modellens<br />

resultater. Kapitlet afsluttes med en diskussion af usikkerheder ved modellen.<br />

13


2<br />

Grundvand- og<br />

overfladevandsinteraktion<br />

I det følgende kapitel vil der blive set på den interaktion, der er mellem grundvand og<br />

overfladevandet i en ådal. Det vil blive gennemgået, hvilken sammenhæng der er<br />

mellem et områdes geomorfologi og den udveksling af vand, der foregår mellem<br />

vandløbet og oplandet. Det er vigtigt at have et kendskab til de faktorer, der påvirker<br />

grundvandets interaktion med overfladevandet, da det er formålet med denne rapport at<br />

klarlægge, hvilken påvirkning en indvinding af grundvand fra Volsted Plantage vil have<br />

på Hølbækken. Med baggrund i den teori som bliver præsenteret i dette kapitel, er det<br />

muligt at bestemme hvilke parametre, der skal undersøges nærmere for at kunne opstille<br />

en grundvandsmodel. Kapitlet skal også give en forståelse af de processer, der foregår i<br />

ådalene, som er nødvendig for at kunne analysere de resultater, der kommer frem i<br />

undersøgelserne.<br />

Grundvands- og overfladevandsinteraktion (GOI) er et komplekst område og forbindes<br />

ofte med en mængde forskellige naturvidenskabelige fagretninger, blandt andet geologi,<br />

biologi, naturgeografi og hydrologi. Grundlaget for gennemgangen er den GOItypologi,<br />

DMU har udviklet (Nilsson et al., 2005). Der vil udelukkende blive lagt vægt<br />

på de elementer, der er relevante for denne.


Kapitel 2 Grundvand- og overfladevandsinteraktion<br />

2.1 Ådalen<br />

I Danmark er ådale udstrømningsområder for grundvand. Der kan både være tale om<br />

overfladenær udstrømning og udstrømning fra dybere lag (Nilsson et al., 2004: 22).<br />

Interaktionen mellem grundvand og overfladevand kan derfor betragtes som den<br />

udveksling, der sker mellem vandløbet, ådalsmagasinet og den tilgrænsende hydrologiske<br />

enhed, figur 2.1 (Nilsson et al., 2003: 19). En ådal dannes, når der sker en<br />

nedskæring af ældre geologiske aflejringer og hermed dannes en dal som delvis fyldes<br />

op af nyere aflejringer. I Danmark er dannelsen ofte forbundet med gletscheres<br />

virksomhed. I ådalen findes der et ådalsmagasin bestående af yngre geologiske<br />

aflejringer, som langs siderne afgrænses af en skræntfod. Det er forskelligt, hvor let det<br />

er at skelne ådalsmagasinet fra den tilgrænsende hydrologiske enhed. I nogle tilfælde<br />

består ådalsmagasinet f.eks. af den samme type aflejringer som den tilgrænsende<br />

hydrologiske enhed, hvilket vanskeliggør en distinktion. Hvordan vandudvekslingen<br />

mellem den tilgrænsende hydrologiske enhed og ådalen er, afgøres af størrelsen,<br />

topografien og formen af oplandet og den tilgrænsende hydrologiske enheds<br />

hydrauliske ledningsevne.<br />

Figur 2.1. Begreber der knytter sig til en ådal (Nilsson et al., 2004: 22).<br />

16<br />

Definitioner<br />

Geologi – Læren om Jordens fysiske struktur og Jordens historie<br />

om forskellige istider, og hvordan de ændrede landskabet.<br />

Hydrologi – Læren om vand og vands bevægelse på jordens<br />

overflade.<br />

Hydrogeologi – Den del af hydrologien, der beskæftiger sig med<br />

vandets strømning i de geologiske lag.<br />

(Wikipedia, 2006)<br />

Opland Ådal Opland<br />

Skrænt<br />

Skræntfod<br />

Vandløb<br />

Ådalsmagasin<br />

Tilgrænsende hydrologisk enhed<br />

Relief<br />

Opfyldning


2.2 Ådales geomorfologi<br />

I det følgende vil der blive set på ådalenes udformning og der vil blive set nærmere på,<br />

hvordan ådale kan inddeles i forskellige typer, afhængig af deres geomorfologi. Der er<br />

stor sammenhæng mellem landskabstypen og interaktionen mellem grundvand og<br />

overfladevand og at det derfor er centralt at analysere det landskab, ådalen findes i, hvis<br />

det er ønsket at sige noget om vandudvekslingen i ådalen (Langhoff et al., 2005: 15).<br />

Der tages udgangspunkt i nogle overordnede inddelinger og der vil efterfølgende blive<br />

gået i dybden med den klassificering, der anvendes til danske forhold i forbindelse med<br />

GOI. Der vil blive lagt vægt på at præsentere en geomorfologisk inddeling, der senere<br />

kan kobles med de hydrologiske processer, der finder sted i ådalene. Grunden til at det<br />

er interessant at kende et områdes geomorfologi er, at der er en stærk sammenhæng<br />

mellem geomorfologien og hydrologien i en ådal og det er ofte lettere at skaffe<br />

geomorfologiske data (Larkin & Sharp, 1992: 1608).<br />

Overordnet set kan vandløbsstrækninger karakteriseres som enten aggraderende, hvor<br />

der sker en sedimentation eller degraderende, hvor der sker en erosion (Dahm et al.,<br />

1998: 430). Om en vandløbsstrækning er aggraderende eller degraderende bestemmes<br />

af, hvor hurtigt vandet strømmer og hvilke typer af sedimenter, der findes i vandet. Det<br />

er derfor en inddeling, der kan ændres over tid, eftersom strømningshastigheden og<br />

mængden af suspenderet materiale varieres over tid.<br />

Vandløbet kan enten være begrænset eller ubegrænset. En begrænset vandløbsstrækning<br />

er kendetegnet ved, at ådalen er mindre end fire gange så bred som åen, mens en<br />

ubegrænset vandløbsstrækning er kendetegnet ved, at ådalen er bredere end fire gange<br />

vandløbets brede (Dahm et al., 1998: 430). En ubegrænset vandløbsstrækning kan både<br />

være aggraderende og degraderende, mens en begrænset vandløbsstrækning<br />

sædvanligvis kun er degraderende. På strækninger, hvor vandløbet er aggraderende, er<br />

interaktionen mellem grundvand og overfladevand særligt høj, hvilket også er tilfældet<br />

på strækninger, hvor begrænsede vandløb åbner op og bliver ubegrænsede. Derimod er<br />

strækninger, hvor vandløbet er degraderende kendetegnet ved, at have en mindre<br />

interaktion mellem grundvandet og overfladevandet.<br />

I Nilsson et al. (2004) opdeles de danske ådale i fem forskellige typer, afhængende af<br />

områdernes geomorfologiske og geologiske historie.<br />

1. Morænelandskab: Et morænelandskab består af både sandede og lerede<br />

aflejringer, og repræsenterer den landskabstype i Danmark, der har det største<br />

relief. Udviklingen af en ådal er ofte ikke så udpræget i denne type landskab og<br />

indholdet af organisk materiale langs åen ikke er så stort. Derfor vil<br />

17


Kapitel 2 Grundvand- og overfladevandsinteraktion<br />

18<br />

tilstrømningen til et vandløb i et morænelandskab ofte være lokal og vil variere<br />

meget gennem året.<br />

2. Smeltevandsslette: En smeltevandsslette består af relativt homogene sandede<br />

aflejringer med en svag hældning. Vandløb på en smeltevandsslette er ofte dybt<br />

nedskårne i det omgivne landskab og vandføringen varierer ofte meget i løbet af<br />

året, særligt på det øvre løb, hvor vandløbet ofte mister kontakten med den<br />

tilgrænsende hydrauliske enhed om sommeren.<br />

3. Havaflejring: Havaflejringerne stammer fra Yoldia- og Litorinahavet og er<br />

derfor ofte fintsandede eller lerede. Disse er særlig afsat i de områder, hvor<br />

gletschere fra tidligere tider har fjernet materiale. Ådalen i disse typer af<br />

aflejringer har ofte et udpræget relief og der findes her et højt indhold af<br />

organisk materiale.<br />

4. Fladedækkede ferskvandsaflejringer: Fladedækkende ferskvandsaflejringer<br />

har en udbredelse, der er over to km bred, består af tørv og kan være afsat oven<br />

på havaflejringer. Ådalen er derfor bred og ådalsmagasinet vil have en dybde af<br />

flere meter. Ådalen er karakteriseret ved at have et meget lavt relief og vil derfor<br />

ofte oversvømmes om vinteren.<br />

5. Klitlandskab: Klitlandskaberne består af fintsandede flyvesandsaflejringer og<br />

findes ofte langs kysten. Derfor er der i denne landskabstype ingen eller kun lille<br />

udvikling af en ådal og vandløb i denne landskabstype er ofte meget små med<br />

stor variation i vandføringen.<br />

En inddeling af et område efter den geologiske og geomorfologiske dannelseshistorie<br />

giver et overordnet billede af, hvordan de hydrologiske forhold er i området og<br />

endvidere et indtryk af, hvor stort et område resultaterne af en undersøgelse er<br />

anvendelige for. I forhold til interaktionen mellem grundvandet og overfladevandet er<br />

relieffet og de materialer, der findes i ådalen, væsentlige. Materialet er med til at afgøre<br />

den hydrauliske ledningsevne, mens relieffet er med til at afgøre potentialegradienten.<br />

Burt et al. (2002) har i et studie af ådale i 8 forskellige europæiske lande vist, at ådalens<br />

form også spiller en stor rolle for vandudvekslingen. De anvendte ådalene i studiet var<br />

lokaliseret på forskellige typer af smeltevandsaflejringer. I studiet kom det frem, at<br />

vandløb i en bred ådal med flade skrænter har en mere jævn vandføring end en smallere<br />

ådal med stejlere skrænter. Årsagen til dette er, at grundvandsspejlets overflade vil følge<br />

konturerne i landskabet og ådale med forskellig relief vil derfor også have forskellige<br />

grundvandsgradienter. En ådal med stejle skrænter har derfor også et grundvandsspejl


med en stejl gradient, der varierer meget gennem året, mens en ådal med fladere<br />

skrænter vil have en fladere grundvandsgradient, der variere mindre<br />

Et andet geomorfologisk karakteristika der kan anvendes til at karakterisere<br />

interaktionen mellem grundvand og overfladevand er åens slyngningsgrad. Slyngninger<br />

betyder, at vandet skal transporteres en længere strækning end hvis åen var lige, hvilket<br />

fører til, at vandløbsgradienten bliver lavere, jo flere slyngninger vandløbet har.<br />

Vandløbets slyngningsgrad bruges derfor som et mål for at karakterisere vandløbets<br />

interaktion med ådalsmagasinet. Slyngningsgraden er defineret som vandløbets længde,<br />

divideret med den kortest mulige afstand mellem to punkter på vandløbet (Langhoff et<br />

al., 2005: 5). Jo højere vandløbets slyngningsgrad er, jo større er interaktionen mellem<br />

grundvand og overfladevand.<br />

Det er endvidere værd at bemærke, at et vandløb ikke kun strømmer igennem en<br />

landskabstype, men at landskabstypen ofte skifter fra vandløbets udspring til udløbet i<br />

havet, og at de fem overordnet landskabstyper, kan underinddeles. Mennesket har også i<br />

stor grad påvirket ådalenes fysiske form, hvilket besværliggør en karakterisering af<br />

vandløbene, da det udseende de har i dag, ofte ikke afspejler deres naturlige form. Det<br />

kan derfor være nødvendigt at anvende historiske data over vandløbene, hvis det er<br />

ønsket sige noget om den naturlige situation i området.<br />

2.3 Vandbalanceligningen<br />

Bevægelsen af vand i undergrunden styres af potentialeforskelle, hvor vandet bevæger<br />

sig fra områder med et højt potentiale til områder med lavt potentiale. Vandtransport i<br />

undergrunden er en del af det hydrologiske kredsløb og for at have styr på vandtransporten,<br />

er det nødvendigt at have overblik over hele kredsløbet. Til dette formål er<br />

vandbalanceligningen anvendelig. Vandbalancen bruges til, at bestemme vandmængden<br />

der kommer ind og ud af et afgrænset område. I ligning 1.1. kan der ses en udgave af<br />

vandbalanceligningen, hvor tilførslen af vand, dvs. nedbøren er lig den mængde, der<br />

fjernes samt en ændring i grundvandreservoirets magasinering.<br />

Et vandbalanceregnskab kan også stilles op over en ådal og kan her anvendes til at<br />

undersøge, hvordan vandet bevæger sig i ådalen, figur 2.2. I forhold til denne rapports<br />

problemstilling er det nødvendigt, at der opstilles et vandbalanceregnskab over<br />

Hølbækkens ådal. Dette skyldes, at oppumpningen af vand i oplandet til vandløbet er et<br />

indgreb i Hølbækkens vandbalance, og formålet med denne rapport er at undersøge,<br />

hvor stor påvirkning denne oppumpning vil få for Hølbækkens vandføring, (Ao og Au).<br />

19


Kapitel 2 Grundvand- og overfladevandsinteraktion<br />

Ligning 1.1. Vandbalanceligningen (Burcharth & Jørgensen, 1976: 5).<br />

Vandudvekslingen mellem grundvandet og overfladevandet, kan i ådalen ske på fire<br />

forskellige måder, som hver især kan betragtes som et led i et vandbalanceregnskab<br />

(Nilsson et al., 2004). Det er værd at være opmærksom på, at en enkelt vandløbsstrækning<br />

ikke kun får tilført vand fra den ene af disse strømningsveje, men at der er et<br />

stærkt sammenspil mellem de forskellige strømningsveje, altså at vandløbet får tilført<br />

vand både via overfladisk, overfladenær, via dræn og dyb grundvandsafstrømning.<br />

Strømningsvejene er illustreret på figur 2.2. og forklares herunder.<br />

Figur 2.2. Vandudvekslingen mellem en ådal og den tilgrænsende hydrauliske enhed (Nilsson et al.,<br />

2003: 23).<br />

20<br />

Lokal<br />

overfladenær<br />

grundvandstilstrømning<br />

Nedbør Fordampning<br />

Overfladenær<br />

grundvandsafstrømning<br />

N = E + Au + Ao + Q + M<br />

Hvor:<br />

N : Nedbør Ao : Overfladisk afstrømning<br />

E : Fordampning Q : Oppumpet vandmængde<br />

Au : Underjordisk afstrømning M : Ændring i jordens magasinering<br />

Overfladisk afstrømning<br />

Dyb grundvands<br />

afstrømning<br />

Dyb regional grundvandstilstrømning<br />

Afstrømning fra dræn<br />

Dyb grundvandsafstrømning: Den dybe grundvandsafstrømning stammer fra<br />

regionalt dannet grundvand, som perkolerer op gennem vandløbsbunden.<br />

Vandet, der afstrømmer denne vej, har derfor kun begrænset kontakt med<br />

ådalsmagasinet. Dette drejer sig særligt om vandløb, hvor den hydrauliske<br />

ledningsevne i bundens sedimenter er høj. Den forekommer også, hvis vandløbet<br />

er blevet kraftigt reguleret, således at der er blevet gravet gennem ådalens<br />

sedimenter og der nu er direkte kontakt mellem vandløbet og den tilgrænsende


hydrologiske enhed. Vandets opholdstid i ådalen forventes derfor at være fra<br />

timer til dage.<br />

Overfladenær afstrømning: Den overfladenære afstrømning sker ved, at<br />

grundvandet i den tilgrænsende hydrologiske enhed strømmer gennem<br />

ådalsmagasinet og diffust strømmer ud i vandløbet. Er sedimenterne i ådalen<br />

højpermeable, vil grundvandet lettere strømme af som overfladenær<br />

afstrømning. Vandet, der strømmer af overfladenært kan være grundvand, der er<br />

dannet både lokalt og regionalt og vandets opholdstid i ådalsmagasinet er derfor<br />

mellem uger og år.<br />

Overfladisk afstrømning: Den overfladiske afstrømning sker oven på<br />

sedimenterne i ådalen. Dette sker hovedsagelig, når ådalens sedimenter er<br />

lavpermeable og ikke kan transportere så meget vand. Derudover sker det, hvis<br />

der findes kilder i ådalen, der tilfører vand til åen. Denne form for afstrømning<br />

forekommer derfor ofte, hvis der er et tykt dække af organisk materiale eller<br />

andre finkornede sedimenter i ådalen. Det grundvand, der afstrømmer<br />

overfladisk, vil ofte være dannet tæt på ådalen, men hvis det stammer fra kilder,<br />

kan det også være dannet regionalt. Vandets opholdstid i ådalsmagasinet er<br />

mellem timer til dage.<br />

Afstrømning fra dræn: Mange ådale er blevet drænlagte og har fået gravet<br />

grøfter. Vandets opholdstid i ådalen reduceres, når der er nedgravet dræn i<br />

forhold til den naturlige situation og kan forventes at være fra timer til dage<br />

(Nilsson et al., 2005: 52).<br />

Udover dette kan vandet også strømme ud af vandløbsbunden og strømme i ådalen<br />

langs vandløbet. Larkin & Sharp (1992) kalder denne form for strømning for underflow.<br />

I hydrologiske modeller bliver der ofte ikke taget højde for underflow i ådalen, hvilket<br />

kan være årsag til usikkerheder, da der derfor ikke bliver gjort rede for denne del af<br />

afstrømningen (Larkin & Sharp, 1992: 1608). Ifølge Larkin & Sharp er underflow<br />

dominerende, hvis vandløbets slyngningsgrad er under 1,3 og vandløbets hældning<br />

overstiger 0,0008 og ådalens sedimenter er lavpermeable (Larkin & Sharp, 1992: 1619).<br />

Der er en sammenhæng mellem hvilken type landskab et vandløb findes i og hvordan<br />

interaktionen mellem grundvand og overfladevand er. På smeltevandssletten, og i<br />

klitlandskabet, hvor sandede aflejringer dominerer, vil der ofte ske en høj grad af diffus<br />

udstrømning fra det overfladenære og det dybe grundvand. I morænelandskabet, på<br />

havaflejringer og de fladedækkende ferskvandsaflejringer vil der i højere grad ske en<br />

overfladisk afstrømning, hvilket skyldes, at sedimenterne i disse typer af landskaber<br />

21


Kapitel 2 Grundvand- og overfladevandsinteraktion<br />

ikke er nær så permeable. I morænelandskabet forstærkes dette yderligere af det<br />

udviklede relief.<br />

Den overfladenære afstrømning var tidligere den mest udbredte strømningsvej fra<br />

ådalen til vandløbet, men på grund af den omfattende dræning og grøftning er det i dag i<br />

højere grad afstrømningen fra dræn og dyb grundvandsafstrømning, der dominerer.<br />

Udover disse fire nævnte strømningsveje, kan grundvandet også forlade ådalen som<br />

fordampning, hvilket både sker i ådalen og fra vandløbet. I nogle ådale kan fordampningen,<br />

især i sommermånederne, være betragtelig, hvilket kommer af, at grundvandsspejlet<br />

står højt i ådalen og at der i ådalen findes et højt indhold af organisk materiale<br />

med en høj kapilaritet, hvor de kapilære kræfter trækker vandet fra den mættede zone og<br />

op til ådalens overflade, hvor det fordamper. Ådale tørrer derfor sjældent ud, da der er et<br />

stort potentiale for transport af vand (Burt et al., 2002: 130). Burt et al. (2002)<br />

karakteriserer således ådale som temporære vådområder, hvilket henviser til, at ådale<br />

ofte er præget af lavpermeable sedimenter og en lav højdegradient. Derfor bliver mange<br />

ådale fuldstændigt vandmættede i perioder, hvor der kommer meget nedbør. I nogle<br />

perioder, især på steder, hvor der ikke er så stor forskel på grundvandsstanden i ådalen<br />

og vandstanden i vandløbet, kan grundvandsspejlet falde så meget, at vand fra åen<br />

trænger ud på de omkringliggende områder (Burt et al., 2002: 130).<br />

I et studie af oplandet til Karup Å, der løber gennem et morænelandskab og en<br />

hedeslette, er det blevet vist, at det våde område i ådalen kan bruges som indikator for,<br />

hvor meget vand der strømmer ud i åen som underjordisk diffus udstrømning (Langhoff<br />

et al., 2005: 1-17). I dette opland blev det fundet, at omkring halvdelen af<br />

udstrømningen sker ved udsivning fra bunden. Hvis det våde område er lille, vil en stor<br />

del af udstrømningen ske via diffus udstrømning fra bunden. Er det våde område bredt,<br />

vil udstrømningen til vandløbet i stedet ske som overfladisk eller overfladenær<br />

afstrømning.<br />

En af de faktorer, der har stor indflydelse på strømningen fra grundvandsmagasinet til<br />

ådalen er, hvor permeabel bunden af vandløbet er. Dette er et produkt af, hvilke typer af<br />

aflejringer, der findes i ådalen og hvordan de er lejrede. Hvor god kontakten mellem<br />

grundvandsmagasinet og vandløbet er, bestemmes af bundens lækagekoefficient, der er<br />

den hydrauliske ledningsevne for bundens lag af sedimenter divideret med tykkelsen af<br />

lagene. Er vandløbet aggregerende, vil nogle af de partikler, der transporteres med<br />

vandløbet, blive afsat på bunden og begrænse (clogge) den hydrauliske ledningsevne for<br />

vandløbsbunden. Vandløbsbundens lækagekoefficient er derfor et mål, der varierer over<br />

selv små afstande. Den varierer med strømningshastigheden og vil dermed variere over<br />

tid, da strømningshastigheden varierer med tiden (Nilsson et al., 2003: 22). Endvidere<br />

22


har lækagekoefficienten betydning for, hvor stor en del af ådalens vand der strømmer<br />

som baseflow og hvor stor en del der strømmer som underflow.<br />

2.4 Opsummering<br />

I det foregående er der blevet set på, hvordan det komplekse sammenspil mellem<br />

grundvand og overfladevand er i en ådal, og nogle af de faktorer, der har indflydelse på<br />

denne interaktion er blevet gennemgået. For at besvare rapportens problemformulering<br />

er det, som vist i dette kapitel nødvendigt at kende områdets geologi og geomorfologi,<br />

idet det giver et overblik over hele oplandet og er baggrunden for den vandudveksling,<br />

der finder sted mellem ådalen og den tilgrænsende hydrologiske enhed. Derudover er<br />

det nødvendigt at have et mere detaljeret kendskab til geologien i selve ådalen, da den er<br />

baggrund for de hydrologiske processer der finder sted.<br />

Derudover er det nødvendigt, at der gennemføres undersøgelser af, hvordan<br />

udstrømningen af vand til vandløbet sker, da dette har en direkte indflydelse på,<br />

hvordan påvirkningen af Hølbækken kommer til at være. Dette skyldes, som tidligere<br />

nævnt, at vandet der strømmer ad de forskellige strømningsveje har forskellige<br />

opholdstid i ådalen og vandets opholdstid er betydende for, hvordan Hølbækken bliver<br />

påvirket i løbet af året. Det er således væsentligt, at den undersøgelse, der gennemføres,<br />

ikke kun fokuserer på hydrologien, som den ser ud på et bestemt tidspunkt, men også<br />

tager højde for, at denne ændres over tid. Dette skyldes, som tidligere nævnt, at ådale er<br />

kendetegnet ved, at der er en stor sæsonvariation i, hvordan vandudvekslingen sker.<br />

23


3<br />

Konceptuel model af<br />

projektområdet<br />

I dette kapitel udarbejdes der en konceptuel model over projektområdet, som skal danne<br />

input til grundvandsmodellen, der opsættes i det næste kapitel. Udarbejdelse af en<br />

konceptuel model indebærer, at der dannes et overblik over de vigtigste geologiske lag i<br />

projektområdet, hvorudfra der laves en geologisk model. Den geologiske model udgør<br />

udgangspunktet for opbygning af grundvandsmodellen i lag og linser. I forhold til den<br />

geologiske model er det nødvendigt at koble værdier for hydraulisk ledningsevne på de<br />

fundne geologiske lag, da strømningerne i projektområdet således bedst muligt kan<br />

gengives. Desuden er det nødvendigt at undersøge hydrologiske parametre for<br />

projektområdet, hvilket i dette kapitel vil gøres med udgangspunkt i vandbalanceligningen,<br />

igennem hvilken der ligeledes kommer input til grundvandsmodellen.<br />

Det er valgt at inddele dette kapitel i fire dele. Indledningsvis vil fokusområdet blive<br />

præsenteret ud fra en kortanalyse og de observationer, der blev gjort under feltturene.<br />

Dette danner grundlag for et kendskab til projektområdet, der både anvendes direkte i<br />

grundvandsmodellen og desuden bruges til at tolke på interaktionen mellem<br />

grundvandet i projektområdet og Hølbækken. Herefter vil projektområdets geologi blive<br />

undersøgt, både ud fra egne håndboringsundersøgelser og GEUS’ database. Dette<br />

munder ud i en geologisk model over området. Efter dette vil fokusområdets hydrologi<br />

undersøges i forhold til vandbalanceligningen, ud fra egne målinger af afstrømning og<br />

NJAs målinger af infiltrationen. Afslutningsvis vil projektområdets hydrogeologi blive<br />

undersøgt. Det drejer sig først og fremmest om værdier for de hydrauliske lednings-


Kapitel 3 Konceptuel model af projektområdet<br />

evner i ådalsmagasinet ved Hølbækken. Undervejs i kapitlet vil de forsøg, gruppen har<br />

udført, blive inddraget. I dette kapitel inddrages udelukkende forsøgenes resultater,<br />

mens metoder og fejlkilder findes som bilag bag i rapporten. Udregningerne der ligger<br />

til grund for resultaterne findes på den vedlagte Cd-rom.<br />

3.1 Præsentation af fokusområdet<br />

Først vil Hølbækkens historiske udvikling analyseres, på baggrund af kort fra tre<br />

forskellige tidsperioder. Dette illustrerer, hvilke menneskelige indgreb området har<br />

været udsat for, og giver noget af forklaringen på områdets nuværende udformning.<br />

Gruppen udførte flere feltture til området ved Hølbækken, for at indsamle data til<br />

grundvandsmodellen. Under disse feltture blev der desuden foretaget observationer af<br />

området, hvoraf de væsentligste i forhold til det videre arbejde med rapporten vil<br />

præsenteres i dette kapitel. Samlet set skal afsnittet sætte læseren bedre ind i rapportens<br />

fokusområde og danne baggrund for at de observationer gruppen gjorde under feltturene<br />

kan inddrages i det videre arbejde med den konceptuelle model.<br />

3.1.1 Kortanalyse<br />

For at få indblik i den menneskelige påvirkning af Hølbækken og Hølbækkens naturlige<br />

ændring over tid, er der foretaget en kortanalyse af området. Der er anvendt kort lavet i<br />

tre perioder, 1878-1887, 1924-1945 og det sidste blev færdiggjort i 1994. Kortene har<br />

samme målestok, centerkoordinat og orientering. Kortene kan ses på figur 3.1, 3.2 og<br />

3.3.<br />

På kortene ses Hølbækken midt i billedet, løbende mod sydøst til Lindenborg Å, der ses<br />

i nederste højre hjørne. Sammenlignes kortene ses det, at bækkens slyngning har ændret<br />

sig. Det midterste stykke af Hølbækken har en række slyngninger på det ældste kort,<br />

figur 3.3. Disse slyngninger er blevet mindre på kortet fra 1924-1945, hvor der også er<br />

oprettet 3 dræn, og i 1994 er slyngningerne helt væk. Dette viser, at bækken gradvist er<br />

blevet rettet ud på dette stykke. På den nedre del af Hølbækken ses det, at der på alle tre<br />

kort er omfattende dræning. Der har således fra senest engang sidst i 1800-tallet, har<br />

været et ønske om at dræne jorden. Dette har sandsynligvis været med henblik på at<br />

anvende arealerne til dyrkning eller afgræsning. Det kan herudfra tolkes, at der er en<br />

god kontakt mellem det overfladenære grundvandsmagasin og Hølbækken. Udløbet til<br />

Lindenborg Å har også gennem tiden ændret sig. På det ældste kort havde Lindenborg<br />

Å to parallelle åløb ved Hølbækkens udmunding. På kortet fra 1924-1945 er det<br />

nordligste af disse parallelle løb nærmest blevet til et dræn, og i 1994 er der kun et åløb.<br />

26


Figur 3.1. Kort over Hølbækken. Data fra DTK, nationalt kortværk færdigtegnet i 1994 (NJA, 2006).<br />

Kortets bredde er 2,2 km. Nord er opad.<br />

Figur 3.2. Historisk kort over Hølbækken, tegnet perioden 1924-1945 (NJA, 2006). Kortets brede er 2,2<br />

km. Nord er opad.<br />

Figur 3.3. Historisk kort over Hølbækken, tegnet perioden 1878-1887 (NJA, 2006). Kortets brede er 2,2<br />

km. Nord er opad.<br />

27


Kapitel 3 Konceptuel model af projektområdet<br />

Ved sammenligning af kortene, kan det ses, at vegetationen i området har ændret sig. I<br />

slutningen af 1800-tallet var området vest og nord for Hølbækken dækket af hede, mens<br />

der senere er kommet store skovarealer til og 1994 var der kun skov. Derudover er der<br />

sket en ændring i placeringen af dræn i Hølbækkens ådal i den periode kortene er fra.<br />

Der er derfor mulighed for, at der steder i Hølbækkens ådal findes ældre rørlagte dræn,<br />

der muligvis stadig fungerer, hvorfor området nu kan være drænet andre steder, end det<br />

er angivet på det nyeste kort. Denne antagelse understøttes af de observationer, der blev<br />

gjort under feltturen, hvor der flere steder blev fundet udløb fra drænrør langs<br />

vandløbet, der ikke var registreret på kortet. Dette vidner om, at der er en stor afstrømning<br />

ad denne vej. Under naturlige forhold ville dette have været overfladisk eller<br />

overfladenær afstrømning.<br />

Måles bredden af Hølbækken og dens ådal, ses det at ådalen er mere end fire gange så<br />

bred som åen. I forhold til GOI termologien vil det sige, at Hølbækken er ubegrænset og<br />

den kan dermed både kan være aggraderende og degraderende. I samme forbindelse er<br />

hældningen på vandløbet blevet målt til at være 0,0012 og slyngningsgraden er 1,17, og<br />

der er derfor mulighed for, at der i ådalen findes en komponent af underflow.<br />

3.1.2 Observationer<br />

Området omkring Hølbækken kan opdeles i to overordnede dele, ud fra relieffet i ådalen<br />

og i oplandet, figur 3.4. I den øvre del, hvor Hølbækken har sit udspring, løber<br />

vandløbet gennem en markeret ådal afgrænset af stejle skrænter. Der findes skov på den<br />

vestlige side af ådalen, mens den nordøstlige side består af afgræssede marker.<br />

Hølbækken starter som en smal grøft og bliver gradvist større og bredere. Stik øst for<br />

dens udspring findes en lille sø, figur 3.4, og der kan ses et rør i Hølbækken lige<br />

nedenfor, der sandsynligvis er overløb fra søen. Ved røret går Hølbækken fra at være<br />

smalt og nærmest uden vandføring, til at have en synlig vandføring. Vandløbet er her i<br />

underkanten af 0,5 m bredt, ca. 15 cm dybt. Ved kilde A på kortet bliver vandføringen<br />

markant større og vandløbet går fra at være en halv meter til en hel meter i bredden.<br />

Kilde A er en samling af små kildefremspring i åbunden. I kildefremspringene A, B, C,<br />

F og G ses lyst sand, der bobler på grund af vandets udstrømning, mens D, E og H er<br />

diffuse udstrømninger, der samles til små vandløb. Kilderne blev observeret flere steder<br />

i den øvre del, og afmærkningen på kortet dækker kun de største. Der er ikke fundet<br />

kilder på den nedre del. Vandet til kilderne stammer sandsynligvis fra det dybe<br />

regionale grundvandsmagasin, og er blevet transporteret til Hølbækken via sprækker i<br />

kalken, der her ligger tæt på jordoverfladen.<br />

28


Den nedre del af Hølbækkens ådal er også en del af Lindenborg Ås ådal og her flader<br />

området ud. I denne del er der, som nævnt i kortanalysen, indlagt en række dræn og<br />

markerne er inddraget til landbrugsformål. Den nedre del er præget af dræn, figur 3.4. I<br />

den nedre del af området blev der stedvis observeret okkerfarvning på brinkerne og i<br />

dræn på Hølbækkens østlige bred. Dette tyder på oxidering af jernmineraler som følge<br />

af dræningen. Vegetationen langs åen bestod i dette område af kraftig høj sivbevoksning,<br />

der bredte sig ud i åen.<br />

Figur 3.4. Hølbækkens ådal opdelt i øvre og nedre del. Opdelingen af Hølbækkens ådal i en øvre og<br />

nedre del vil blive brugt i den resterende del af rapporten. Data fra TOP10DK (2000), skræntfod fra<br />

(NIRAS, 2005).<br />

Langs hele Hølbækken blev der observeret flere våde områder, et særligt stort vådt<br />

område er observeret i ådalen, på den øvre del, hvilket er markeret på figur 3.4. Her<br />

stod grundvandet så højt, at der mellem tuerne kunne observeres et frit vandspejl i<br />

ådalen. Nedstrøms for dette område fortsatte et mere tørt område og hele den nedre del<br />

af ådalen var af denne karakter. Hølbækken er, som nævnt i kapitel 1 (Indledning),<br />

klassificeret som gyde- og yngelopvækstområde for laksefisk. Flere steder i åen kunne<br />

der ses indgreb for at forbedre levevilkårene for fiskene, blandt andet er der på flere<br />

strækninger udlagt sten i åen og åløbet er blevet manipuleret med udlægning af<br />

halmballer, som fremmer en slyngning. Der blev også observeret stimer af fiskeyngel<br />

flere steder.<br />

29


Kapitel 3 Konceptuel model af projektområdet<br />

3.2 Projektområdets geologi<br />

Formålet med denne rapport er at undersøge påvirkningen af en grundvandsindvinding i<br />

Volsted Plantage på Hølbækken. Det er i den forbindelse valgt at sætte op en numerisk<br />

grundvandsmodel for at undersøge indvirkningen. Før selve grundvandsmodellen bliver<br />

opstillet, er det vigtigt at kende den undergrund, modellen laves for. Ud fra en tolkning<br />

af undergrunden laves en geologisk model, der er en forenklet version af undergrundens<br />

opbygning. I dette kapitel bliver undergrunden i projektområdet undersøgt. Projektområdet<br />

bliver beskrevet ud fra aflejrings- og erosionsmiljøer, der har været dominerende<br />

siden øvre kridt og frem til i dag. Områdets geomorfologi vil blive beskrevet, da<br />

denne, som kapitel 2 (Grundvand- og overfladevandsinteraktion) viste, indirekte kan<br />

sige noget om interaktionen mellem grundvandet og overfladevandet. Hølbækkens ådal<br />

bliver beskrevet særskilt ud fra undersøgelser udført af gruppen, da dette område er af<br />

særlig interesse i forhold til problemformuleringen. Herefter vil grundvandsspejlets kote<br />

kortlægges, for at kunne placere grundvandsspejlets beliggenhed i forhold til de<br />

geologiske lag. Grundvandsspejlets beliggenhed er væsentligt i forhold til den<br />

geologiske model, da det sætter fokus på de lag, hvori grundvandsstrømningerne<br />

foregår. Som en opsummering af kapitlet bliver der opstillet en geologisk tolkningsmodel<br />

over projektområdet, der vil være udgangspunkt for det videre arbejde med<br />

grundvandsmodellen.<br />

3.2.1 Geologisk historie<br />

Øvre kridt<br />

Øvre kridt er den geologiske periode fra ca. 97 mio. år siden til ca. 65 mio. år siden.<br />

Dengang var Danmark, med undtagelse af Bornholm og Skagens Odde, gennem<br />

kridttiden dækket af et hav. Klimaet i øvre kridt var tørt og varmt og havet var præget af<br />

stille forhold. Partikelsedimentationen fra landområderne aftog, og aflejringerne fra<br />

denne periode består hovedsageligt af nedfældede kalkskaller fra organismer, der levede<br />

i havet. Eftersom laget af aflejringer blev tykkere, blev det kompakteret som det<br />

fremstår i dag (<strong>Rasmussen</strong>, 1975: 46-57; Fredericia, 1989: 60-67). Kalkoverfladens kote<br />

(m. o. h.) i projektområdet er afbildet på figur 3.5. Det ses, at kalkoverfladen ligger<br />

mellem 25 og 50 m. o. h. i hoveddelen af projektområdet og at overfladen falder ved<br />

Lindenborg Å i syd og Skriveren i vest. Kalkoverfladens relief blev udformet i de<br />

efterfølgende geologiske perioder både ved tektonisk aktivitet og ved erosion.<br />

Sænkningsfelterne i den prækvartære overflade ligger i de samme områder som under<br />

istiden blev udgravet af smeltevandsfloder, og det er sandsynligvis en sammenhæng<br />

mellem disse to forekomster, som beskrevet af Berthelsen (1987), i forhold til Aalborgs<br />

30


geologi. Den vestlige del af projektområdet domineres af en højtliggende kalkoverflade,<br />

også i forhold til terræn, mens den dykker til mellem 20-60 m under terræn i den østlige<br />

del. En højtliggende kalkoverflade er typisk for denne del af Nordjylland.<br />

Projektområdet ligger i den fennoskandiske randzone, hvor kalklaget generelt er<br />

højtliggende og desuden kommer op i en tykkelse på over 2.000 m (Thomsen, 1995:<br />

39). Der er i 2004 blevet foretaget en mikropalæontologisk undersøgelse af kalken i en<br />

boring midt i projektområdet (Hedeselskabet, 2006). Kalken i boringen stammede i hele<br />

dybden (165 m) fra den øvre del af del af kridttiden og var af typen skrivekridt, hvilket<br />

er den type, der sædvanligvis træffes overfladenært i Nordjylland (NJA, 2004).<br />

Figur 3.5. Den prækvartære overflades kote, denne kan anvendes som udtryk for kalkoverfladens<br />

beliggenhed da de tertiære aflejringer er borteroderet. Den prækvartære overflade er hentet fra GEUS<br />

(1994) og bygger på data fra boringer.<br />

Tertiær<br />

I den efterfølgende periode, tertiær, ca. 65 mio. år siden til ca. 2 mio. år siden, blev<br />

klimaet fugtigere og partikelsedimentation fra landområderne øgede igen. Sand, ler og<br />

mergel blev aflejret i denne periode, men på grund af isens senere erosion af<br />

aflejringerne, er tertiæraflejringerne i projektområdet blevet eroderet væk (Thomsen,<br />

1995: 32-43), og de vil derfor ikke blive beskrevet nærmere her. Under tertiærtiden<br />

opstod der imidlertid en række salthorste, især lokaliseret i et øst-vest gående bånd<br />

gennem Jylland, hvor kalklaget er af betydelig tykkelse. Salthorste opstår når<br />

evaporitlaget fra permtiden (ca. 270 - 230 mio. år siden) kommer under højt tryk fra<br />

31


Kapitel 3 Konceptuel model af projektområdet<br />

overliggende lag med højere densitet, bliver plastiske og presses op (Galsgaard, 1998:<br />

32). Tættest på projektområdet ligger en salthorst ved Suldrup, ca. 16 km vest for<br />

Volsted (<strong>Rasmussen</strong>, 1975: 78-79). Denne salthorst kan have haft betydning for<br />

kalklaget i projektområdet ved, da den kan have medført sprækkedannelse i kalken.<br />

Kvartær og nyere tid<br />

I kvartærperioden ændrede klimaet sig igen, og perioden var præget af hurtigt skiftende<br />

klimaforhold fra kolde perioder, hvor gletschere dækkede hele eller dele af Danmark, til<br />

milde interglaciale perioder, hvor isen smeltede (<strong>Rasmussen</strong>, 1975: 80). Isen formede<br />

landskabselementer som tunneldale, bakkeøer og moræneaflejringer og satte sit tydelige<br />

præg på landskabet, som det ser ud i dag. På figur 3.6 ses de forskellige aflejringsmiljøer,<br />

der i dag kan genfindes i projektområdet og det omkringliggende område.<br />

Hoveddelen af projektområdet består af yngre moræneaflejringer afsat af gletschere.<br />

Moræne er usorteret materiale med stor spredning i partikelstørrelsen. Langs<br />

projektområdets sydlige og vestlige grænse, langs Lindenborg Å og Skriveren, ses der<br />

tunneldale med aflejringer fra Yoldiahavet og Litorinahavet. Tunneldale dannes af<br />

smeltevand, der strømmer ud under gletscheren. Dette vand kan strømme med stor<br />

hastighed og således have en effektfuld eroderende kraft. Aflejringerne i tunneldale er<br />

sorterede og størrelsen på partiklerne er bestemt af smeltevandsflodens strømningshastighed.<br />

I de interglaciale perioder steg havets niveau grundet den fra isen frigivne<br />

vandmængde, og landområder blev sat under vand. Yoldiahavet bredte sig ind over<br />

Nordjylland for ca. 15.000 år siden (Larsen, 1989: 37). Efter isen havde smeltet fra<br />

Danmark fra sidste gang bredte Litorinahavet, også kaldet Stenalderhavet, sig ind over<br />

Danmark for ca. 7.000 – 4.000 år siden. Aflejringer fra Litorinahavet findes i Lindenborg<br />

Å og Østerå ådale i projektområdet og for øvrigt også nord og øst for<br />

projektområdet. På figur 3.6 ses det, at projektområdet ligger på en ø af moræneaflejringer,<br />

omgivet af aflejringer fra Litorinahavet. Disse to typer af aflejringer er<br />

kendetegnet ved en høj grad af overfladisk afstrømning, jf. kapitel 2 (Grundvand- og<br />

overfladevandinteraktion), hvilket også blev observeret i fokusområdet. Sammenholdes<br />

figur 3.6 med GOI typografien ses det, at Hølbækken dels strømmer gennem et<br />

morænelandskab og på havaflejringer. Disse to aflejringstyper er kendetegnet af en høj<br />

grad af overfladisk og overfladenær afstrømning, hvilket også blev observeret i<br />

området.<br />

32


Figur 3.6 Landskabselementerne i projektområdet og omegn (DJF-geodata,, 2006).<br />

3.2.2 GEUS boredata<br />

Med baggrund i beskrivelsen af projektområdets geologi er undergrunden blevet<br />

undersøgt ud fra boringer fra GEUS’ Jupiter boredatabase. Tolkningen er gjort ud fra<br />

250 boringer, hvoraf 110 boringer findes inden for projektområdet, figur 3.7. Som<br />

figuren viser, er der en stor samling boringer i den vestlige del af projektområdet, men<br />

desværre fokuserer de mest på kalklaget og de overliggende lag er således dårligt<br />

beskrevet. Der er ligeledes mange boringer i den nordøstlige del af projektområdet, men<br />

disse er for hovedpartens vedkommende ikke særlige dybe, men har til gengæld en<br />

detaljeret beskrivelse af undergrunden. I den midterste del af projektområdet, tæt på<br />

fokusområdet, er der kun enkelte boringer. På trods af dette giver boringerne en god<br />

mulighed for at få indblik i projektområdets undergrund. For at øge overskueligheden af<br />

dataene og lette tolkningen af dem, er projektområdet inddelt i 22 profillinjer med 200<br />

m buffer til hver side. Samme profiler bliver brugt i digitaliseringen af projektområdet i<br />

forbindelse med grundvandsmodellen.<br />

33


Kapitel 3 Konceptuel model af projektområdet<br />

Figur 3.7. Boreprofiler for projektområdet. Det udvalgte profil analyseres herunder. Boringer fra GEUS 2<br />

(2006).<br />

I figur 3.8 er vist det fremhævede profil i figur 3.7. På profilet er terrænoverfladen vist<br />

sammen med den prækvartære overflade. Ved at sammenligne dataene fra de enkelte<br />

boringer med den prækvartære overflade og terrænoverfladen ses visse uoverensstemmelser.<br />

Dette kan til dels skyldes, at den prækvartære overflade er tegnet direkte<br />

efter profillinien, mens borehullerne kan ligge op til 200 m derfra. Terrænoverfladen<br />

kan ligeledes variere. Uoverensstemmelserne kan også skyldes unøjagtigheder i<br />

dataene. Kalken er i nogle boringer registreret som skrivekridt, mens andre blot er<br />

defineret som kalk, men det antages, at al kalken i projektområdet er skrivekridt.<br />

Profilet går fra morænen i det nordøstlige hjørne af projektkområdet mod Lindenborg Å<br />

ådal. Det ses, at profilet hovedsageligt udgøres af boringer med lille dybde, hvilket er et<br />

generelt problem med boredata fra den nordøstlige del af området. Ud fra viden om<br />

dannelsesmiljø og kalkoverfladens kote antages det, at boringerne i morænen giver et<br />

rimeligt billede af de kvartære aflejringer. I morænen findes hovedsagelig sand. Der er<br />

registreret hyppige skift mellem sand og grus med forekomster af glacialt<br />

smeltevandssand. Dette tyder på, at smeltevand har strømmet over projektområdet og<br />

har overlejret noget af morænen med smeltevandssand. Der forekommer ler i den<br />

nordøstlige del af projektområdet. På figur 3.9 ses i hvilke boringer, der er fundet ler.<br />

Lerforekomsterne er koncentreret i det nordøstlige hjørne af projektområdet, hvilket<br />

formentlig er moræneler.<br />

34


Figur 3.8. Boreprofil. Boredata fra GEUS 2 (2006), prækvartær overflade fra GEUS (1994), terrænmodel<br />

fra DDH (1995).<br />

Figur 3.9. Områder hvor det skønnes at der findes lerlinser. Boredata fra GEUS 2 (2006).<br />

35


Kapitel 3 Konceptuel model af projektområdet<br />

På den del af boreprofilet, der viser Lindenborg Ås ådal (ved ca. 4.000 – 6.000 m) ses et<br />

sænkningsfelt i den prækvartære overflade, figur 3.8. Sænkningsfeltet genfindes langs<br />

hele ådalen. I ådalen er der opfyldningen registreret som tørv, mergel og sandforekomster,<br />

men der er få boringer at tolke geologien i ådalen ud fra, hvorfor det er<br />

svært at lave en generel tolkning af ådalen ud fra disse boringer.<br />

Fejlkilder i forbindelse med tolkning af boredata<br />

Boredata er punktdata, og når undergrunden tolkes herudfra er det nødvendigt at være<br />

opmærksom på hvilke generaliseringer, der bliver gjort i den forbindelse. Undergrundens<br />

opbygning kendes kun, hvor boringerne er foretaget. Mellem borepunkterne<br />

bliver der interpoleret. Dette kan i områder med få boringer medføre betragtelige<br />

fejlkilder i tolkningen. Det er derfor nødvendigt at tolke, hvorvidt en forekomst i<br />

en/flere boringer er et sammenhængende lag, der kan interpoleres eller om forekomsten<br />

er en linse, og dermed har begrænset horisontal udbredelse. Ligeledes må resultatet af<br />

interpoleringen overvejes i forhold til det antal boringer, den er lavet ud fra. En måde at<br />

imødekomme disse fejlkilder på er, at benytte geofysiske metoder til at undersøge<br />

området. Især i den nordøstlige del af projektområdet er tolkningen foretaget på<br />

baggrund af få boringer og er af denne grund specielt unøjagtig.<br />

3.2.3 Håndboringer i fokusområdet<br />

Ådalene i projektområdet adskiller sig fra morænen i opbygning ved, at det er en anden<br />

type aflejringer, der udgør opfyldningen end aflejringerne i resten af projektområdet.<br />

Fokusområdet i denne rapport er Hølbækken, hvor der ikke findes GEUS boringer. Det<br />

er væsentligt at undersøge ådalens geologi, og danne input til en mere præcis geologisk<br />

model for dette område, da det er i Hølbækkens ådal, at vandudveksling mellem<br />

grundvandet og overfladevandet sker. Undersøgelsen tager udgangspunkt i boredata<br />

indhentet af NIRAS (Jensen et al., 2005). Efter undersøgelse af disse data har gruppen<br />

udført en række boringer med håndbor, hvor lagdelingen i området og udbredelsen af de<br />

forskellige lag er fundet. NIRAS’ undersøgelse viste, at den nedre del af området er<br />

forholdsvis homogent, hvorfor der er indsamlet få prøver herfra. I den øvre del er<br />

undergrunden mere heterogen, bl.a. er der registreret afslutningen på et gytjelag. Da det<br />

var forventet, at gytjelaget har en lille hydraulisk ledningsevne, og dermed er væsentlig<br />

i forhold til grundvandsstrømningerne i området, blev boringerne koncentreret i denne<br />

del.<br />

Ud fra boringerne er der opstillet fire profiler, tre tværgående og et langsgående. For<br />

overskuelighedens skyld opstilles det langsgående profil først kun for den øvre del af<br />

området, mens den nedre del undersøges særskilt. Placeringen af boringerne og<br />

profillinierne kan ses på figur 3.10.<br />

36


Figur 3.10. Placering af boringerne og profillinierne i projektområdet.<br />

Langsgående profil<br />

På figur 3.11 ses det langsgående profil fra den øvre del. Profilet er standset ved<br />

grænsen mellem den øvre og nedre del af området. Det ses, at kalklaget ligger i højest<br />

kote i den øvre del af Hølbækken og daler nedstrøms. Dette kunne tyde på, at kalklaget<br />

på et tidspunkt er blevet eroderet af smeltevand fra gletscheren, som har dannet en ådal,<br />

hvilket kunne være oprindelsen til den ådal, som findes i dag. Samtidigt er der sandsynligvis<br />

sket en indsynkning af ådalene. I de fleste borehuller findes et sandlag ovenpå<br />

kalken. Sandets kornstørrelsesfordeling varierede, hvilket tyder på forskellige aflejringsmiljøer.<br />

I forhold til de geologiske perioders aflejringsmiljøer, beskrevet ovenfor, er det<br />

sandsynligt, at der både findes smeltevandsaflejringer og aflejringer fra Litorinahavet<br />

blandt sandet fundet i ådalen. Der er ikke i boringerne skelnet mellem disse forskellige<br />

typer sand og fremover vil der blive henvist samlet til disse forekomster som sand.<br />

Sandet har sin største vertikale udbredelse ca. 400 m fra profilets start og stopper<br />

tilsyneladende ca. 150 m fra profilets start. Dette kan skyldes en øvre afgrænsning af<br />

sandets aflejringsmiljø. I 4 af de øvre boringer er der fundet et stenlag, sandsynligvis<br />

afsat af smeltevand med stærkere strømningsforhold under aflejring i forhold til sandet.<br />

Der er ikke fundet kalk i disse 4 boringer, sandsynligvis fordi det har været svært at<br />

bore forbi stenlaget.<br />

37


Kapitel 3 Konceptuel model af projektområdet<br />

Ca. 400 m fra profilets start, starter et gytjelag, i de fleste boringer ovenpå sandet. Gytje<br />

er en marin aflejring, der består af nedbrudt organisk materiale med et højt indhold af<br />

uorganisk materiale (Larsen et al., 1988: 52). Dette viser, at området ved Hølbækken<br />

har været delvis dækket af stillestående vand i en periode, for at det finkornede<br />

materiale har kunnet sedimentere. Gytjelaget bliver gradvist tykkere i den nedre del af<br />

ådalen. I den øvre del er der ikke fundet gytje i alle boringerne, hvilket sandsynligvis<br />

skyldes, at gytjelaget ikke har en lineær afgrænsning og nogle boringer er derfor lavet<br />

lige udenfor gytjelaget. Der er også mulighed for, at der kan være blevet aflejret gytje i<br />

hele ådalen og at laget efter aflejringen er blevet borteroderet. Øverst i boringerne ligger<br />

tørvelaget, hvis overflade hælder svagt nedstrøms mod Hølbækkens udløb mod<br />

Lindenborg Å. Tørv består af dårligt nedbrudt organisk materiale, der er dannet i fugtige<br />

eller sumpede områder, hvor stillestående vand har forhindret ilttilførsel (Fobian, 1984:<br />

60). Fugtige områder blev under feltturen observeret flere steder langs Hølbækken og<br />

det er sandsynligvis i dag den samme type aflejringsmiljø.<br />

I afsnit 3.1.2 (Observationer) blev det beskrevet, at der kun på den øverste del forekom<br />

kilder. Kilderne tyder på, at der findes præfererede strømningsveje i dette område.<br />

Sandsynligvis er der sprækker i kalken, hvor vandet kan strømme. Stenlaget fundet i<br />

den øvre del af profilet kan ligeledes give mulighed for præferentiel strømning. Da der<br />

forekommer kildefremspring i områder, hvor der er registreret gytje, kunne dette tyde<br />

på, at der findes strømningsveje gennem gytjen, eller at vandet i kilderne er overfladenær<br />

diffus udstrømning, der strømmer ovenpå gytjen.<br />

38<br />

Kote (m)<br />

7<br />

6<br />

5<br />

4<br />

3<br />

2<br />

1<br />

0 100 200 300 400 500 600 700<br />

Længde (m)<br />

Muld/tørv<br />

Figur 3.11. Langsgående profil på den øvre del. Den horisontale akse starter øverst i profilet og følger<br />

Hølbækkens strømningsretning. Der er en overhøjde på 100.<br />

Gytje<br />

Sand<br />

Kalk<br />

Sten


Tværprofil 1<br />

Lagopdelingen i borehullerne i profil 1 kan ses i figur 3.12. Kalklaget danner en dal,<br />

hvilket harmonerer med antagelsen om, at kalken under ådalen er sunket ned og at<br />

smeltevand har eroderet bunden og at denne gamle dal, danner basis for den nuværende<br />

ådal. Ovenpå på kalken er der fundet sand i fire af boringerne, men i den ene boring (nr.<br />

4 fra venstre) er det kun fundet et tyndt lag. Dette kunne tyde på, at der har sket en<br />

erosion af sandet her, der har efterladt en fordybning i terrænet. I boring 3 og 4 fra<br />

venstre er der fundet gytje, hvilket viser, at dette område på et tidspunkt har været<br />

dækket af stillestående vand. Gytjelaget i boring 4 er tykkere, mens det i boring 3 er<br />

tyndere og ligger lavere, hvilket kunne tyde på, der er sket en erosion af gytjen efter<br />

aflejringen.<br />

Tværprofil 2<br />

Lagopdelingen i borehullerne i profil 2 kan ses i figur 3.13. To af boringerne i dette<br />

profil blev stoppet før kalklaget blev nået, men resultaterne fra de andre boringer tyder<br />

alligevel på, at kalklaget også her danner en dal, som måske ligger dybere. Det er<br />

vanskeligt at vurdere udbredelsen af sandlaget, da to af boringerne er stoppet netop som<br />

dette lag er nået. Der ses dog en stor forskydning af overfladen på sandlaget, hvilket<br />

muligvis skyldes, at løbende vand har eroderet i midten af ådalen, og skåret ned i<br />

sandlaget. Som ved tværprofil 1 er der ikke fundet sand i boring nr. 4 fra venstre,<br />

sandsynligvis pga. erosion i denne del af tværsnittet.<br />

Tværprofil 3<br />

Lagopdelingen i borehullerne i profil 3 kan ses i figur 3.14. Profil 3 ligger umiddelbart<br />

før overgangen til den nedre del af Hølbækken, figur 3.10. Kalklaget danner en dal,<br />

ligesom i de to beskrevne profiler. Sandlaget har forskellig tykkelse og overfladen af<br />

dette danner en stejl dal. Gytjen er i dette profil væsentlig tykkere end den var i de to<br />

foregående profiler, hvilket tyder på, at gytjelagets tykkelsen øges nedstrøms i ådalen<br />

som forventet.<br />

Nedre del<br />

Figur 3.15 viser de fem boringer udført langs Hølbækken i den nedre del af området.<br />

Boringernes overflade viser, at topografien her flader ud. Den midterste boring er taget<br />

ca. 200 m fra Hølbækken, mens de andre boringer ligger tættere på vandløbet, figur<br />

3.10. Der er kun boret forbi gytjelaget i to af boringerne, mens de sidste tre boringer er<br />

stoppet i gytjelaget. Der er fundet sand under gytjen i begge de to boringer. Det ses, at<br />

gytjelagets tykkelse er tiltaget. Det antages, at disse boringsresultater er repræsentative<br />

for ådalen i denne del af Hølbækken og for Lindenborg Ås ådal, da der i alle boringerne<br />

er fundet ca. 2 m tørv, der overlejrer gytjen og at gytjen i alle boringerne findes som en<br />

væsentlig del af opfyldningen.<br />

39


Kapitel 3 Konceptuel model af projektområdet<br />

40<br />

Kote (m)<br />

7<br />

6<br />

5<br />

4<br />

3<br />

2<br />

1<br />

0<br />

0 20 40 60 80<br />

-1<br />

Længde (m)<br />

Figur 3.12. Tværprofil 1. Profilet ses fra vest til øst. Der er en overhøjde på 10.<br />

Kote (m)<br />

7<br />

6<br />

5<br />

4<br />

3<br />

2<br />

1<br />

0<br />

0 20 40 60 80<br />

-1<br />

Længde (m)<br />

Figur 3.13. Tværprofil 2. Profilet ses fra vest til øst. Der er en overhøjde på 10.<br />

Muld/tørv<br />

Gytje<br />

Sand<br />

Kalk<br />

Hølbækken<br />

Muld/tørv<br />

Gytje<br />

Sand<br />

Kalk<br />

Hølbækken


Kote (m)<br />

7<br />

6<br />

5<br />

4<br />

3<br />

2<br />

1<br />

0<br />

0 20 40 60 80<br />

-1<br />

Længde (m)<br />

Figur 3.14. Tværprofil 3. Profilet ses fra vest til øst. Der er en overhøjde på 10.<br />

Kote (m)<br />

5<br />

4<br />

3<br />

2<br />

1<br />

0<br />

0 10 20 30 40 50 60<br />

-1<br />

-2<br />

-3<br />

-4<br />

-5<br />

Længde (m)<br />

Muld/tørv<br />

Gytje<br />

Sand<br />

Kalk<br />

Hølbækken<br />

Muld/tørv<br />

Figur 3.15. Langsgående profil på den nedre del. Den horisontale akse følger Hølbækkens<br />

strømningsretning. Der er en overhøjde på 7,5.<br />

Gytje<br />

Sand<br />

Kalk<br />

41


Kapitel 3 Konceptuel model af projektområdet<br />

Gytjeafgrænsning<br />

På figur 3.16 er der markeret, hvilke boringer, der indeholdt gytje. Herudfra tolkes det,<br />

at der er en afgrænsning af gytjelaget. Analysen af profilerne har vist, at gytjelaget<br />

bliver tyndere længere oppe i ådalen, og det ses, at det stopper midt i ådalen. Dette<br />

vidner om at udbredelsen af det stillestående vand, hvori gytjen blev aflejret ikke har<br />

gået helt op i den øvre del af ådalen, men hovedsageligt har været lokaliseret på den<br />

nedre del. Gytjelagets udbredelse har betydning for grundvandsstrømningerne i ådalen,<br />

da gytje har en lav hydraulisk ledningsevne og således fungerer som en aquitard. Det<br />

må forventes, at kildefremspringene kan sættes i forbindelse med gytjelaget, idet der er<br />

registrerede kilder i de områder, hvor der er registreret brud på gytjelaget. Det kan<br />

derfor tænkes, at bruddet på gytjelaget muliggør strømning af vand op gennem opfyldningen<br />

i ådalen. Udover de kilder, der kan sættes i forbindelse med brud på gytjelaget,<br />

findes der et par mindre kilder, som muligvis skyldes, at der sker en overfladenær diffus<br />

udstrømning oven på gytjen. Gytjen er derfor i Hølbækkens ådal en vigtig parameter i<br />

forhold til strømningen af vand, idet laget begrænser kontakten mellem vandløbet og det<br />

regionale grundvandsmagasin i de områder, hvor laget findes.<br />

Figur 3.16. Boringer, hvor der er fundet gytje er indtegnet og afgrænsningen af gytjelaget er herudfra<br />

blevet skønnet.<br />

Det kan også tænkes, at gytjelaget ikke er blevet aflejret ens i hele ådalen grundet<br />

kilder, der medfører strømninger i det eller stillestående vand. Dette forudsætter, at de<br />

42


kilder, der findes i Hølbækkens ådal, strømmer gennem lag, der er højpermeable,<br />

eksempelvis sprækker i kalken. Strømningen gennem kilden vil forhindre sedimentationen.<br />

Der er også den mulighed, at noget af opfyldningen i ådal er blevet fjernet<br />

efterfølgende og at dette er en medvirkende årsag til de mange kilder. Dette understøttes<br />

af, at ådalen gennem de seneste hundrede år er blevet drænlagt. Særlig i den øvre del af<br />

ådalen, hvor gytjelagte ligger højt og er tyndt, kan dette have ført til, at der er skabt<br />

strømningsveje gennem gytjelaget.<br />

3.2.4 Grundvandsspejlets beliggenhed<br />

På figur 3.17 er grundvandskoten sammenlignet med den prekvartære overflade. Det<br />

ses, at grundvandsspejlet i størstedelen af projektområdet ligger under kalkoverfladen.<br />

Kalken udgør således en væsentlig del af grundvandsmagasinet i projektområdet. I<br />

ådalene langs Lindenborg Å og Skriveren ligger kalkoverfladen dybere end grundvandsspejlet.<br />

Grundvandsspejlet ligger i den vestlige del af projektområdet, bl.a. i Volsted<br />

Plantage, i kalken. Fokus i grundvandsmodellen er grundvandsmagasinet, hvorfor en<br />

præcis kortlægning af morænen her ikke er af så stor betydning. I den østlige del af<br />

projektområdet ligger grundvandsspejlet i morænen, hvorfor lerlinser her er af større<br />

betydning og tolkningen her er gjort mere omhyggeligt. Kalk er et porøst materiale, der<br />

kan indeholde store mængder vand, men porestørrelsen er lille, hvorfor den hydrauliske<br />

ledningsevne for sammenpakket kalk er lille. Sprækkedannelse i kalken vil derfor have<br />

en kraftig indflydelse på kalkens samlede hydrauliske ledningsevne. Sprækkerne vil<br />

fungere som makroporer og vandet vil kunne transporteres hurtigt gennem disse.<br />

Beliggenheden og udbredelsen af sprækkerne er derfor en væsentlig faktor til at<br />

bestemme, hvordan grundvandet strømmer i kalk og dermed også i projektområdet. I<br />

den boring Hedeselskabet foretog for Nordjyllands Amt i 2004, blev der udført en<br />

flowlog, som er en måling over hvor hurtigt vandet strømmer ud i en boring, hvis der<br />

sker en oppumpning af vand og i hvilke lag, strømningen foregår. I den udførte flowlog<br />

kom det frem, at der er en jævn udstrømning af vand til boringen til en dybde af 80<br />

meter, hvorefter den aftager og bliver ubetydelig, bilag A. Kalken i projektområdet<br />

antages derfor at være opsprækket de øverste ca. 80 meter.<br />

Det er også interessant i forhold til grundvandsmodelleringen at kende den nedre<br />

grænse for grundvandsmagasinet. Den nedre grænse af gytjelaget kan angives ved at<br />

bestemme dybden til det salte grundvand. I 2004 blev der gennemført 23 TEM målinger<br />

i Volsted Plantage (Watertech, 2005), hvor konklusionen var, at det salte grundvandsspejl<br />

befandt sig i en dybde af omkring 140 meter, hvorfor dette kan anvendes som den<br />

nedre afgrænsning af grundvandsmagasinet.<br />

43


Kapitel 3 Konceptuel model af projektområdet<br />

Figur 3.17. Grundvandsspejlets overflade i forhold til den prækvartære overflade. Positive tal angiver, at<br />

grundvandsspejlet ligger over den prækvartære overflade og negative angiver, at grundvandsspejlet ligger<br />

under den prækvartære overflade. Prekvartæroverflade fra GEUS (1994) potentiale kort fra (NJA, 2005).<br />

3.2.5 Geologisk model<br />

Ud fra de foregående afsnit, hvor undergrunden i projektområdet generelt og særskilt i<br />

Hølbækkens ådal er blevet undersøgt, vil den geologiske tolkningsmodel opsættes. Det<br />

er valgt at opdele undergrunden i tre hoveddele: kalklaget, morænen og ådalene, figur<br />

3.18. Opdelingen er begrundet i analysen af undergrunden, hvor det kommer frem, at<br />

der er markant forskellige hydrogeologiske forhold, i de tre dele. Endvidere er datagrundlaget<br />

for de to områdetyper forskelligt, som beskrevet ovenfor, hvorfor det er<br />

hensigtsmæssigt at adskille dem. Ådalen er angivet som det område af modellen, der<br />

ligger langs Lindenborg Å og Hølbækken. Det er blevet afgrænset ved hjælp af<br />

højdekurver og områder med dræn, figur 3.19.<br />

44


Tørv<br />

Gytje<br />

Ådal Moræne<br />

Figur 3.18. Principiel opbygning af undergrunden, anvendt i modelopsætningen.<br />

Figur 3.19. Opdeling af projektområdet i ådal og moræne.<br />

Morænemateriale<br />

Kalk<br />

Kalk (saltvandsgrænse)<br />

45


Kapitel 3 Konceptuel model af projektområdet<br />

46<br />

Kalklaget - Grundet dannelsesmiljøet, bliver kalklaget tolket til at være et<br />

sammenhængende lag. Kalklaget er repræsenteret i hele projektområdet og<br />

udgør desuden grundvandsmagasinets nedre afgrænsning. Kildeforekomsten i<br />

den øvre del af Hølbækkens ådal tyder på, at der her findes præfererede<br />

strømningsveje for grundvandet, muligvis er det her tale om sprækker i kalken.<br />

Kalkens dybde er sat til kote -140, da der her findes salt grundvand og<br />

strømningerne derfor ikke relevante i forhold til denne rapports grundvandsmodel.<br />

Morænen - Ud fra GEUS boringer og de præsenterede kort er det blevet tolket,<br />

at morænen hovedsageligt består af sand, med forekomster af ler, silt og grus,<br />

hvilket er forventeligt, når det er moræne. For at repræsentere denne tolkning er<br />

det valgt at tolke morænen som bestående af et sandlag. Dette skal repræsentere<br />

både sandet og gruset, hvor der ikke er så stor forskel i den hydrauliske<br />

ledningsevne. For at angive forekomster af ler og silt, der har en lavere<br />

hydraulisk ledningsevne end sand, er det valgt at indlægge linser. Da<br />

grundvandsspejlet i hoveddelen af området ligger i kalken, er der lavere krav til<br />

detaljeringsgraden af morænetolkningen.<br />

Ådalen - Der er udført håndboringer i Hølbækkens ådal, hvilket giver et<br />

detaljeret tolkningsgrundlag for dette område. Det er nødvendigt at generalisere<br />

lagfølgen for Lindenborg Ås ådal herudfra, da der kun er et ringe antal boringer i<br />

GEUS’ database i dette område. I fokusområdet blev det fundet, at der ligger<br />

tørv i det øverste 2 meter, hvorunder der er et gytjelag af varierende tykkelse. I<br />

nogle af boringerne blev der fundet marint sand, hvorfor det langs ådalen over<br />

kalken, antages, at der ligger et sandlag.<br />

3.3 Projektområdets hydrologi<br />

Som del af den konceptuelle model er det nødvendigt at finde parameteren der indgår i<br />

vandbalanceligningen for området. Parametrene skal bruges som input til<br />

grundvandsmodellen, der opsættes senere i denne rapport. De parametre, der indgår i<br />

vandbalanceligningen er nedbøren, fordampningen, oppumpning, ændring i magasinering<br />

og den overjordiske og underjordiske afstrømning. Det antages, at der kan ses<br />

bort fra magasineringen i projektområdet når dette betragtes over en årrække, da<br />

magasineringen er minimal i forhold til den resterende vandudveksling. Det er ikke<br />

muligt at måle den underjordiske afstrømning. De opsatte nulflux grænser i den nordlige<br />

del af projektområdet sikrer imidlertid, at vandudveksling ikke sker her og<br />

trykgrænserne i den resterende del sikrer, at vandudvekslingen i dette område kan


eregnes i grundvandsmodellen. Derfor vil nedbør, fordampning, oppumpning og<br />

overfladisk afstrømning i det følgende blive undersøgt for projektområdet. På grund af<br />

den tætte sammenhæng mellem nedbør og fordampning, og problemer med at skaffe<br />

data, vil disse to bliver undersøgt sammen.<br />

3.3.1 Overfladisk afstrømning<br />

Først vil afstrømningen fra Hølbækken blive undersøgt, for at kortlægge interaktionen<br />

mellem grundvandet og overfladevandet i projektområdet og for at give et bud på den<br />

samlede afstrømning fra Hølbækken. Hølbækkens vandføring dækker ikke hele<br />

projektområdet og samtidig er der ikke kontinuerte målinger for et helt år. Derfor<br />

undersøges afstrømningen fra oplandet til Lindenborg Å bagefter, da Hølbækkens<br />

opland er en del af dette opland. Afslutningsvis foretages en sammenligning mellem de<br />

to resultater.<br />

Oplandet til Hølbækken<br />

Vandføringen i Hølbækken udgør kun en del af den samlede vandføring i<br />

projektområdet. For at få et overblik over, hvor stor andel vandføringen i Hølbækken<br />

udgør af den samlede vandføring, er det hydrologiske og det topografiske opland til<br />

Hølbækken blevet afgrænset, figur 3.20.<br />

Figur 3.20. Det topografiske og hydrologiske opland til Hølbækken, optegnet på et kort over<br />

grundvandskoten. Grundvandskoten er interpoleret ud fra data fra GEUS boredata (GEUS 2, 2006) og<br />

koter på åer og dræn. Åer og dræns placering er taget fra TOP10DK (2000), og har fået koter vha. DHH<br />

(1995).<br />

47


Kapitel 3 Konceptuel model af projektområdet<br />

Forskellen mellem det topografiske og det hydrologiske opland kan dels skyldes<br />

usikkerhed i de anvendte data, hvor f.eks. potentialekortet er interpoleret ud fra en<br />

række punkter og derfor muligvis ikke er så præcist. Punkterne potentialekortet er lavet<br />

på baggrund af, er der også usikkerhed omkring, da de har forskellig alder, den ældste<br />

boring er eksempelvis fra 1930. Derudover varierer potentialet ca. 2 m i løbet af året<br />

(Søndergaard, 2003: B4). På trods af disse usikkerheder tyder det på, at der er en forskel<br />

på de to oplande, hvilket må skyldes undergrundens opbygning i området. Eksempelvis<br />

findes der lerlinser i projektområdet, der kan være med til at ændre grundvandsstrømningen,<br />

så den nedbør der falder i det topografiske opland strømmer til<br />

Hølbækken.<br />

Vandføringen i Hølbækken<br />

For at undersøge den samlede vandføring og variationen i vandføringen langs Hølbækken<br />

er vandføringen målt i en række tværsnit langs åen over en periode på to dage.<br />

Placering af vandføringsmålingerne kan ses på figur 3.21. Forud for målingerne var der<br />

opholdsvejr og ca. 0 °C, og derfor kan det antages, at målingerne tilnærmelsesvis beskriver<br />

en baseflowsituation. For at få et indtryk af variationerne i vandføringen, blev<br />

vandstanden også målt over en periode på 52 dage. Dette er ikke nok til at vurdere<br />

vandføringen over et år, men giver alligevel et indblik i variationer i bækkens vandføring.<br />

Udførelsen af vandføringsmålingerne, vejrforhold i tiden op til målingerne samt<br />

fejlkilder forbundet med forsøget kan ses i bilag B. Udregninger findes på den vedlagte<br />

Cd-rom.<br />

Figur 3.21. Vandføringsmålingernes placering i Hølbækken.<br />

48


Der blev foretaget vandføringsmålinger ved 9 punkter i Hølbækken, hvor punkt 6 og 7<br />

blev målt igen efter ca. to måneder. Derudover blev vandføringen målt i kilde E, figur<br />

3.4, mellem måling 4 og 5. Resultaterne for vandføringsmålingerne i de forskellige<br />

punkter kan ses på figur 3.22.<br />

Vandføring (l/s)<br />

70<br />

60<br />

50<br />

40<br />

30<br />

20<br />

10<br />

1<br />

2<br />

3<br />

4<br />

5<br />

Målinger 27.02 og 28.02 Målinger 20.04<br />

0<br />

0 200 400 600 800 1000 1200<br />

Målingernes indbyrdes afstand (m)<br />

Figur 3.22. Vandføring ved målepunkterne i Hølbækken. Usikkerhedsintervallet er sat til 10 %, hvilket er<br />

angivet for at dække over de usikkerheder, der er forbundet med forsøget, bilag B.<br />

Vandføringsmålingerne viser, at Hølbækken bidrager med ca. 60 l/s til Lindenborg Å.<br />

Vandføringen fordobles fra øverste til nederste vandføringsmåling, figur 3.22, dvs. at<br />

næsten halvdelen af vandføringen stammer fra den øverste del af vandløbet. Der blev<br />

registreret flere kilder opstrøms for målepunkt 1, afsnit 3.1.2 (Observationer), der<br />

givetvis står for en del af vandføringen i måling 1.<br />

Analyse af resultaterne<br />

De første to målepunkter (måling 1 og 2) er taget på hver sin side af den største kilde,<br />

for at få et indtryk af afstrømningen fra denne. Forskellen mellem disse to målepunkter<br />

er på ca. 6 l, og da den samlede udstrømning til Lindenborg Å er på ca. 60 l, kommer<br />

tæt på 10 % af Hølbækkens vandføring fra denne kilde. Målepunkt 1 og 2 er ligger<br />

opstrøms gytjelaget. Mellem målepunkt 2 og 3 falder vandføringen med 2,4 l/s. Dette<br />

fald kunne teoretisk skyldes, at grundvandsspejlet ligger under vandspejlet i åen og<br />

vandet således strømmer væk fra åen og ud i den omkringliggende jord. Dette anses dog<br />

ikke for sandsynlig, da grundvandsspejlet i hele området var højtliggende, hvorfor det<br />

derfor ikke synes sandsynligt, at vandet på dette sted kan strømme ud af vandløbet.<br />

Gytjelagets udbredelse stopper mellem måling 2 og 3, hvilket kunne tænkes at have en<br />

indflydelse på nedgangen i vandføring mellem måling 2 og 3, idet gytjen som aquitard<br />

begrænser vandets strømningshastighed ud i Hølbækken. Målingerne i målepunkt 2 og<br />

3 blev taget to forskellige dage, hvilket også kunne være medvirkende til nedgangen i<br />

målt vandføring. Imidlertid indikerer vejrforholdene, bilag B, at vandføringen i givet<br />

6<br />

7<br />

8<br />

9<br />

49


Kapitel 3 Konceptuel model af projektområdet<br />

fald burde være større i målepunkt 3. Vandføringen falder ligeledes mellem målepunkt<br />

8 og 9. Forskellen mellem de to målinger er i begge tilfælde under 5 l/s og har<br />

overlappende usikkerhedsintervaller, hvorfor det antages, at faldet kan tilskrives<br />

måleusikkerhed.<br />

Mellem målepunkt 3 og 4 blev der lavet en vandføringsmåling i kilde E, der strømmede<br />

til Hølbækken umiddelbart opstrøms for målepunkt 4. Vandføringen blev her målt til<br />

1,78 l/s. Kildens vandløb var tilnærmelsesvis kvadratisk, hvorfor det var let at foretage<br />

målingen og udregningen. Imidlertid blev hastigheden kun målt midt i kilden, og da<br />

vandføringen er langsommere langs vandløbets kanter grundet friktion, må det forventes,<br />

at den målte vandføring er et overestimat. Målingen illustrerer, at tilstrømning i<br />

dette område derfor må kunne tilskrives diffus udstrømning, sandsynligvis fra<br />

brinkerne. Mellem målepunkt 4 og 5, som ligger forholdsvis tæt på hinanden, sker der<br />

en forøgelse af vandføringen på ca. 9 l/s. Som nævnt i afsnit 3.1.2 (Observationer) blev<br />

der registreret et sumpet område langs vandløbsbredden mellem disse to målinger, hvor<br />

der var mindst en kilde. Dette kunne tyde på, at der på dette sted er et brud i gytjelaget,<br />

der muliggør opstrømning af vand fra de omkringliggende lag. I afsnit 3.2.3<br />

(Håndboringer i fokusområdet) blev det beskrevet, at der ved dette målepunkt blev<br />

registreret et brud på gytjelaget i undergrunden på Hølbækkens vestlige bred, hvor der i<br />

stedet var et dybere tørvelag. Området var desuden sumpet, hvilket tyder på et<br />

højtliggende grundvandsspejl. I dette våde område blev der observeret kilder på begge<br />

sider af åen. Stigningen i vandføring kan derfor skyldes, at vandet her kan strømme<br />

lettere til bækken, både gennem vådområdet og kilderne i dette område.<br />

Mellem målepunkt 5 og 9 har Hølbækken en jævn stigning i vandføringen, og alle<br />

usikkerhedssintervaller er overlappende. Den jævne vandføring skyldes sandsynligvis,<br />

at det nederste område af ådalen mod Lindenborg Å er forholdsvis fladt, hvor<br />

grundvandsspejlets gradient er lille. Derudover er der heller ikke observeret kilder på<br />

den nederste del, men der findes en del dræn på den nedre del af Hølbækken, som det<br />

må antages ikke bidrager noget særligt til vandføringen i den baseflowsituation, der her<br />

er lavet undersøgelser på.<br />

Samlet set kan det siges, at stort set alt vandet, der strømmer i Hølbækken, stammer fra<br />

den øvre del, idet vandføringsmåling 6 udgør grænsen mellem de to dele. Dette<br />

indikerer igen, at gytjelaget, der er mest udbredt på den nedre del, er begrænsende for<br />

udstrømningen af grundvand til Hølbækken, da omkring halvdelen af udstrømningen fra<br />

ådalen til vandløbet er sket før gytjelaget. Derfor er især kilderne på den øvre del vigtige<br />

i forhold til at bidrage med vand til Hølbækken.<br />

50


Ændringer i Hølbækkens vandføring<br />

For at undersøge ændringerne i Hølbækkens vandføring over tid, blev trykniveauet i<br />

bækken og i atmosfæren målt over en periode på 52 dage mellem den 27.02 og den<br />

20.04. Beskrivelsen af forsøget kan ses i bilag B. Data findes på den vedlagte Cd-rom.<br />

Forskellen mellem de to trykniveauer, hhv. i Hølbækken og atmosfæren, viser relative<br />

ændringer i vandstanden og er vist på figur 3.23. Det ses, at Hølbækkens trykniveau<br />

ligger minimum mellem 5 og 10 cm, og ved nedbørshændelser øges vandføringen<br />

efterfølgende. I den sidste del af perioden viser divermålingerne, at trykniveauet lå 2-3<br />

cm højere end i starten af perioden, hvilket kan tyde på, at der er sket en stigning i<br />

tilstrømningen fra overfladen og fra grundvandet, på grund af nedbøren.<br />

Den 20.04 blev vandføringsmålinger ved målepunkt 6 og 7 gentaget, med henblik på, at<br />

have vandføringsmålinger fra 2 forskellige dage at tolke målingerne af vandstanden ud<br />

fra, figur 3.22. Det ses, at vandføringen ved begge målepunkter tilnærmelsesvis var den<br />

samme, som ved den første måling. Divermålingerne tyder derfor på at Hølbækken er<br />

grundvandsfødt, da der ikke er de store udsving i vandstanden. Den målte vandføring<br />

repræsenterer således baseflowet i Hølbækken for en vintersituation. Baseflow<br />

vandføring for Hølbækken er beregnet til 5.160,5 m 3 /d, hvis det målte baseflow ganges<br />

op. Dette er formodentlig mere end baseflow i sommermånederne, men en<br />

underestimering af den gennemsnitlige årlige vandføring pga. af at nedbørshændelser<br />

ikke er repræsenteret. Resultatet vil i den resterende del af rapporten repræsentere den<br />

gennemsnitlige vandføring, da det ikke har været muligt at finde mere præcise data for<br />

denne. Dette giver en vandføring på ca. 1,89 mio. m 3 /år.<br />

Vandtryk (cm)<br />

25<br />

20<br />

15<br />

10<br />

5<br />

0<br />

28.02 6.03 13.03 19.03 25.03 31.03 7.04 13.04 19.04<br />

Figur 3.23. Relative ændringer i Hølbækkens vandstand fundet ved forskellen mellem atmosfærens tryk<br />

og trykket i åen. Målingerne mellem 6.03 og 19.03 er ikke taget med, da disse målinger betragtes som<br />

fejlmålinger, bilag B.<br />

Dato<br />

51


Kapitel 3 Konceptuel model af projektområdet<br />

Afstrømning fra Lindenborg Å<br />

I forhold til formålet med denne rapport er det væsentlig at finde et estimat for<br />

medianminimumsvandføringen for Hølbækken. Den målte vandføring repræsenterer et<br />

øjebliksbillede af baseflowet i en vintersituation, og det må forventes, at dette er højere<br />

end medianminimumsvandføringen, idet fordampningen er lille om vinteren. Der er<br />

ikke lavet længerevarende tidsserier over vandføringen i Hølbækken, hvorfor det er<br />

nødvendigt at opskalere og undersøge et større område, for derefter at bruge de fundne<br />

værdier til at vurdere størrelsen af medianminimumsvandføringen i Hølbækken. Det er<br />

valgt at bruge vandføringsmålinger fra Lindenborg Å til dette formål, da<br />

projektområdets opland er et delopland til Lindenborg Å. Vandføringen i Lindenborg Å<br />

er målt ved målestationen Lindenborg Bro, der ligger ca. 4 km nedstrøms ad<br />

Lindenborg Å fra Hølbækkens udløb og oplandet hertil er på 218 km 2 (NJA, 2005),<br />

figur 3.24.<br />

Figur 3.24. Hydrologiske opland til vandløbsstationen Lindenborg Bro. Data fra (NJA, 2005) og<br />

TOP10DK (2000).<br />

Undergrunden i oplandet til denne målestation er blevet undersøgt ud fra GEUS<br />

boredata, og herudfra er det fundet, at undergrunden i Lindenborg Bros opland ligner<br />

undergrunden i projektområdet i opbygning. Over kalken ligger morænesand og –ler.<br />

Vandføringen ved Lindenborg bro kan ses på Figur 3.25.<br />

52


Vandføring (m 3 /s)<br />

7<br />

6<br />

5<br />

4<br />

3<br />

2<br />

1<br />

0<br />

j-03<br />

f-03<br />

m-03<br />

a-03<br />

m-03<br />

j-03<br />

j-03<br />

a-03<br />

s-03<br />

o-03<br />

n-03<br />

d-03<br />

j-04<br />

f-04<br />

Måned-år<br />

Figur 3.25. Hydrograf for Lindenborg Å, Lindenborg Bro (DMU 2, 2006).<br />

Lindenborg Å har ved Lindenborg Bro en gennemsnitlig vandføring på 2,5 m 3 /s, et<br />

baseflow på 2 m 3 /s og en medianminimumsvandføring på 1,7 m 3 /s. Dermed udgør<br />

baseflowet 80 % af den samlede vandføring, hvilket betyder at Lindenborg Å kan<br />

klassificeres som grundvandsfødt. Medianminimumsvandføringen udgør ca. 71 % af<br />

gennemsnitsvandføringen. På hydrografen ses, at vandføringen er lidt større om<br />

vinteren end om sommeren, hvilket kan tilskrives større nedbørsmængder i<br />

vinterhalvåret. Der ses også forskellige regnhændelser, hvor de største er i efteråret. Til<br />

tider kan der komme regnhændelser, der tredobler afstrømningsmængden, f.eks. februar<br />

2004. De store udsving viser, at der er god sammenhæng mellem det overfladenære<br />

grundvand og Lindenborg Å.<br />

Afstrømningen fra Hølbækken<br />

Da oplandet til Hølbækken er et delopland til Lindenborg Å og da de begge er<br />

grundvandsfødte med nogenlunde samme type undergrund, er der mulighed for at drage<br />

paralleller mellem dem. Hydrografen for Lindenborg Bro viste, at<br />

medianminimumsvandføringen her var på ca. 71 % af den gennemsnitlige vandføring.<br />

Det kan derfor forventes, at medianminimumsvandføringen for Hølbækken ligger i<br />

nogenlunde samme størrelsesorden. Det kan forventes, at medianminimumsvandføringen<br />

påvirkes mere af grundvandsindvindingen end gennemsnitsvandføringen. Det<br />

skyldes, at indvindingen af grundvand i Volsted Plantage vil ske fra det dybere<br />

grundvand og den dybe grundvandsafstrømning udgør en større andel af medianminimumsvandføringen<br />

end af gennemsnitsvandføringen. Beregningerne i denne<br />

m-04<br />

a-04<br />

m-04<br />

j-04<br />

j-04<br />

a-04<br />

s-04<br />

o-04<br />

n-04<br />

d-04<br />

j-05<br />

53


Kapitel 3 Konceptuel model af projektområdet<br />

rapport vil ske på baggrund af den formodede gennemsnitsvandføring, men det antages,<br />

at påvirkningen af medianminimumsvandføringen vil være højere.<br />

3.3.2 Fordampning og nedbør<br />

Hvor stor grundvandsdannelsen er i et område, afhænger af nedbøren, jordbundstypen,<br />

relieffet og evapotranspiration, der igen afhænger af vegetationen. For at have styr på<br />

grundvandsdannelsen i et område, er det nødvendigt både at kende nedbøren og<br />

fordampningen. Forskellen mellem disse to kaldes for nettonedbøren. Nettonedbøren<br />

kan både infiltrere ned til grundvandet eller løbe af på jordoverfladen. Det er forsøgt at<br />

udregne nettonedbøren ud fra vandløbsstationen Lindenborg Bro, men det er vurderet,<br />

at dette resultat er for upræcist, da der er for store usikkerheder i bestemmelsen af<br />

arealet af oplandet og oppumpningen heri. Da Hølbækken er grundvandsfødt, antages<br />

det, at størstedelen af nettonedbøren infiltrerer ned til grundvandet og at den mængde<br />

der afstrømmer fra Hølbækken repræsenteres af infiltrationen i oplandet. Den<br />

nettonedbør der her angives må derfor være i underkanten af det den i virkeligheden er,<br />

men betragtes som det bedste bud. Fordampningen kan derfor med rimelig nøjagtighed<br />

bestemmes som forskellen mellem nedbør og infiltration.<br />

Den geografiske fordeling af nedbør og infiltrationen i projektområdet ses på figur<br />

3.26. Det ses, at der falder ca. 725 mm/år i projektområdet. Der falder lidt mere i<br />

sydvest og lidt mindre i nordøst. Infiltrationen i projektområdet er et estimat udarbejdet<br />

af NJA (2003), der bygger på typetal for vegetation og jordbund i området. Typetallene<br />

giver et estimat over infiltrationen baseret på standardværdier, og repræsentationen i<br />

kortet er i 500m * 500 m, hvilket giver visse unøjagtigheder i resultatet. I gennemsnit er<br />

infiltrationen i projektområdet på ca. 308 mm/år. Forskellen mellem nedbøren og<br />

infiltrationen er på ca. 400 mm/år, der udgør fordampningen i området.<br />

3.3.3 Oppumpning og afledning<br />

Placering af kildepladser for vandindvinding i eller tæt på projektområdet ses på figur<br />

3.27. Det ses, at inden for projektområdet findes vandindvinding i Volsted og<br />

Essendrup. Udenfor området ligger Oppelstrup og Gunderup, som begge er forbundet<br />

med Essendrup. Da Essendrup er forbundet med indvindinger udenfor området, kan der<br />

både føres vand frem og tilbage mellem de forskellige boringer, hvilket gør det<br />

vanskeligt at holde styr på vandstrømmene.<br />

54


Figur 3.26. Tv. Nedbørens geografiske fordeling. Klimanormal for 1961 til 1990 (DJF-geodata, 2006).<br />

Th. Infiltrationen fordelt over projektområdet (NJA, 2003).<br />

Figur 3.27. Placering af vandindvindinger i og tæt på projektområdet. Data fra TOP10DK (2000) og AK<br />

(2005).<br />

55


Kapitel 3 Konceptuel model af projektområdet<br />

På figur 3.28 ses de indvundne vandmængder for de fire vandværker. Det ses, at<br />

Essendrups indvinding er lille i forhold til de andre vandværker. Samtidig forbruger<br />

Essendrup ca. 2.000 m 3 /år (AK, 2005) og dermed bidrager denne kildeplads ikke<br />

væsentligt til import eller eksport af vand. Volsted indvinder mellem 25.000 og 15.000<br />

m 3 /år, hvor den laveste er fra 2002. Generelt har der været en faldende tendens i<br />

vandforbruget siden midten af 1980’erne. I forhold til Hølbækkens afstrømning på ca.<br />

1,89 mio. m 3 /år, udgør oppumpningen en ganske lille del af den samlede grundvandsdannelse<br />

i projektområdet.<br />

56<br />

Indvinding 1.000 m 3<br />

45<br />

40<br />

35<br />

30<br />

25<br />

20<br />

15<br />

10<br />

5<br />

0<br />

1985 1997 2002<br />

Årstal<br />

Figur 3.28. Vandindvindingen ved projektområdet (AK, 2006).<br />

Essendrup<br />

Gunderup<br />

Oppelstrup<br />

Volsted<br />

Det oppumpede vand vil være af størst betydning for vandbalancen, hvis det<br />

transporteres ud af projektområdet. Bebyggelsen i projektområdet er hovedsageligt<br />

spredte gårde, der formentlig udleder deres vand i nærheden af gården. Volsted er det<br />

eneste område, hvor der samles spildevand til rensning. Rensningen foregår på<br />

Rensningsanlæg Vest i Aalborg (AK, 2006). Spildevandet forlader altså projektområdet.<br />

Derfor antages det, at alt vandet, der oppumpes ved Volsted, forlader projektområdet,<br />

mens den øvrige oppumpning bliver returneret som spildevand. Da nedgangen i<br />

vandforbruget ikke er entydig antages det, at oppumpningen af vand i projektområdet,<br />

svarer til det nyeste tal for oppumpning i Voldsted boringen dvs. ca. 15.000 m 3 /år.<br />

3.4 Projektområdets hydrogeologi<br />

I dette afsnit vil projektområdets hydrogeologi blive undersøgt. Dette gøres for at<br />

bestemme nogle af de faktorer, der har indflydelse på grundvandsstrømningerne i<br />

projektområdet, der skal anvendes til opsætningen af den numeriske grundvandsmodel.


I forbindelse med feltturen blev der gennemført en række undersøgelser i ådalen, figur<br />

3.29. På den første felttur blev der udtaget intaktprøver af de forskellige lag i<br />

Hølbækkens ådal. Disse blev indsamlet for at bestemme lagenes hydrauliske<br />

ledningsevne for de lag, der i den geologiske model blev tolket til at være af indflydelse.<br />

Dette er et væsentligt input til grundvandsmodellen. Efter prøvernes hydrauliske<br />

ledningsevner var bestemt ved gennemstrømningsforsøg, blev der foretaget<br />

fasebestemmelse af prøverne. For at kunne vurdere sammenhængen mellem prøverne<br />

og de hydrauliske ledningsevner, vil resultaterne fra fasebestemmelserne først blive<br />

gennemgået, efterfulgt af resultaterne fra gennemstrømningsforsøget. Herefter vil<br />

resultaterne fra udførte slugtests, som også giver et bud på den hydrauliske<br />

ledningsevne, blive præsenteret. For at bestemme interaktionen mellem vandløbet og<br />

ådalsmagasinet er der lavet et forsøg med måling af åbundens hydrauliske ledningsevne<br />

og et forsøg med at bestemme den diffuse udstrømning gennem åbunden.<br />

Figur 3.29. Lokaliteterne hvor der er udført undersøgelser og udtaget prøver i afsnittene. Åbunden<br />

henviser til de steder, hvor der blev lavet bestemmelse af åbundens hydrauliske ledningsevne.<br />

3.4.1 Fasebestemmelse<br />

Under feltturen blev der udtaget 20 intaktprøver fra de 4 aflejringer, der udgør<br />

opfyldningen i ådalen, fem prøver fra hvert af lagene: tørv, sand samt øverste og<br />

nederste del af gytjelaget, tabel 3.1. Tre prøver fra hvert lag blev analyseret. Prøverne<br />

57


Kapitel 3 Konceptuel model af projektområdet<br />

blev alle udtaget i den øvre del af Hølbækkens ådal, figur 3.29. Der blev udtaget to<br />

prøver fra gytjen, da det i håndboreundersøgelsen var observeret, at der skete en<br />

ændring i gytjelagets sammensætning. Prøverne blev undersøgt for vandindhold,<br />

porøsitet og indhold af organisk materiale. Endvidere er det undersøgt, hvor stor en del<br />

af materialet i aflejringerne, der er under 0,125 mm, dvs. mindre end finsand (Finneren<br />

et al., 1996: 132). Resultaterne kan ses i tabel 3.2. Udførlig beskrivelse af forsøgene<br />

kan ses i bilag C.<br />

58<br />

Tabel 3.1. Aflejringstype og dybde for udtagelse.<br />

Prøve Dybde (cm)<br />

Tørv 75-90<br />

Gytje (øverst) 180-240<br />

Gytje (nederst) 240-280<br />

Sand 60-75<br />

Tabel 3.2. Uddrag af resultaterne fra fasebestemmelsen. Tallene er gennemsnit fra de tre prøver for hver<br />

lokalitet, da afvigelser mellem prøverne var ubetydelige. Det var ikke muligt at finde volumen for hhv.<br />

uorganisk og organisk materiale.<br />

Tørv<br />

Øvre<br />

gytje<br />

Vægtprocent<br />

(gennemsnit)<br />

Vand 90,5 %<br />

Porøsitet<br />

(min. – maks.)<br />

Andel af materiale<br />

< 0,125 mm<br />

Uorganisk materiale 1 %<br />

88,5-90 %<br />

Ingen data<br />

Organisk materiale 8,5 %<br />

Vand 78,3 %<br />

Uorganisk materiale 16,3 % 42,6 %<br />

85-87 %<br />

Organisk materiale 5,3 %<br />

13,6 %<br />

Vand 73 %<br />

Nedre<br />

Uorganisk materiale 22 %<br />

81-86 %<br />

42,6 %<br />

gytje Organisk materiale 5 %<br />

11,2 %<br />

Sand<br />

Vand 29,5 %<br />

Uorganisk materiale 66,5 % 7,4 %<br />

51-53 %<br />

Organisk materiale 4 %<br />

1,1 %<br />

Tabel 3.3 viser tal for porøsitet og effektiv porøsitet, fundet af Nilsson et al. (2003).<br />

Den effektive porøsitet er den del af det totale porevolumen, hvori vandgennemstrømningen<br />

foregår (Nilsson et al., 2003: 63). Tallene anvendes som sammenligningsgrundlag<br />

for resultaterne fra fasebestemmelsen samt til at finde et udtryk for den<br />

effektive porøsitet for sedimentprøverne.


Tabel 3.3. Tabellen viser forskellige værdier for porøsiteten og den hydrauliske ledningsevne (Nilsson et<br />

al., 2003:63).<br />

Materiale Porøsitet % Effektiv porøsitet %<br />

Grovkornet sand 25 – 30 25 – 30<br />

Mellemkornet sand 30 – 40 20 – 25<br />

Finkornet sand 40 – 50 15 – 20<br />

Silt 35 – 50 10 – 15<br />

Ler 40 – 70 5 – 10<br />

Gytje 80 – 90 20 – 30<br />

Svagt humificeret tørv 90 – 98 40 – 90<br />

Moderat humificeret tørv 85 – 90 15 – 40<br />

Stærkt humificeret tørv 60 – 85 10 – 15<br />

Tørv: Tørven havde et højt vandindhold og bestod derudover hovedsagelig af<br />

organisk materiale (over 90 % af den samlede masse). Porøsiteten blev fundet til<br />

at ligge inden for intervallet svagt til moderat humificeret tørv, tabel 3.3. Ud fra<br />

tabellen fremgår det, at strømningen foregår i 15 – 90 % af porerne i tørven.<br />

Tørvprøverne blev ikke sigtet, da det var umuligt at opdele materialet. Tydelige<br />

fibre på flere centimeters længde blev observeret i prøverne. Dette tyder på, at<br />

det organiske materiale er dårligt omsat, hvilket indikerer, at der er tale om svagt<br />

humificeret tørv.<br />

Øvre gytje: Porøsiteten for gytjen lå i det forventede interval i forhold til tabel<br />

3.3. Dette svarer til en effektiv porøsitet på 20 – 30 %, hvilket viser, at<br />

hovedparten af porerne i gytjen er så små, at vandet bliver bundet til partiklerne<br />

og derfor giver en lav samlet hydraulisk ledningsevne. Over 50 % af det<br />

organiske materiale er < 0,125 mm, hvilket bekræfter, at det er et organisk<br />

sediment og giver noget af forklaringen på, hvorfor den effektive porøsitet<br />

formodentlig er kun er 20 – 30 % af den samlede porøsitet. Under forsøget blev<br />

det observeret, at det fandtes større fibre af organiske materiale i den øvre gytje<br />

end i den nedre.<br />

Nedre gytje: Det organiske indhold i den nedre gytje er det samme som i den<br />

øvre gytje, men den nedre gytje har en lavere porøsitet, og forventes derfor at<br />

have en lavere effektiv porøsitet. Grunden til denne forskel er efter al<br />

sandsynlighed, at den nedre gytje er mere kompakteret, hvorfor porerummene er<br />

presset mere sammen. Andelen af organisk materiale < 0,125 mm er desuden en<br />

anelse større for den nedre gytje end for den øvre, hvilket vidner om, at det<br />

organiske materiale er mere omsat i det nedre lag.<br />

59


Kapitel 3 Konceptuel model af projektområdet<br />

60<br />

Sand: Indholdet af organisk materiale er i sandet betragtelig mindre end i gytjen<br />

og tørven, hvilket bekræfter, at jorden i højere grad består af minerogene<br />

sedimenter. Porøsiteten svarer til finkornet sand, hvilket må betyde, at sandets<br />

kornstørrelser findes i den finere ende af intervallet 0,125 – 2 mm.<br />

3.4.2 Hydraulisk ledningsevne<br />

I det følgende vil resultaterne af gennemstrømningsforsøget blive præsenteret og<br />

analyseret. Undersøgelsen blev udført på de samme prøver som der blev lavet<br />

fasebestemmelse på. Resultaterne vil blive sat i relation til fasebestemmelsen af<br />

prøverne. Resultaterne fra gennemstrømningsforsøget kan ses i tabel 3.4. Generelt ses<br />

det, at der ikke er så stor forskel mellem prøver fra samme lag, men at der er ret stor<br />

forskel lagene imellem. Den højeste hydrauliske ledningsevne er fundet i tørven, mens<br />

den laveste er fundet i det nederste gytjelag. I tabel 3.5 ses værdier for hydraulisk<br />

ledningsevne for forskellige sedimenttyper, fundet af Nilsson et al. (2003).<br />

Tabel 3.4. Hydraulisk ledningsevne for de materialer der er blevet lavet forsøg på.<br />

Materiale Forsøg Hydraulisk ledningsevne<br />

(cm/s)<br />

Rør 1 2,62 x 10 -4<br />

Rør 2 4,35 x 10 -4<br />

Tørv<br />

Rør 3 5,35 x 10 -4<br />

Øvre gytje<br />

Nedre<br />

gytje<br />

Sand<br />

Rør 1 5,29 x 10 -6<br />

Rør 2 6,64 x 10 -6<br />

Rør 3 Ingen data<br />

Rør 1 1,42 x 10 -6<br />

Rør 2 1,19 x 10 -6<br />

Rør 3 0,48 x 10 -6<br />

Rør 1 2,36 x 10 -5<br />

Rør 2 1,43 x 10 -5<br />

Rør 3 Ingen data<br />

Gennemsnit<br />

(cm/s)<br />

4,11 x 10 -4<br />

5,97 x 10 -6<br />

1,03 x 10 -6<br />

1,89 x 10 -5<br />

Tørven har en hydraulisk ledningsevne, der svarer til svagt til moderat<br />

humificeret tørv, tabel 3.5, hvilket understøttes af den fundne porøsitet for<br />

tørven, afsnit 3.4.1 (Fasebestemmelse). I dette lag er der derfor god overensstemmelse<br />

mellem tabelværdi og resultat fra gennemstrømningsforsøg.<br />

Gytjen ligger for begge de gennemførte forsøg inden for intervallet af, hvad der<br />

kan forventes for et gytjelag. Intervallet er stort og den hydrauliske ledningsevne<br />

ligger for begge lag i den lave ende af skalaen. I forhold til tørvens hydrauliske


ledningsevne understreger dette, at gytjens begrænsende effekt på<br />

grundvandsstrømningerne i ådalen. Det ses af resultaterne, at den hydrauliske<br />

ledningsevne for det øverste gytjelag var fem gange højere end for det nederste<br />

gytjelag. En forklaring på denne forskel er, at det nederste lag er blevet mere<br />

kompakteret. De to prøver havde forholdsvis ens kornstørrelsesfordeling, men<br />

porøsiteten var lavere i den nedre gytje, afsnit 3.4.1 (Fasebestemmelse). I det<br />

nederste gytjelag er der en faktor tre i forskel mellem prøven med den højeste og<br />

laveste hydrauliske ledningsevne. Det var ikke forventeligt, at der skulle være så<br />

stor forskel mellem disse prøver, da de er indhentet fra samme lokalitet. Det blev<br />

fundet gennem fasebestemmelsen, at denne forskel skyldes forskel i porøsitet,<br />

rør 1 har højere porøsitet end rør 3, lige som rør 3 har et større indhold af små<br />

partikler, der kan begrænse vandets strømning gennem jorden.<br />

Sandet har en lavere hydraulisk ledningsevne end det var forventet for en sandet<br />

prøve. Dette hænger sandsynligvis sammen med det høje indhold af små<br />

partikler, der er blevet observeret under fasebestemmelsen.<br />

Tabel 3.5. Intervaller for hydraulisk ledningsevne i forhold til sedimenttype (Nielsen et al., 2004: 63).<br />

Materiale Hydraulisk<br />

ledningsevne (cm/s)<br />

Grovkornet sand 10 -2 – 10 -1<br />

Mellemkornet sand 10 -3 – 10 -2<br />

Finkornet sand 10 -5 – 10 -3<br />

Silt 10 -8 – 10 -4<br />

Ler 10 -11 – 10 -8<br />

Gytje 10 -8 – 10 -5<br />

Svagt humificeret tørv 10 -4 – 10 -1<br />

Moderat humificeret tørv 10 -6 – 10 -4<br />

Stærkt humificeret tørv 10 -7 – 10 -8<br />

I ovenstående diskussion af resultaterne er det kommet frem, at der for alle de<br />

undersøgte prøver var god overensstemmelse mellem den fundne hydrauliske<br />

ledningsevne i gennemstrømningsforsøget og den standardværdi, der er angivet i<br />

litteraturen. Derudover er det kommet frem, at der var en forskel i den hydrauliske<br />

ledningsevne, der blev bestemt for den øverste og den nederste del af gytjen og at denne<br />

forskel muligvis kan forklares ud fra, at den nederste del af laget er mere kompakteret<br />

end den øverste del.<br />

61


Kapitel 3 Konceptuel model af projektområdet<br />

3.4.3 Slugtest<br />

Målet med at gennemføre slugtesten var at bestemme et direkte udtryk for den samlede<br />

hydrauliske ledningsevne i de lag, der findes i fokusområdet. Dette vil blive<br />

sammenholdt med resultaterne fra undersøgelsen af den hydrauliske ledningsevne i<br />

sidste afsnit. Slugtest udføres ved at hæve eller sænke grundvandsspejlet i et borehul og<br />

måle højden på grundvandsspejlet over tid, indtil det er tilbage i udgangssituationen.<br />

Fremgangsmåde for udførelse af forsøget er beskrevet i bilag E. Placeringen af<br />

borehuller, hvor der er udført slugtests findes på figur 3.29, resultaterne fra de udførte<br />

slugtests kan ses i tabel 3.6.<br />

Tabel 3.6. Resultater fra de udførte slugtests.<br />

62<br />

Konduktivitet ved<br />

nedsænkning af slug (cm/s)<br />

Konduktivitet ved optagning<br />

af slug (cm/s)<br />

Hul 1 Hul 2 Hul 3 Hul 4<br />

2,00 x 10 -3<br />

0,65 x 10 -3<br />

0,71 x 10 -3<br />

0,69 x 10 -3<br />

Materiale Tørv Tørv<br />

1,50 x 10 -3 0,79 x 10 -3<br />

1,20 x 10 -3 0,28 x 10 -3<br />

Tørv og<br />

gytje<br />

Sand<br />

Slughul 1<br />

I dette hul blev der boret gennem tørv i de øverste 160 cm og da hullet til slugtesten er<br />

135 cm dybt, er der hovedsageligt blevet målt hydraulisk ledningsevne på tørv. Den<br />

hydrauliske ledningsevne var ved hævning af grundvandsspejlet (nedsænkning af<br />

sluggen) en faktor 3 større end ved sænkning af grundvandsspejlet (optagning af<br />

sluggen), hhv. 2,00 x 10 -3 cm/s og 0,65 x 10 -3 cm/s, tabel 3.6. Denne forskel er<br />

bemærkelsesværdig, da det var forventet, at de to delforsøg ville give tilnærmet samme<br />

hydrauliske ledningsevne. En mulig begrundelse for denne afvigelse er, at der var<br />

forskel i hullets gennemsnitlige radius over og under vandspejlet. Dette kommer af, at<br />

hullet faldt noget sammen under grundvandsspejlet. Dette resulterede så i, at<br />

højdeforskellen mellem vandspejlet og udgangssituationen var større, når sluggen blev<br />

lagt i hullet end når den blev taget op. Grunden til dette var, at radius i den øverste del<br />

af hullet var mindre end i den nedre del. I princippet burde der ikke være forskel mellem<br />

den hydrauliske ledningsevne der beregnes, da der i beregningen er taget højde for dette<br />

forhold (der regnes ikke på den absolutte stigning og det absolutte fald, men på den tid<br />

det tager for 37 % af vandet, at strømme hhv. ind og ud af hullet). I beregningen bliver<br />

der imidlertid ikke taget forbehold for, at hullets diameter i forsøget varierede. Betragtes<br />

h/h0 – kurven for nedsænkningen af sluggen ses det, at den buer nedad, bilag E. h/h0kurven<br />

angiver forholdet mellem den højde vandet steg eller faldt til og udgangshøjden<br />

over tid. Har hullet den samme radius hele vejen ned gennem vandsøjlen, skal h/h0-


kurven være ret, dette selvfølgelig forudsat, at materialet ikke ændrer karakter ned<br />

gennem hullet. Når h/h0-kurven derfor buer nedad, skyldes det, at hullets radius er<br />

mindre over grundvandsspejlet, idet der ikke var de store ændringer i materialet ned<br />

gennem vandsøjlen. Problemet er imidlertid ikke det samme når sluggen tages op af<br />

hullet, hvilket ses af den tilnærmelsesvis rette h/h0-kurve, bilag E. Resultatet fra<br />

sænkningen af grundvandsspejlet vurderes derfor til at være brugbare. Den hydrauliske<br />

ledningsevne for tørv fundet i gennemstrømningsforsøget var en faktor 10 lavere end<br />

slugtestens resultater. Denne forskel er vigtig at huske på når grundvandsmodellen skal<br />

kalibreres. Forskellen mellem de to resultater kan skyldes, at de måles på forskellige<br />

skalaer (Jensen 2002: 33-34).<br />

Slughul 2<br />

I dette hul blev der udelukkende fundet tørv. Den hydrauliske ledningsevne for dette<br />

slughul kan derfor betragtes som et udtryk for den hydrauliske ledningsevne for tørv.<br />

Resultatet fra slughul 2 viser, modsat slughul 1, at der er god overensstemmelse mellem<br />

værdierne for hævning og sænkning af grundvandsspejlet, hhv. 0,71 x 10 -3 cm/s og 0,69<br />

x 10 -3 cm/s. Derudover er der god overensstemmelse mellem den hydrauliske<br />

ledningsevne for tørven, bestemt i dette forsøg og i gennemstrømningsforsøget, hvor<br />

gennemsnittet er 0,411 x 10 -3 cm/s. h/h0-kurven er rimelig retliniet for både optagningen<br />

og nedsænkningen af sluggen og der er ikke så stor forskel på hullets radius over og<br />

under grundvandsspejlet, hvorfor resultaterne fra denne test er anvendelige.<br />

Slughul 3<br />

Resultatet af den tredje slugtest viser, ligesom slughul 2, små forskelle mellem den<br />

beregnede hydrauliske ledningsevne for nedsænkning og optagning af sluggen, hhv.<br />

1,50 x 10 -3 cm/s og 1,20 x 10 -3 cm/s. De geologiske lag i dette hul består i de øverste<br />

160 cm af tørv, mens det herunder findes gytje. Det er derfor umiddelbart overraskende,<br />

at den hydrauliske ledningsevne, der blev bestemt for dette hul er højere end for slughul<br />

2, da gytje er mindre permeabel end tørv. En af grundene til denne forskel, kan være, at<br />

hullet var boret noget dybere end den dybde diveren sank ned på. Hullet var boret<br />

omkring 2 meter dybere end den dybde, der blev anvendt i beregningerne, og dette kan<br />

forklare noget af forskellen. Det er imidlertid ikke hele forklaringen, idet<br />

genberegninger med andre dybder ikke alene kan forklare disse afvigelser. Det må<br />

derfor konkluderes, at de hydrauliske ledningsevner beregnet i forsøget dækker over<br />

lokale variationer.<br />

Slughul 4<br />

Ligesom i slughul 1 blev der i slughul 4 beregnet forskellig hydraulisk ledningsevne ved<br />

nedsænkning og optagning af sluggen, hhv. 0,79 x 10 -3 cm/s og 0,28 x 10 -3 cm/s. Dette<br />

resultat er omkring 20 gange højere, end den hydrauliske ledningsevne sandlaget blev<br />

63


Kapitel 3 Konceptuel model af projektområdet<br />

bestemt til i gennemstrømningsforsøget. En del af forklaringen på dette kan være, at der<br />

kunne observeres to sandlag i hullet, et mere leret end det andet og at<br />

gennemstrømningsforsøget er blevet udført på det mest lerede lag, mens slugtesten er et<br />

gennemsnit over den hydrauliske ledningsevne for begge lag. Ligesom i slughul 1, var<br />

der også her en forskel i radius af det borede hul over og under vandspejlet, hvilket kan<br />

være årsag til denne forskel. Det blev observeret, at der skete nogle udskridninger af<br />

materialet under boringen af hullet, hvilket kan forklare forskellen.<br />

4.4.4 Bestemmelse af åbundens hydrauliske ledningsevne<br />

I det følgende vil det blive undersøgt, hvordan kontakten mellem vandløbsbunden og<br />

ådalsmagasinet er. Dette er en vigtig faktor at kende i forhold til grundvandsmodelleringen,<br />

da åbundens hydrauliske ledningsevne siger noget om hvor let det er for<br />

vandet at strømme fra ådalsmagasinet og ud i vandløbet eller modsat. Åbundens<br />

hydrauliske ledningsevne blev bestemt ved at presse et ca. 1 m gennemsigtigt plasticrør<br />

ned gennem åbunden, hæve vandspejlet og måle den tid, det tog for det naturlige<br />

vandspejl at indstille sig, bilag F. Grunden til at der laves et særskilt forsøg til at<br />

bestemme åbundens hydrauliske ledningsevne, når den hydrauliske ledningsevne for de<br />

geologiske lag i ådalen allerede er bestemt via slugtesten og gennemstrømningsforsøget<br />

er, at der i åen sker forskellige processer, der kan ændre den hydrauliske ledningsevne<br />

for de geologiske lag. Eksempelvis kan der i bunden af vandløbet sedimenteres små<br />

partikler på strækninger, hvor vandløbet er aggraderende, eller vandløbet kan være<br />

degraderende og have skåret vandløbsbunden ned i en anden type sedimenter, der har en<br />

anden hydraulisk ledningsevne. Resultaterne kan ses i tabel 3.7.<br />

Åbundens hydrauliske ledningsevne blev bestemt to steder i vandløbet, figur 3.29. De<br />

øverste to meter af ådalen består af tørv, hvorfor dette materiale også muligvis også<br />

udgør bunden de steder, hvor undersøgelsen er lavet.<br />

64<br />

Tabel 3.7. Resultater fra bestemmelsen af åbundens hydrauliske ledningsevne.<br />

1. forsøg 2. forsøg<br />

Dybde af åen (cm) 25 35<br />

Rørets dybde i bunden (cm) 17,5 42<br />

Vandspejlets overhøjde ved start (cm) 42 18<br />

Hydraulisk ledningsevne (cm/s) 1,4 x 10 -4<br />

5,8 x 10 -4<br />

Den hydrauliske ledningsevne for åbunden er fire gange mindre i 1. forsøg end i 2.<br />

forsøg, tabel 3.7. En forklaring kan være, at tørvelaget har forskellig hydraulisk<br />

ledningsevne, hvilket også kunne ses af slugtesten. En anden forklaring er, at


vandspejlets overhøjde ved start var større i det første forsøg, hvilket har presset<br />

sedimenterne en smule sammen og dermed gjort den hydrauliske ledningsevne lavere.<br />

3.4.5 Måling af diffus udstrømning fra bunden<br />

Formålet med at måle udstrømningen fra bunden var at undersøge, hvor stor en del af<br />

åens vandføring, der stammer fra den diffuse udstrømning fra bunden og hvor meget der<br />

stammer fra andre strømningsvarianter. Forsøget blev udført ved hjælp af et<br />

seepagemeter. Forsøget blev udført ved at presse seepagemetret ned i åbunden og<br />

opsamle det vand der strømmede ud i en pose påmonteret i den ende røret er lukket,<br />

bilag G.<br />

Målingen blev foretaget mellem vandføringsmåling 5 og 6, hvor det kunne konstateres<br />

en stigning i vandføringen og hvor der kun var en mindre kilde. Stedet blev endvidere<br />

valgt, fordi ådalen her er ret smal og det derfor var forventeligt, at der var en større<br />

grundvandsgradient her, end på de fladere områder. Denne antagelse blev understøttet<br />

af, at der på dette sted i ådalen er et sumpet område, dvs. at grundvandet her står en<br />

anelse højere end i resten af ådalen. Det blev derfor forventet, at den diffuse<br />

udstrømning fra åbunden her ville være højere end den nedre del af vandløbet. Det blev<br />

forventet, at den diffuse udstrømning ville være lille, da den hydrauliske ledningsevne<br />

som tidligere nævnt er lav for opfyldningen i ådalen og da grundvandsgradienten i<br />

ådalen er lav.<br />

For at kunne bestemme den samlede diffuse udstrømning fra bunden af vandløbet til<br />

Hølbækken, blev bundens areal beregnet, arealet blev her beregnet til 2464 m 2, på<br />

baggrund af målinger af åbundens brede og længden af vandløbet, for beregning se<br />

bilag G.<br />

Resultat<br />

Antagelsen om, at den diffuse udstrømning fra vandløbsbunden til vandløbet ville være<br />

ubetydelig, viste sig at holde stik, idet den samlede udstrømning fra hele vandløbsbunden<br />

blev beregnet til 0,026 l/s, tabel 3.8. I forhold til den samlede vandføring i<br />

vandløbet på omkring 60 l/s, er dette meget lavt og må siges at være uden betydning.<br />

Tabel 3.8. Målinger af den diffuse udstrømning.<br />

Areal af cylinderen (cm 2 ) 1193<br />

Strømning (l/s) 1,28 x 10 -6<br />

Diffus udstrømning fra bunden (l/m 2 x s -1 ) 1,07 x 10 -5<br />

Samlet udstrømning fra bunden (l/s) 0,026<br />

65


Kapitel 3 Konceptuel model af projektområdet<br />

Ifølge ådalstypologien vil denne udstrømning stamme fra en kombination af den<br />

overfladenære diffuse udstrømning og den dybe regionale grundvandafstrømning. Af<br />

denne undersøgelse kan det derfor konkluderes, at kontakten mellem det regionale<br />

grundvandsmagasin og vandløbet er dårlig, da udstrømningen ville være højere, hvis<br />

kontakten havde været god. Dette stemmer godt overens med de boreprofiler, der er<br />

lavet i ådalen, som viser, at der i området findes et tykt gytjelag, der har en lav<br />

hydraulisk ledningsevne og begrunder derfor antagelsen om, at gytjelaget fungerer som<br />

aquitard. Det stemmer også godt overens med vandføringsmålingerne, hvor det kommer<br />

frem, at hovedparten af Hølbækkens vandføring kommer fra kilder ovenfor<br />

afgrænsningen af gytjelaget. Den overfladenære diffuse udstrømning strømmer i højere<br />

grad ud gennem brinkerne i siden af vandløbet, jf. kapitel 2 (Grundvand- og<br />

overfladevandsinteraktion). Det er ud fra denne undersøgelse derfor ikke til at<br />

konkludere, hvor stor den samlede diffuse udstrømning til vandløbet er.<br />

Det er tidligere blevet nævnt, at der i Hølbækkens ådal er mulighed for, at der vil være<br />

en komponent af underflow i vandløbet, da slyngningsgraden og vandløbets gradient,<br />

giver mulighed for dette. I dette afsnit er den diffuse udstrømning fra bunden blevet<br />

undersøgt, og det kom her frem, at den samlede udstrømning fra bunden var lav, samt at<br />

den hydrauliske ledningsevne der blev målt for både åbunden og ådalens opfyldning var<br />

lav. Af denne grund skønnes det, at der ikke er en betydelig komponent af underflow i<br />

Hølbækken. Dette understøttes endvidere af, at der i den nedre del af ådalen ligger et<br />

flere meter tykt gytjelag med en lav hydraulisk ledningsevne, der giver meget dårlige<br />

forhold for strømningen af vand i ådalen.<br />

3.5 Opsummering<br />

I dette kapitel er de undersøgelser, der er blevet lavet af projektområdet, blevet<br />

præsenteret og den konceptuelle forståelse af projektområdet, der vil blive arbejdet ud<br />

fra i forbindelse med opsætningen af grundvandsmodellen, er blevet præsenteret.<br />

Der er blevet udarbejdet en geologisk model for projektområdet. Geologien er en vigtig<br />

faktor i forbindelse med at forstå, hvordan grundvandet strømmer i projektområdet. Det<br />

er blevet gennemgået, at der i hele projektområdet findes et kalklag, som er det dybest<br />

liggende lag, hvor grundvandet strømmer. De fleste steder i projektområdet er kalklaget<br />

overlejret af morænesand, med forekomster af lerlinser. Nogle steder er kalklaget<br />

imidlertid sunket sammen eller er blevet borteroderet og der er dannet en ådal, ved<br />

delvis opfyldning af den eksisterende dal. I Hølbækkens og Lindenborgs ådal findes et<br />

lavpermeabelt gytjelag, der er overlejret af et ca. to meter tykt tørvelag. På den øvre del<br />

af Hølbækkens ådal bliver gytjelaget gradvist tyndere, for til sidst at forsvinde. Dette<br />

66


gytjelag har stor indflydelse på vandstrømningerne i ådalen, idet det fungerer som<br />

aquitard og begrænser kontakten mellem det dybe regionale grundvandsmagasin og<br />

Hølbækken. Både placeringen af kilderne og vandføringsmålingerne bekræfter, at<br />

udstrømningen fra grundvandet til Hølbækken hovedsageligt sker udenfor gytjelagets<br />

udbredelse. Imidlertid er der indlagt en række dræn i Hølbækkens ådal, der sikrer en<br />

hurtig afstrømning af vand. Dette betyder, at gytjelagets begrænsning af kontakten<br />

mellem grundvandet og Hølbækken modvirkes, således, at kontakten bliver bedre.<br />

Gennem dette kapitel er der fundet nogle af de parametre, der er input til<br />

grundvandsmodellen, kapitel 5 (Opsætning af grundvandsmodellen). Projektområdets<br />

nettonedbør er fundet ud fra infiltrationsdata. Hølbækkens afstrømning til Lindenborg Å<br />

er endvidere målt over 2 måneder, hvorudfra det er vurderet, at resultatet fra<br />

vandføringsmålingen kan repræsentere en gennemsnitsvandføring og et overestimat af<br />

medianminimumsvandføringen. Endelig er den hydrauliske ledningsevne for lagene i<br />

ådalen fundet. Dette er rapportens fokusområde, og er således grundigt blevet undersøgt<br />

i denne henseende. Parametre for kalken har det ikke været muligt at lave undersøgelser<br />

på, hvorfor det har været nødvendigt, at basere inputparametrene for grundvandsmodellen<br />

på sekundære kilder. Hydraulisk ledningsevne for morænesand og ler er<br />

ligeledes fundet i litteraturen, da det blev vurderet, ud fra rapportens fokus på<br />

Hølbækken og grundvandsspejlets beliggenhed i kalken, at det ikke var nødvendigt med<br />

så nøjagtige data herfra.<br />

Tabel 3.9. Inputparametre til grundvandsmodellen. Maksimum og minimum værdier for tallene kommer<br />

fra Nielsen et al. (2004) og Bach et al. (2001). Startværdier for ler og morænesand er midtpunktet af<br />

intervallet mellem maks. og min., mens startværdien for kalk er anslået og resten af startværdierne er<br />

taget fra afsnit 3.4.2 (Hydraulisk ledningsevne).<br />

Input Fastsat værdi<br />

Nettonedbør (m/dag) 0,00087<br />

Afstrømning (m 3 /dag) 5160<br />

Hydrauliske ledningsevner<br />

(m/dag)<br />

Startværdi Maksimum Minimum<br />

Tørv 0,335 86,4 8,64 * 10 -6<br />

Gytje øverst 5,16 * 10 -3<br />

8,64 * 10 -3<br />

8,64 * 10 -6<br />

Gytje nederst 8,90 * 10 -4<br />

8,64 * 10 -3<br />

8,64 * 10 -6<br />

Sand i ådal 0,0163 8,64 8,64 * 10 -3<br />

Morænesand 4,32 8,64 8,64 * 10 -3<br />

Ler 4,32 * 10 -6<br />

8,64 * 10 -6<br />

8,64 * 10 -9<br />

Kalk 10 86 0,86<br />

67


4<br />

Opsætning af<br />

grundvandsmodellen<br />

Som beskrevet i indledningen er formålet med denne rapport at kvantificere den<br />

påvirkning, der er af Hølbækkens vandføring, hvis der oprettes en grundvandsindvinding<br />

i Volsted Plantage. Dette undersøges ved at opstille en numerisk<br />

grundvandsmodel over projektområdet, der bruges til at simulere indvindingen af<br />

grundvand fra Volsted Plantage. I dette kapitel beskrives først, hvordan<br />

grundvandsmodellen opbygges på baggrund af den opsatte konceptuelle model, kapitel<br />

3 (Konceptuel model af projektområdet). Herefter beskrives kalibreringen af modellen<br />

og de ændringer, det har været nødvendigt at foretage for at tilpasse modellen.<br />

Resultater fra grundvandsmodellen præsenteres og det bliver vurderet, hvilken<br />

indflydelse indvindingen af vand i Volsted Plantage vil have på Hølbækkens vandføring<br />

og opland.<br />

4.1 Grundvandsmodellering generelt<br />

En numerisk grundvandsmodel opsættes ved at inddele det ønskede område i en række<br />

celler og opstille en ligning, der beskriver strømningen mellem cellerne. Ligningen<br />

opstilles ved hjælp af vandbalanceligningen og Darcys lov, ligning 4.1.


Kapitel 4 Opsætning af grundvandsmodellen<br />

Ligning 4.1. Styrende ligning (Jensen, 2005).<br />

h<br />

S s K x<br />

t<br />

2<br />

h<br />

K 2 y<br />

x<br />

2<br />

2<br />

h h<br />

K z Q R<br />

2<br />

2<br />

y<br />

z<br />

Hvor:<br />

Ss : Specifikt magasintal (m -1 )<br />

Kx, Ky og Kz : Hydrauliske ledningsevner for de tre akser, X, Y og Z (m/s)<br />

h : Højden af potentialet<br />

Q : Oppumpning m 3 /s<br />

R : Opfyldning m 3 /s<br />

En af parametrene i den styrende ligning er det specifikke magasintal, der angiver, hvor<br />

meget vand en 1 m 3 celle frigiver, hvis vandspejlet sænkes 1 m. Denne frigivelse af<br />

vand skyldes, at sænkningen af vandspejlet mindsker opdriften af de mineralkorn, der<br />

overlejrer grundvandsmagasinet, hvorved der sker en sammenpresning af porerne, så<br />

magasinet ikke kan indeholde så meget vand. Derudover udvider vandet sig ved et<br />

lavere tryk, men dette er af mindre betydning. Det specifikke magasintal er væsentligt i<br />

en dynamisk model, hvor grundvandsspejlet kan variere, mens det normalt sættes til 0 i<br />

stationære modeller (Jensen, 2005).<br />

Beskrivelsen af strømningen fra en celle til den næste forudsætter, at hver celle tildeles<br />

en hydraulisk ledningsevne og et potentiale. Da der her er tale om en tredimensionel<br />

model vil der til hver celle være 6 tilgrænsende celler, hvorimellem strømningen kan<br />

foregå. Da de tilgrænsende celler ikke nødvendigvis har de samme hydrogeologiske<br />

egenskaber, beregnes en samlet hydraulisk ledningsevne for strømningen mellem to<br />

celler. Den samlede hydrauliske ledningsevne beregnes som et harmonisk gennemsnit<br />

over både den hydrauliske ledningsevne for de to celler og længden fra cellekanten til<br />

cellens midtpunkt, ligning 4.2.<br />

Ligning 4.2. Harmonisk gennemsnit (Spitz & Moreno, 1996: 134).<br />

70<br />

Hvor:<br />

Kz<br />

z1 og z2<br />

K1, og K2<br />

K z<br />

z1<br />

z2<br />

<br />

z1<br />

z2<br />

<br />

K K<br />

: Gennemsnitlig hydraulisk ledningsevne (m/s)<br />

: Længden mellem cellens beregningspunkt og cellekanten (m)<br />

: Hydrauliske ledningsevner de to celler (m/s)<br />

1<br />

2


4.2 Grundvandsmodellen<br />

I denne rapport vil grundvandsmodelleringen blive foretaget ved hjælp af programmet<br />

GMS 5.1 (Groundwater Modeling System). Ved opsætningen af grundvandsmodellen er<br />

der truffet en række valg i forhold til modeltype, inddeling af projektområdet og<br />

repræsentation af lag og linser i modellen. Disse vil i det følgende fremlægges.<br />

Modeltype<br />

Grundvandsmodellering kan overordnet foretages på 2 måder, enten med en numerisk<br />

model eller en analytisk model. Det er valgt at anvende en numerisk model i denne<br />

rapport, da den giver mulighed for at bestemme detaljerede variationer i undersøgelsesområdet<br />

ved inddeling i mange celler. En numerisk model bliver opstillet på baggrund<br />

af flere parametre end en analytisk model og giver et mere detaljeret resultat. Den<br />

numeriske model kan desuden bedre repræsentere mindre forekomster i undergrunden,<br />

som det skønnes har stor betydning, eksempelvis gytjelaget i denne rapports model.<br />

Tidsperspektiv<br />

Grundvandsmodeller kan være enten dynamiske, hvor tidsafhængige variable kan<br />

inddrages i udregningerne, eksempelvis varierende oppumpning, vandføring og grundvandsdannelse<br />

eller de kan være stationære, hvor udregningerne sker på baggrund af<br />

gennemsnitlige værdier (Refsgaard & Henriksen, 2005). Den stationære grundvandsmodel<br />

giver således en gennemsnitssituation for modelområdet (Brun, 2005). Det er<br />

valgt at opsætte en stationær grundvandsmodel, da formålet med modellen er at vise,<br />

hvordan vandføringen i Hølbækken gennemsnitligt bliver ændret ved en ændring af<br />

oppumpning af vand i Volsted Plantage. Derudover er der ikke er tilstrækkeligt<br />

datamateriale over området til, at der kan opsættes en dynamisk model, eksempelvis<br />

mangler der tidsserier over vandføringen i Hølbækken, grundvandsdannelsen, og<br />

variation af grundvandsspejlets højde over tid.<br />

Antal celler<br />

For at løse den styrende ligning er det nødvendigt, at sætte den op i en matrice og<br />

løsningen kræver derfor, at der er mindst seks celler. For en grundvandsmodel med<br />

eksempelvis 1000 celler bliver det til en matrice med 1000 ligninger, hvor der er 1000<br />

ubekendte (Spitz & Moreno, 1996: 138). I GMS bliver dette matricesystem løst ved<br />

hjælp af en løsningsalgoritme, der beregner en stationær tilstand for matricesystemet.<br />

Beregningen sker iterativt, hvilket vil sige, at beregningerne gentages til der ikke er<br />

større forskel mellem løsningerne end et angivet interval. Løsningen af matricesystemet<br />

er forbundet med en vis fejl, der opstår på baggrund af den måde, algoritmen løser det<br />

opsatte matricesystem på. Den stationære tilstand, der beregnes, angiver koten på<br />

71


Kapitel 4 Opsætning af grundvandsmodellen<br />

grundvandsspejlet og flowet gennem cellerne. I en stationær grundvandsmodel vil det<br />

specifikke magasintal derfor ikke være væsentligt.<br />

Horisontal diskretisering<br />

Ved diskretiseringen af et modelområde er det vigtigt at gøre sig overvejelser over,<br />

hvordan projektområdet inddeles, så den opsatte grundvandsmodel bliver præcis og<br />

konvergent. Det er nødvendigt at foretage en forsimpling af det naturlige system, da det<br />

ellers ikke er muligt, at formulere transportproblemet numerisk (Sonnenborg, 2005).<br />

Samtidig skal der indgås et kompromis mellem modellens detaljeringsgrad og<br />

beregningstiden ud fra de tilgængelige data (Spitz & Moreno 1996: 216). Der findes<br />

overordnet to måder projektområdet kan diskretiseres på, enten ved finite differens<br />

metoden, hvor projektområdet inddeles i celler med kvadratisk grundflade og varierende<br />

højde på de enkelte cellers sider eller ved finite element metoden, hvor projektområdet<br />

inddeles i trekantede celler eller polygoner (Brun, 2005). Med finite differens metoden<br />

bliver potentialet beregnet i cellens midtpunkt, mens det med finite element metoden<br />

bliver beregnet i hjørnerne. I områder hvor det er vurderet, at undergrundes geologi er<br />

kompleks, kan finite element metoden have en fordel, men sædvanligvis er der ingen<br />

forskel i det resultat, der beregnes med de to metoder (Brun, 2005). Til denne<br />

grundvandsmodellering er det valgt at anvende finite differens metoden, da der gennem<br />

semestret er opnået erfaring i denne metode.<br />

72<br />

Valg af cellestørrelse - Mængden af vand, der bevæger sig gennem en celle,<br />

afhænger af den størrelse, cellen er sat til. Det betyder, at modellens beregninger<br />

bliver mere upræcise og opløsningen på simuleringerne bliver mere grovkornet,<br />

jo større cellestørrelse, der vælges at arbejde med. I forhold til valg af cellestørrelse<br />

er det vigtigt, at opløsningen er finest i de områder, der har særlig<br />

interesse i modelleringen, såsom områder med vandløb, indvindingsboringer<br />

eller en stor grundvandsgradient (Spitz & Moreno, 1996: 220).<br />

Som tidligere nævnt, består undergrunden i projektområdet af opsprækket kalk,<br />

hvilket giver et specielt problem i forhold til cellestørrelsen. Problemet opstår,<br />

da strømningen i kalk vil foregå hurtigt i sprækkerne og langsomt i matrixet.<br />

Hvis et lille udsnit af kalken betragtes, kan der enten foregå hurtig strømning<br />

eller slet ingen strømning, afhængig af om dette udsnit indeholder en sprække.<br />

Derfor er det nødvendigt i modellen at vælge en cellestørrelse, der er stor nok til,<br />

at sprækkerne er nogenlunde ensartet fordelt mellem cellerne, så strømningen<br />

ikke afgøres af enkelte sprækker. I denne model er det nogle steder valgt, at<br />

kalken skal inddeles i 100 m x 100 m celler, hvorfor sprækkestørrelsen i forhold<br />

til cellestørrelsen ikke antages at have den store indflydelse på den måde,<br />

strømningen gengives på, men teoretisk set er der mulighed for, at denne


cellestørrelse ikke kan gengive strømningen helt. Denne antagelse underbygges<br />

endvidere af den flowlog, der er blevet optaget i projektområdet og som viser, at<br />

udstrømningen i 80 meter stort set er jævn, bilag A.<br />

Hølbækkens ådal er fokusområdet for denne rapport og da der ud fra egne forsøg<br />

findes detaljerede data herfra, er det besluttet at diskretisere dette område med<br />

en finere cellestørrelse end i den resterende del af projektområdet, figur 4.1.<br />

Den grovere inddeling i resten af området er valgt pga. detaljeringsgraden af de<br />

geologiske data, hvor en finere inddeling ikke ville give et mere nøjagtigt<br />

resultat, men kun resultere i en længere beregningstid. Den overordnede<br />

cellestørrelse i projektområdet er sat til 100 m, mens den i Hølbækkens ådal er<br />

sat til 5 m. Der må ikke være for stor forskel mellem to tilstødende celler, da<br />

beregningen herved bliver upræcis (Spitz & Moreno, 1996: 220). Derfor stiger<br />

cellestørrelsen stiger gradvist med 10 % for hver celle fra Hølbækkens ådal,<br />

indtil de når 100 m, som resten af projektområdet. Der kan opstå problemer med<br />

disse små celler i forhold til sprækker i kalken, men det antages, at placeringen<br />

af de store sprækker i ådalens kalk fremgår af kilderne og derfor er simuleret i<br />

modellen.<br />

Orientering af cellerne - I beregningerne af grundvandsstrømningen er det<br />

væsentligt, hvordan cellerne er orienteret i forhold til den grundvandsstrømning,<br />

de skal simulere. Det skyldes, at den fremherskende strømning splittes op i to<br />

vektorer, hvis ikke cellerne er orienteret parallelt med og vinkelret på den<br />

fremherskende strømningsretning. Denne opsplitning resulterer i numerisk<br />

spredning af grundvandsflowet, hvilket ikke repræsenterer situationen i<br />

virkeligheden. Beregningscellerne er i den opsatte grundvandsmodel blevet<br />

orienteret parallelt med og vinkelret på Lindenborg Å og Hølbækken, figur 4.1.<br />

73


Kapitel 4 Opsætning af grundvandsmodellen<br />

74<br />

Figur 4.1. Diskretiseringen af projektområdet i den finite differens model, der er opsat i<br />

rapporten.<br />

Diskretisering af de geologiske lag<br />

Den vertikale diskretisering kan overordnet foregå på to måder (Henriksen et al., 2001:<br />

26-27):<br />

Vektormetoden – cellerne følger de geologiske lag. Således opdeles cellerne<br />

efter lag med samme eller lignende hydrogeologiske egenskaber. Denne tilgang<br />

kan med fordel bruges på grundvandsmagasiner som er tolket til at have en<br />

forholdsvis enkelt opbygget undergrund.<br />

Pixelmetoden – cellerne er (forholdsvis) rektangulære. Denne metode er bedre<br />

egnet til områder med en kompleks sammensat undergrund. For hver celle<br />

findes det dominerende materiale. Ulempen er, at små lag, der kan være af stor<br />

betydning for grundvandsstrømningen, eksempelvis på grund af lille hydraulisk<br />

ledningsevne, ikke bliver inddraget i modellen.<br />

Den geologiske model af projektområdet viste, at undergrunden er opdelt i to områder:<br />

morænen og ådalen. Begge dele består i høj grad af sammenhængende lag, hvorfor<br />

vektormetoden vil blive anvendt i grundvandsmodelleringen. Den vertikale<br />

diskretisering af projektområdet følger princippet i figur 4.2. Figuren viser at kalken er


opdelt i fire lag. Denne inddeling blev gjort, idet der var en forventning om, at kalkens<br />

hydrauliske ledningsevne blev mindre med dybden pga. kompaktion. Resultaterne fra<br />

flowlogundersøgelsen nord for Hølbækken viste, at der er vertikale zoner med forskellig<br />

hydraulisk ledningsevne i kalken, men at ledningsevnen ikke nødvendigvis er<br />

aftagende. Ud fra dette viste det sig fornuftigt at inddele kalken i zoner. Den nedre<br />

grænse af kalklaget er sat til kote -140, hvor TEM målinger over området viste, at der<br />

findes saltvand, og derfor antages det, at vandudvekslingen i denne dybde er ubetydelig.<br />

Dette underbygges af, at der ikke sker en blanding af grundvandet og saltvandet, fordi<br />

saltvand har en højere densitet end ferskvand og da den hydrauliske ledningsevne for<br />

kalklaget er lav i denne dybde.<br />

Tørv<br />

Gytje 1<br />

Gytje 2<br />

Ådal<br />

Sand<br />

Moræne<br />

Figur 4.2. Principskitse af den vertikale diskretisering af projektområdet.<br />

Lerlinse<br />

Sand 1<br />

Sand 2<br />

Sand 3<br />

Sand 4<br />

Kalk 1<br />

Kalk 2<br />

Kalk 3<br />

Kalk 4<br />

Saltvand<br />

De lag der findes over kalken er opdelt i ådal og moræne. I ådalen er der under tørven<br />

indsat to gytjelag med ens højde og forskellig hydraulisk ledningsevne, da gennemstrømningsforsøgene<br />

viste, at gytjen havde en variation i hydraulisk ledningsevne.<br />

Herunder findes sandet. Morænen er opdelt i fire lag med samme hydrauliske<br />

ledningsevne, for at kunne indpasse de lerlinser og siltforekomster, der er fundet i<br />

området. Lerlinserne er placeret i det lag, hvor de passer bedst og tykkelsen på laget er<br />

korrigeret tilsvarende. Inddelingen er valgt for at kunne repræsentere lerlinserne, med<br />

75


Kapitel 4 Opsætning af grundvandsmodellen<br />

lavere hydraulisk ledningsevne, hvilket kan ændre strømningsvejene. Tykkelsen på<br />

lagene går fra 1 til 10 m, dvs. at ler- og siltforekomster på under ca. 1 m tykkelse ikke er<br />

taget med og ler og silt forekomsterne er slået sammen. Den horisontale udbredelse af<br />

linserne er interpoleret ud fra boringerne og understøttes af tolkningen af terrænet og<br />

dannelsesmiljøet.<br />

Vandløb, dræn og kilder<br />

Både vandløbene og drænene i projektområdet er modelleret ved at anvende<br />

drænfunktionen i GMS, mens vandløbsfunktionen ikke anvendes. Dette skyldes, at<br />

vandløbsfunktionen tillader strømning både ind og ud af vandløbet, mens<br />

drænfunktionen kun fjerner vand, hvilket betyder, at drænfunktionen bedst repræsenterer<br />

situationen i projektområdet. I modellen er vandløbene og drænene blevet<br />

indlagt, så de dækker hele ådalen. Dette er gjort for at afspejle situationen i ådalen, der<br />

gennem en lang årrække er blevet drænet, så der sandsynligvis er god kontakt mellem<br />

de øverste lag i ådalen og vandløbet. De kilder, der er observeret i ådalen, er modelleret<br />

ved at indsætte en søjle af celler med en høj hydraulisk ledningsevne ind i fokusområdet.<br />

Herved fremkommer der en fri strømningsvej mellem det dybe regionale<br />

grundvandsmagasin og ådalen, hvilket sandsynligvis svarer til situationen i virkeligheden,<br />

hvor kilderne formodentlig skyldes sprækker i kalken. De kilder, der er taget<br />

med i modellen, er kilde B, C, D, G og søen, jf. figur 3.4, hvilket er skønnet til at være<br />

de mest betydende.<br />

4.3 Kalibrering af grundvandsmodellen<br />

En vigtig del af det at arbejde med grundvandsmodellering er at tilpasse modellen de<br />

virkelige forhold, hvilket sker gennem en kalibrering af grundvandsmodellens<br />

parametre. Kalibreringen foregår ved, at variere modellens parametre bestemme den<br />

forskel i potentialet, der beregnes i forhold til det observerede. Kalibreringen skal<br />

således komme frem til den kombination af parameterværdier, hvor modellens afvigelse<br />

i forhold til virkeligheden er mindst mulig (Henriksen, 2005). Dette er nødvendigt, da<br />

de parameterværdier, der til at begynde med er indsat, bygger på forsøg med fejlkilder<br />

og på gennemsnitsbetragtninger af området.<br />

De første gennemkørsler af modellen gav store afvigelser mellem observeret og<br />

beregnet potentiale. Dette kan skyldes, at de anvendte data for både hydraulisk<br />

ledningsevne og potentiale, ikke er repræsentative, eller at opsætningen af den<br />

geologiske tolkningsmodel afviger fra virkeligheden. I grundvandsmodellen er<br />

værdierne for hydrauliske ledningsevner for de forskellige lag derfor kalibreret, og det<br />

modellerede grundvandsspejl er sammenlignet med pejlinger af grundvandsspejlet<br />

76


fundet i GEUS’ Jupiter database. For at sikre nøjagtigheden af pejledataen til<br />

kalibreringen er der kun anvendt registreringer fra efter 1980, samt to målinger fra 1972<br />

og 1973, i et område, hvor der ikke fandtes nyere data. Pejlingerne kan variere over 2 m<br />

over et år (Søndergaard, 2003: B4). Derfor er målet med kalibreringen at<br />

grundvandsspejlet skal passe med pejlingerne inden for ca. 2 m. Kalibreringen af<br />

grundvandsmodellen er gjort efter trial-and-error-metoden, hvilket vil sige, at<br />

parameterværdierne for den hydrauliske ledningsevne blev varieret til forskellen mellem<br />

observeret og beregnet potentiale var mindst mulig. Som udgangspunkt blev de værdier,<br />

der er fundet igennem kapitel 3 (Konceptuel model af projektområdet), anvendt. En<br />

opsummering af disse kan ses i tabel 3.9. Samtidig er den geologiske model blevet<br />

holdt tilnærmelsesvis konstant bortset fra, at der er ændret lidt på placeringen af<br />

drænene og indlagt nogle kilder. Potentialedataene er anvendt som punktdata, da der er<br />

for store usikkerheder omkring en interpoleret overflade lavet ud fra punkter<br />

(Henriksen, 2005). Dette skyldes, at der i interpoleringen ikke bliver taget højde for de<br />

geologiske strukturer og de højdeforhold, der findes mellem målingerne af potentialet.<br />

Der er 11 parametre i modellen, der er blevet kalibreret. Parametrene er hydrauliske<br />

ledningsevner for de forskellige aflejringer i projektområdet, konduktansen i drænene<br />

og ledningsevne for sprækker til kilder i Hølbækken. For alle lag er det valgt at den<br />

horisontale og vertikal ledningsevne er ens, dvs. anisostropien i alle retninger er sat til 1.<br />

Det ville have taget meget lang tid, hvis alle parametrene skulle kalibreres samtidig.<br />

Derfor er parametrene blevet kalibreret en eller flere af gangen, mens de øvrige<br />

parametre er blevet fastholdt. Kalklagets hydrauliske ledningsevne var mest usikker og<br />

derfor blev dette kalibreret først, da præcisionen forbedres ved at kalibrere mod de mest<br />

sikre parametre. De andre parametre er blevet kalibreret efterfølgende i en rækkefølge<br />

afhængig af, hvor stor deres usikkerhed er vurderet til at være. I sidste ende var det stort<br />

set udelukkende kalken, der har indflydelse på grundvandsstrømningen i projektområdet,<br />

hvorfor kalibreringen fokuserede mest på at få dette lag kalibreret rigtigt.<br />

Herunder følger en beskrivelse af, hvordan kalibreringen af de forskellige aflejringer er<br />

udført. Resultatet af kalibreringen fremgår af tabel 4.1.<br />

Kalklaget – Kalklaget ligger under grundvandsspejlet i det meste af projektområdet<br />

og derfor er kalklaget det vigtigste lag i forhold til vandtransporten, da<br />

det er i dette lag, den største strømning foregår. Som udgangspunkt, blev det<br />

antaget, at sprækkerne i kalken befandt sig i den øvre del af kalklaget. Derfor<br />

blev kalken øverst opdelt i 3 lag på hver 10 m, mens resten af kalklaget ned til<br />

kote -140 m blev sat til lag fire. Ved nærmere undersøgelse af flowloggen viste<br />

det sig, at der var en høj grundvandstrømning dybere ned i kalken, hvilket<br />

indikerede at sprækkerne strakte sig dybere ned i undergrunden. Derfor blev der<br />

lavet en ny inddeling af kalklagene ud fra flowloggen, hvor det øverste lag blev<br />

77


Kapitel 4 Opsætning af grundvandsmodellen<br />

78<br />

sat til 10 m, det næste 30 m, det tredje 40 m, mens det sidste lag igen blev sat til<br />

kote -140 m og varierer således mellem 0 og 125 m tykkelse. Disse to bud på<br />

kalkens opbygning gav væsentligt forskellige hydrauliske ledningsevner under<br />

kalibreringen, men deres resultat i forhold til Hølbækkens vandføring var<br />

nogenlunde ens. Det er derfor valgt at bibeholde begge opdelinger i det videre<br />

forløb, for kunne komme med forskellige tolkninger af vandindvindingens<br />

påvirkning på Hølbækken. De vil herefter blive refereret til som henholdsvis<br />

model 1, hvor tykkelsen på de tre øverste lag er ens og model 2, hvor tykkelsen<br />

på de øverste lag er afgrænset ud fra flowloggen, figur 4.3.<br />

10 m K1<br />

10 m K2 < K1<br />

10 m K3 < K2<br />

K4 < K3<br />

Kote - 140<br />

10 m K * 1<br />

30 m K * 2 > K * 1<br />

40 m K * 3 > K * 2<br />

Model 1 Model 2<br />

K * 4


Morænen – Kalibreringen af sandet i morænen viste, at denne kun havde<br />

betydning for den midterst østlige del af projektområdet, nordøst for<br />

Hølbækken. Leret udgør kun mindre lokale linser og derfor har den hydrauliske<br />

ledningsevne i leret ikke afgørende betydning for den samlede vandføring, der<br />

beregnes i modellen. Af dette tolkes det, at den hydrauliske ledningsevne for ler<br />

kan variere uden at det har stor betydning for modelresultatet.<br />

Ådalene – De første gennemkørsler af modellen, med de hydrauliske<br />

ledningsevner, der blev fundet i laboratoriet på de intakte prøver, viste, at de<br />

målte værdier var for lave, idet der kom for lidt vand ud til åerne. Ved at<br />

sammenholde de laboratoriebestemte ledningsevner med ledningsevnerne fundet<br />

ved slugtestene, blev det fundet, at der var op til to dekaders forskel.<br />

Ledningsevnen for jordtyperne i ådalen blev derfor sat til ti gange højre end de<br />

laboratoriebestemte. Jensen (2002) skriver, at der kan forekomme fejl mellem<br />

parametre fundet i et område i naturen i forhold til laboratoriebestemte. Fejlen<br />

ligger i, at parametrene er målt på forskellig skala (Jensen 2002: 33-34). Dette<br />

kan forklare disse forskelle. Under kalibreringen blev de to gytjelag, og sandet i<br />

ådalen slået sammen for at simplificere kalibreringen. Ledningsevnen herfor kan<br />

betragtes som den gennemsnitlige ledningsevne for disse tre lag.<br />

Dræn - Konduktansen for drænene i projektområdet kendes ikke nøjagtig. Det<br />

er blevet skønnet, at denne vil være af samme størrelse som tørven, da<br />

bestemmelsen af åbundens hydrauliske ledningsevne viste, at åbunden havde<br />

den samme hydrauliske ledningsevne som tørven. Derfor er der valgt at sætte<br />

drænenes konduktans højt, således at tørven vil være begrænsende for, hvor<br />

meget vand der vil strømme ud i disse.<br />

Kilder - Kilderne skal fungere som strømningsveje for grundvandet i kalken op<br />

til jordoverfladen. Derfor er dennes hydrauliske ledningsevne sat højt, for ikke at<br />

være begrænsende.<br />

I figur 4.4 ses en sammenligning af det observerede og beregnede grundvandsspejl for<br />

begge modeller. Det ses, at der for model 2 er større overensstemmelse med<br />

observationspunkterne, end for model 1. Forskellen mellem de to modeller ligger i den<br />

forskellige tolkning af kalken.<br />

79


Kapitel 4 Opsætning af grundvandsmodellen<br />

Tabel 4.1. De hydrauliske ledningsevner, der er bestemt i kalibreringen for de to modeller.<br />

80<br />

Materiale Model 1 Model 2<br />

(m/d) (m/d)<br />

Sand 4 4<br />

Ler 0,00432 0,00432<br />

Tørv 3,55 3,55<br />

Gytje og ådals sand 0,0516 0,0516<br />

Kalklag 1 12 20<br />

Kalklag 2 9 20<br />

Kalklag 3 7 30<br />

Kalklag 4 5 0,5<br />

Kalklaget i det nordøstlige hjørne 5 3<br />

Kilder 100 100<br />

Andre hydrauliske parametre<br />

Konduktans for dræn 10 10<br />

Figur 4.4. Kalibrering af modellerne. For hvert observationspunkt ses usikkerhedsintervallet på 2 m,<br />

størrelsen og farven på søjlen angiver hvor godt det modellerede grundvandsspejl passer med<br />

observationen (grønt er inden for sikkerhedsintervallet, rødt er uden for). De to blå farver på gridet<br />

angiver den horisontale opdeling af de tre øverste lag kalk.


4.4 Placering af indvindingsboringerne<br />

Der foreligger ingen nærmere planer for placering af indvindingsboringerne, blot at<br />

Volsted Plantage er udset til dette formål. I dette afsnit vil der ses nærmere på nogle af<br />

de hensyn, der skal tages ved placeringen af indvindingsboringerne. Først undersøges<br />

kort de juridiske og økonomiske forhold, der har indvirkning på placeringen. Derefter<br />

undersøges placeringen i forhold til potentialeforhold. Dette munder ud i en mulig<br />

placering af indvindingsboringerne.<br />

4.4.1 Juridiske og økonomiske hensyn<br />

Volsted Plantage har ikke været brugt til landbrug mindst de sidste 100 år jf. afsnit<br />

3.1.1 (Kortanalyse) og derfor er der formentlig ikke blevet udledt gylle (nitrat) eller<br />

pesticider, med undtagelse af et mindre område i den sydvestlige del af plantagen, hvor<br />

der dyrkes juletræer. Ydermere er der kun få huse, der udleder spillevand i nærheden.<br />

Plantagen giver således mulighed for en placering væk fra evt. forureningskilder.<br />

Miljøbeskyttelsesloven foreskriver, at der kan fastlægges et beskyttelsesområde, kaldet<br />

kildepladszone, indenfor hvilken det ikke er tilladt at aflede spildevand eller udføre<br />

anden aktivitet, der kan give anledning til at forurene grundvandet. Denne zone er<br />

sædvanligvis 300 m (Miljøbeskyttelsesloven, 2006: § 22; Miljøstyrelsen, 2000: 18). Der<br />

er ingen direkte forbud mod landbrugsdrift med brug af pesticider eller nitrat inden for<br />

kildepladszonen.<br />

Indenfor indsatsområder, der omfatter det hydrologiske opland og dermed også<br />

kildepladszonen, kan der indføres regulering af pesticider og nitrat. Indsatsområder<br />

udpeges, hvor der er behov for at yde en indsats for at sikre den fremtidig drikkevandsforsyning.<br />

Der findes ingen metoder til at udpege pesticidfølsomme områder. Nitratfølsomme<br />

områder udpeges, hvis f.eks. grundvandsmagasinet består af kalk og der ikke<br />

findes noget beskyttende lag (Miljøstyrelsen 2000: s. 13-17). I projektområdet består<br />

grundvandsmagasinet af kalk og der er ingen beskyttende lag, da der kun forekommer<br />

enkelte lerlinser, jf. kapitel 3.2 (Projektområdets geologi). Derfor vil det hydrologiske<br />

opland til indvindingsboringen blive klassificeret som Område med Særlige Drikkevandsinteresser<br />

(OSD), og der skal iværksættes afbødende foranstaltninger til beskyttelse<br />

af drikkevandsressourcen.<br />

Indenfor projektområdet, er kun husene i Volsted koblet til den kommunale<br />

spildevandsopsamling (AK, 2005). Derfor forventes det, at husene nær Voldsted<br />

Plantage har egne nedsivningsanlæg. Dette betyder, at indvindingsboringerne minimum<br />

skal placeres 300 m væk disse nedsivningsanlæg. Det er muligt at ændre nedsivningsanlæggenes<br />

beliggenhed eller koble husstande til spildevandsopsamlingen, hvilket dog<br />

vil medføre omkostninger.<br />

81


Kapitel 4 Opsætning af grundvandsmodellen<br />

I Volsted Plantage findes der mere end ti gravhøje (NJA, 2006), hvor der er en 100 m<br />

beskyttelseszone omkring, inden for hvilken der ikke må foretages ændringer i<br />

tilstanden (Naturbeskyttelsesloven, 2004: § 18). Der kan søges dispensation fra 100 m<br />

zonen, og det vurderes, at indvindingsboringerne ikke vil skæmme gravhøjene, så længe<br />

de ikke placeres oveni dem. Af denne grund vil der ses væk fra gravhøjene under den<br />

forudsætning, at selve gravhøjene bliver inddraget i detailplanlægningen af placeringen.<br />

Tages økonomien i betragtning vil det være fornuftigt at placere boringerne tæt på en<br />

vej, så det er let at komme til med boremaskiner mv. Det er selvfølgelig muligt at lave<br />

en ny vej, men det vil være billigere at benytte eksisterende veje.<br />

4.4.2 Placering af boringer i forhold til grundvandets strømningsretning<br />

Når der indvindes grundvand, vil der omkring boringen ske en sænkning af<br />

grundvandsspejlet. Placeres flere boringer langs strømningsretningen vil disse forstærke<br />

hinandens sænkning, hvilket ikke er tilfældet, hvis de placeres på tværs af<br />

strømningsretningen, se figur 4.5.<br />

Figur 4.5. Placering af boringer i forholdt til strømningsretningen. Tv. Boringerne er placeret på langs af<br />

strømningsretningen, sænkningen forstærkes omkring den øverste boring. Oplandet til den nederste<br />

boring ligger uden for oplandet til den øverste boring. Th. Boringerne er placeret på tværs af<br />

strømningsretningen, og oplandene ligger parallelle uden at de påvirker hinanden.<br />

Grundvandssænkningen begrænses, når der etableres flere boringer, da sænkningen<br />

herved fordeles over et større område, i stedet for at være koncentreret et sted. I<br />

modellen er det valgt at bruge 5 boringer, der hver skal indvinde maksimalt 200.000<br />

m 3 /år. Dette er gjort efter model fra Aalborg Kommunes nye kildeplads Kongshøj, 10<br />

km nord for Volsted, der tilnærmelsesvis indvinder denne mængde. I Kongshøj består<br />

82<br />

Strømningsretning


grundvandsmagasinet ligeledes af kalk og da det ligger forholdsvis tæt på Volsted,<br />

skønnes det rimeligt, at dette maksimum for pumpekapacitet kan overføres direkte.<br />

4.4.3 Valg af Placering<br />

På baggrund af de beskrevne forbehold er indvindingsboringernes placering bestemt.<br />

Der er taget udgangspunkt i, at der skal være en kildepladszone på 300 m. Derfor skal<br />

der være 300 m til nærmeste ejendom. Det er desuden ønsket, at kildepladszonen skal<br />

ligge indenfor Volsted Plantage, hvilket er bestemt som mere end 300 m fra skovkanten.<br />

Endelig skal boringerne være placeret tæt på en eksisterende vej. På figur 4.6 er der<br />

indtegnet bufferzoner på 300 m fra bebyggelse og 300 m fra skovkanten. De ikke<br />

skraverede områder i plantagen udgør således området, hvor det er muligt, at oprette<br />

indvindingsboringerne. Det er valgt at placere indvindingsboringerne langs den<br />

eksisterende vej med 100 m mellemrum. På figur 4.7 ses indvindingsboringernes<br />

placering i forhold til grundvandspotentialet og dermed strømningsretningen.<br />

Figur 4.6. Drikkevandsboringernes placering ud fra krav om mere end 300 meter til beboelse og<br />

skovkanten, samt tæt på en eksisterende vej. Data fra TOP10DK (2000).<br />

De to potentialekort går på tværs af hinanden i plantagen, men det ses, at de oprettede<br />

boringer står i en vinkel på ca. 45º på potentiale kurverne (NJA) og på<br />

potentialekurverne fra egen interpolering (eget) ser det ud som om denne placering er<br />

beliggende på toppen af et vandskel. I forhold til NJAs potentiale er det svært at sige,<br />

83


Kapitel 4 Opsætning af grundvandsmodellen<br />

om sænkningerne vil forstærke hinanden, mens det ikke ser ud til at være tilfældet ved<br />

vores eget interpolerede potentiale. Vores eget interpolerede potentiale angiver det mere<br />

overfladenære grundvandsspejl, idet der også er anvendt koter fra åer og dræn.<br />

Figur 4.7. Drikkevandsboringernes placering i forhold til grundvandspotentialet. Ækvidistance på 1 m.<br />

NJA potentialet blev præsenteret på figur 1.4, eget potentiale blev præsenteret i figur 3.20.<br />

4.5 Anvendelse af grundvandsmodellen<br />

Da formålet med grundvandsmodellerne er at vurdere påvirkningen på Hølbækkens<br />

vandføring, er det væsentligt at vurdere de modellerede vandføringers resultat i forhold<br />

til den formodede gennemsnitsvandføring. Vandføringsmålingerne i afsnit 3.3.1<br />

(Overfladisk afstrømning) viste, hvordan tilstrømningen af vand til Hølbækken er<br />

fordelt. For at kunne sammenligne dette med de beregnede vandføringer, er Hølbækken<br />

delt op i fire områder. Områderne er Hølbækkens ådal, opstrøms vandføringsmåling 2,<br />

5, 6 og 9, figur 3.21, hvor måling 9 er benyttet som mål for den samlede mængde vand,<br />

Hølbækken afleder. I hvert område findes vandføringsdata målt i projektområdet og<br />

vandføringen er ligeledes modelleret for hvert af områderne. Sammenlagt giver<br />

vandføringen i de fire områder den totale vandføring. I tabel 4.2 ses forskellen mellem<br />

den observerede vandføring og de beregnede vandføringer.<br />

84


Tabel 4.2. Forskellen mellem den observerede vandføring og modellernes vandføring.<br />

Område 1 Område 2 Område 3 Område 4 Total<br />

Observeret (m 3 /d) 2.908 876 1.192 186 5.162<br />

Model 1 (m 3 /d) 784 572 985 1.436 3.777<br />

Forskel (%) -73 -35 -17 672 -27<br />

Model 2 (m 3 /d) 427 430 900 1.289 3.048<br />

Forskel (%) -85 -51 -24 593 -41<br />

Begge modeller underestimerer den totale vandføring i Hølbækken med ca. 1/5 i forhold<br />

til den observerede vandføring, der antages at være den gennemsnitlige vandføring.<br />

Område 1, ligger længst opstrøms og er kraftigt underestimeret. Område 2 og 3 også er<br />

underestimeret, dog i mindre grad, og område 4 er væsentligt overestimeret. I område 1<br />

og 2 var der flere kilder, og vandføringsmålingerne indikerede, at disse stod for en<br />

væsentlig af del af tilstrømningen til Hølbækken. Det er forsøgt at simulere disse kilder<br />

i modellerne, men den store forskel i de observerede og målte resultater tyder på, at de<br />

indlagte kilder ikke bidrager med tilstrækkelig vand i modellerne. Modellerne har den<br />

største tilstrømning til Hølbækken i område 3 og 4, hvilket tyder på, at kontakten<br />

mellem grundvandet og Hølbækken er modelleret til at være bedre end den er i<br />

virkeligheden. Da områderne i modellerne viser en anden fordeling end den målte, vil<br />

der herefter kun ses på ændringerne i samlet vandføring, og der vil ikke blive set på en<br />

absolut nedgang i vandføringen, men kun en forholdsmæssig.<br />

4.6 Effekt af grundvandsindvindingen på vandføringen i<br />

Hølbækken<br />

Ud fra de to grundvandsmodeller og med udgangspunkt i den placering af indvindingsboringerne,<br />

der blev præsenteret i foregående afsnit, vil det i dette afsnit blive<br />

undersøgt, hvilken indvirkning grundvandsindvindingen har på Hølbækkens vandføring.<br />

Der bliver taget udgangspunkt i en indvinding på 1.000.000 m 3 /år og herefter ændres<br />

størrelsen på indvindingen i et forsøg på at bestemme, hvor stor indvindingen kan være,<br />

for maksimalt at påvirke Hølbækken med 10 %. Indvindinger på under 500.000 m 3 /år<br />

bliver ikke undersøgt, da det vurderes, at det ikke vil være rentabelt at indvinde<br />

mængder under dette. Efterfølgende bliver der set på, hvordan indvindingen af<br />

grundvand påvirker potentialet i oplandet til Hølbækken.<br />

4.6.1 Ændring i vandføring<br />

I tabel 4.3 ses de modellerede påvirkninger af Hølbækkens vandføring i de to modeller<br />

ved forskellige indvindingsmængder. Hølbækken har i model 1 som udgangspunkt en<br />

vandføring, der er ca. 19 % højere end i model 2. Bortset fra denne forskel mellem<br />

85


Kapitel 4 Opsætning af grundvandsmodellen<br />

modellerne, reagerer de meget ens på de forskellige vandindvindinger. Den geologiske<br />

tolkning spiller altså ikke den store rolle i forhold til den procentvise ændring af<br />

vandføringen i de to modeller.<br />

Tabel 4.3. Vandindvindingens reduktion af Hølbækkens vandføring.<br />

86<br />

Indvundet (m 3 /år)<br />

Vandføring<br />

(m 3 /døgn)<br />

Model 1 Model 2<br />

Reduktion af<br />

vandføring<br />

(%)<br />

Vandføring<br />

(m 3 /døgn)<br />

Reduktion af<br />

vandføring<br />

(%)<br />

1.000.000 3.034,0 19,7 2.470,8 18,9<br />

800.000 3.180,1 15,8 2.583,2 15,3<br />

700.000 3.254,4 13,8 2.639,7 13,4<br />

650.000 3.290,6 12,9 2.667,8 12,5<br />

600.000 3.321,8 12,0 2.693,2 11,7<br />

550.000 3.364,8 10,9 2.727,6 10,5<br />

500.000 3.375,3 10,6 2.757,8 9,5<br />

0 3.776,7 0,0 3.048,4 0,0<br />

Den modellerede procentvise reduktion af Hølbækkens vandføring er illustreret i figur<br />

4.8. Her ses det, at begge modeller reagerer lineært på øget indvinding. Indvirkningen<br />

på vandføringen i model 1 ligger ca. 0,5 procentpoint over model 2. I begge modeller<br />

kan der ikke indvindes mere end ca. 500.000 m 3 /år, hvis den gennemsnitlige vandføring<br />

ikke skal overskride 10 %.<br />

Reduktion (%)<br />

20<br />

15<br />

10<br />

5<br />

0<br />

0 200 400 600 800 1.000<br />

Indvundet (1.000 m 3 /år)<br />

Figur 4.8. Reduktionen af Hølbækkens vandføring i forhold til den indvundne vandmængde.<br />

4.6.2 Hølbækkens opland<br />

Model 1<br />

Model 2


I figur 4.9 er det vist, hvordan potentialet og oplandet ser ud før og efter indvinding af<br />

1.000.000 m 3 /år grundvand fra Volsted Plantage. Det ses, at der er en forskel mellem de<br />

to modellers udregnede udgangspotentiale, hvilket betyder, at de to geologiske<br />

tolkninger har haft en betydning i forhold til udregningen af potentialet. I model 1 ligger<br />

potentialet generelt højere i den nordlige del af projektområdet i forhold til model 2,<br />

mens model 2 i dette område har en større gradient. Denne forskel skyldes, at kalklaget<br />

har forskellige hydrauliske ledningsevner og tykkelser i de to modeller. I model 1 har de<br />

tre øverste kalklag hydrauliske ledningsevner mellem 7-12 m/d og området i det<br />

nordøstlige hjørne ligger på 5 m/d. I model 2 er den hydrauliske ledningsevne for de tre<br />

øverste kalklag sat til mellem 20-30 m/d mens kalklaget i det nordøstlige hjørne er 0,5<br />

m/d. De øverste kalklag i denne model er desuden tykkere. I model 2 er der således<br />

større forskel i den hydrauliske ledningsevne mellem kalklagene, hvilket er en af<br />

forklaringerne på, at der ses denne større gradient på potentialet.<br />

Hølbækkens hydrologiske opland er tidligere defineret, jf. figur 3.20, og dette skiller<br />

sig væsentligt fra oplandene fundet i figur 4.9. Grunden til denne forskel er muligvis, at<br />

der er forskelle i den hydrauliske ledningsevne for kalken, der ikke bliver gjort rede for i<br />

modellen. Det kan også skyldes, at der findes fejl i potentialekortet for det tidligere<br />

definerede hydrologiske opland. Potentialekortet er eksempelvis blevet interpoleret,<br />

hvilket giver usikkerheder, idet der blot udregnes et gennemsnitligt potentiale mellem<br />

en række punkter. Yderligere kan grundvandsspejlet variere ca. 2 m over et år, hvorfor<br />

der kan være en variation i pejlingerne. I grundvandsmodellen bliver potentialet regnet<br />

ud på baggrund af antagelser om, hvordan grundvandet i projektområdet strømmer,<br />

hvilket må være mere præcist, idet der bliver taget højde for områdets skiftende geologi.<br />

Der er også en uoverensstemmelse mellem det hydrologiske opland, der beregnes for<br />

modellen og det topografiske opland, der tidligere er præsenteret, jf. figur 3.20. Da der i<br />

det topografiske opland findes ler i undergrunden, er der mulighed for, at der sker en<br />

overfladenær tilstrømning af vand til Hølbækken, som der ikke bliver gjort rede for i<br />

den opsatte grundvandsmodel. Dette anses imidlertid ikke som et problem, idet de<br />

indvindingsboringer, der skal etableres, indvinder det dybe grundvand.<br />

87


Kapitel 4 Opsætning af grundvandsmodellen<br />

Figur 4.9. Projektområdet med de beregnede potentialer, oplandene til Hølbækken og<br />

indvindingsboringerne, indvindingen er sat til 1.000.000 m 3 /år.<br />

4.6.3 Ændringen af Hølbækkens opland<br />

Størrelsen og formen på indvindingsboringernes opland er blevet beregnet ved hjælp af<br />

partikelbaner i de to modeller. Det ses i figur 4.9, at oplandet til indvindingsboringerne<br />

for begge modeller er langstrakte og smalle, men at oplandet, der udregnes for model 2<br />

er smallere og kortere end for model 1. Grunden til den forskellige størrelse på indvindingsoplandene<br />

er sandsynligvis, at koten på indvindingsboringerne i de to modeller er<br />

88


forskellig og at kalkens hydrauliske ledningsevne er forskellig. I begge modeller er<br />

indvindingen lagt ind i det tredje kalklag, men dette lag starter 30 m ned i kalken i<br />

model 1, mens det i model 2 starter 70 m nede i kalken. Det indvundne vand i model 1<br />

er derfor i høj grad det overfladenære grundvand, mens det i model 2 både er det<br />

overfladenære og det dybe. Af denne grund vil en større del af grundvandet i model 1<br />

kunne strømme under boringerne. Potentialeforholdene ændrer sig kun fra boringens<br />

dybde og op, hvorfor oplandet til boringen i model 1 er større end oplandet i model 2,<br />

ved samme mængde indvundet vand. I begge modeller går oplandet til<br />

indvindingsboringerne nord om Hølbækken, hvilket formodentlig skyldes, at der i<br />

modellen er indlagt en række lavpermeable lag i ådalen, som presser oplandet i denne<br />

retning.<br />

Oplandet til Hølbækken er optegnet på baggrund af de modellerede potentialekort og<br />

det beregnede areal ses i tabel 4.4. Det ses, at Hølbækkens opland er større i model 1<br />

end i model 2, hvilket var forventet, da beregningen af Hølbækkens vandføring i model<br />

1 var større end i model 2. Der er imidlertid ikke stor forskel i den procentvise reduktion<br />

af oplandet til Hølbækken før og efter indvindingen, hvilket stemmer godt overens med,<br />

at der ikke er stor procentvis forskel på udregningerne af indvindingens påvirkning af<br />

vandføringen.<br />

Tabel 4.4. Ændring af arealet af Hølbækkens opland.<br />

Model 1 5,32 km 2<br />

Model 2 4,74 km 2<br />

Før Efter Forskel<br />

4,31 km 2<br />

3,94 km 2<br />

18,8 %<br />

16,8 %<br />

Det er interessant at bemærke, at der er en forskel i responstid inden for de to<br />

indvindingsopland fra de to modeller, udregnet ved hjælp af en partikelbanemodel. I<br />

model 1 tager det over 150 år før vandpartiklerne længst væk når indvindingsboringerne,<br />

mens det for model 2 tager under 10 år. Grunden til denne relativt store<br />

forskel i responstid er sandsynligvis, som nævnt tidligere, at de øverste kalklag i model<br />

2 er tolket til at være tykkere og til at have en lavere hydraulisk ledningsevne end<br />

samme lag i model 1. Det er formodentlig i disse kalklag hoveddelen af strømningen<br />

foregår og jo lavere den hydrauliske ledningsevne er for disse lag, jo længere er<br />

responstiden.<br />

4.6.4 Anden placering af indvindingsboringerne<br />

Indvindingsboringerne blev placeret efter potentialekortet, langs med potentialelinjerne,<br />

i henhold til afsnit 4.4 (Placering af indvindingsboringerne). Efter der blev indlagt<br />

partikelbaner i modellerne er det imidlertid blevet klart, at placeringen virker mod<br />

89


Kapitel 4 Opsætning af grundvandsmodellen<br />

hensigten, idet boringsplaceringen ligger vinkelret på potentialekurverne og således<br />

forstærker sænkningen ved hver boring. Dette er årsag til den smalle oplandsform. I det<br />

følgende vil placeringen af boringerne derfor blive ændret.<br />

Figur 4.10 vises, hvordan boringerne kan placeres, således, at de står langs med<br />

potentiale kurverne og det ses, at oplandet for begge modeller er blevet bredere med den<br />

nye boringsplacering og at sænkningen ikke længere forstærkes. Den nye placering af<br />

boringerne ligger inden for 300 m fra plantagens kant, men stadig mere end 300 m fra<br />

bebyggelsen. Det må i forhold til denne nye placering af boringer forventes, at<br />

sænkningstragten som boringerne skaber ikke bliver så dyb. Forskellen i bredden af<br />

oplandet mellem de to modeller er mindre med denne placering af indvindingsboringerne.<br />

I forhold til oplandet for de to grundvandsmodeller er det interessant, at begge modeller<br />

med begge placeringer af indvindingsboringer beregner et opland, der hovedsageligt<br />

ligger uden for Volsted Plantage. Af denne grund kan det være nødvendigt, at der<br />

iværksættes en beskyttelse af dette område mod nedsivning af pesticider og nitrat, da<br />

arealanvendelsen her er landbrug, der anvender disse stoffer.<br />

Figur 4.10. Projektområdet med de beregnede potentialer, oplandene til Hølbækken og<br />

indvindingsboringerne for den nye placering af boringerne. Indvindingen er sat til 1.000.000 m 3 /år.<br />

90


I forhold til vandføringen sker der ikke så store ændringer, hvis boringernes placering<br />

ændres. Reduktionen af vandføringen bliver 2 procentpoint lavere i model 1, mens den<br />

for model 2 er 0,3 procentpoint lavere, ved en oppumpning på 1.000.000 m 3 /år, tabel<br />

4.5. Grunden til, at der ikke er større forskel i de beregnede vandføringer er<br />

sandsynligvis, at boringernes oplande stadig begrænser Hølbækkens opland med det<br />

samme areal, da det er den samme vandmængde indvindingerne fjerner fra dette opland.<br />

Det vurderes imidlertid, at den nye placering af indvindingsboringerne vil være bedre,<br />

da sænkningstragten formodentlig ikke bliver så dyb. Resultaterne for reduktionen af<br />

vandføringen ved forskellige pumpningsrater, der tidligere er blevet præsenteret, vil<br />

imidlertid stadig blive anvendt, idet afvigelsen i vandføringen efter ændringen af<br />

placeringen af boringerne er minimal.<br />

Tabel 4.5. Vandføring beregnet ved de to placeringer af indvindingsboringerne ved indvinding på<br />

1.000.000 m 3 /år.<br />

Vandføring<br />

(m 3 /døgn)<br />

Model 1 Model 2<br />

Reduktion af<br />

vandføring (%)<br />

Vandføring<br />

(m 3 /døgn)<br />

Reduktion af<br />

vandføring (%)<br />

Ny placering 3.105,0 17,7 2.462,8 19,2<br />

Tidligere placering 3.034,0 19,7 2.470,8 18,9<br />

4.6.5 Opsummering<br />

Det er kommet frem i dette afsnit, at en vandindvinding på 1.000.000 m 3 /år påvirker<br />

Hølbækken for kraftigt i forhold til dens B1 målsætning. Ved en halvering af<br />

indvindingen til ca. 500.000 m 3 /år ligger påvirkningen tættere på det acceptable niveau.<br />

Ved undersøgelse af to forskellige placeringer af indvindingsboringerne i Volsted<br />

Plantage, har det vist sig, at placeringen ikke har den store indflydelse på Hølbækkens<br />

vandføring.<br />

4.7 Usikkerheder ved grundvandsmodellering<br />

I det følgende vil de usikkerheder, der har været i forbindelse med modellering af<br />

grundvandsstrømningerne, blive gennemgået. Først beskrives de generelle usikkerheder,<br />

der er opstået som følge af nødvendige simplificeringer og detaljeringsgraden af de<br />

tilgængelige data og det vil i den forbindelse blive gennemgået, hvad der er gjort for at<br />

sikre, at de nævnte usikkerheder ikke i så høj grad er kommet til at påvirke de opsatte<br />

modeller. Derefter præsenteres en følsomhedsanalyse på de to modeller, for at<br />

91


Kapitel 4 Opsætning af grundvandsmodellen<br />

undersøge deres stabilitet og betydningen af de forskellige parametre i modellerne for<br />

Hølbækkens vandføring.<br />

4.7.1 Usikkerheder i data<br />

Randbetingelser, stationære forhold og diskretisering er begreber, der kun er en del af<br />

modellerne og som ikke er en del af situationen i naturen, men er en række forsimplinger,<br />

der må foretages, for at få modellen tilpasset virkeligheden (Vahman, 2005: 23).<br />

Projektområdet er afgrænset ud fra et potentialekort og et kort over åer og dræn i<br />

området. Kortet over åer og dræn har en høj nøjagtighed, da disse ligger på overfladen<br />

og er nemme at bestemme. Nøjagtigheden af potentialekortet er mindre sikkert, da dette<br />

bygger interpolationer mellem en række punkter. Endvidere varierer potentialet over tid,<br />

og da punktmålingerne af potentialet er taget på forskellige årstider og med års<br />

mellemrum, øger dette usikkerheden på det anvendte potentialekort. Hølbækken ligger<br />

tæt på Lindenborg Å, der fungerer som en trykgrænse, og derfor kan det forventes, at<br />

potentialet har en høj nøjagtighed i dette område. Der er imidlertid risiko for at<br />

pumpningen af vand i Volsted Plantage kan ændre afgrænsningen af projektområdet,<br />

således, at området efter oppumpningen får en anden form. Dette antages imidlertid ikke<br />

at have påvirket resultatet væsentligt, fordi oplandene til både Hølbækken og<br />

indvindingsboringerne ikke går ud til grænserne af oplandet og derfor sandsynligvis<br />

ikke ændrer vandudvekslingen mellem projektområdet og det omgivende areal.<br />

Den overordnede geologiske model for området er opstillet ved tolkning af punktdata,<br />

dvs. boringer som er foretaget i området. Der er derfor ingen sikkerhed for, at geologien<br />

mellem punkterne lineært kan interpoleres, som det er blevet gjort i denne rapport.<br />

Interpolationen er imidlertid foretaget under hensyntagen til det kendskab, der er af<br />

områdets geologiske historie og er derfor et kvalificeret bud på geologien. Desværre er<br />

kendskabet til geologien i plantagen ikke er ret godt, da der kun er få boringer i dette<br />

område, hvilket mindsker nøjagtigheden i gengivelsen af dette område. Da der er god<br />

overensstemmelse mellem potentialet i de opsatte modeller og det observerede<br />

potentiale, tyder dette dog på, at tolkningen af geologien er foretaget med en<br />

nøjagtighed, der er tilfredsstillende i forhold til modellens anvendelse.<br />

For at simplificere modellen var det nødvendigt at antage, at de forskellige<br />

hydrogeologiske forekomster havde en ens hydraulisk ledningsevne i hele<br />

projektområdet. Undersøgelserne af lagenes hydrogeologiske egenskaber viste tydeligt,<br />

at dette ikke er tilfældet i virkeligheden, da der selv indenfor mindre afstande var<br />

forskelle i de hydrauliske ledningsevner. Det fremgik imidlertid også af forsøgene, at<br />

der er væsentlig større forskelle mellem den hydrauliske ledningsevne for de forskellige<br />

lag end der er indenfor det samme lag. Derfor betragtes denne simplificering ikke som<br />

en væsentlig fejlkilde. Værdierne for de hydrauliske ledningsevner blev taget fra egne<br />

92


forsøg, men det var nødvendigt at supplere med tabelværdier, hvilket også medfører en<br />

vis usikkerhed.<br />

Da det under opstillingen af den geologiske model var nødvendigt at sammenligge flere<br />

lag, betyder dette sammen med antagelsen om, at lagene har samme hydrogeologiske<br />

egenskaber over hele området, at projektområdets geologi bliver simplificeret i forhold<br />

til virkeligheden. Dette var også tydeligt under kalibreringen, hvor de hydrauliske<br />

ledningsevner for de forskellige lag blev ændret væsentligt. Dette var en nødvendig<br />

simplificering for at kunne opstille modellerne og da de hydrauliske ledningsevner blev<br />

kalibreret mod potentialet i området, formodes det, at de væsentligste<br />

grundvandsstrømninger i området er repræsenteret. De hydrauliske ledningsevner, der<br />

blev fundet optimale i kalibreringen af grundvandsmodellen, ligger også alle inden for<br />

det interval, der er angivet for lagene og som er præsenteret i kapitel 3 (Konceptuelle<br />

model af projektområdet). Dette antages derfor ikke at have ført til væsentlige fejlkilder.<br />

I forhold til de tilgængelige data, var det nødvendigt at opstille en stationær model, da<br />

årstidsvariationerne for området ikke er kendt. Dermed fremkommer der en<br />

gennemsnitsituation for hele året, men en del af formålet med projektet var at bestemme<br />

påvirkninger i forhold til en minimumssituation, og dette kan derfor føre til<br />

usikkerheder, da påvirkningen på minimumssituationen ikke er kendt. Den påvirkning,<br />

der angives i denne rapport, må derfor forventes at ligge i den lave ende af påvirkningen<br />

på medianminimumsvandføringen, da den gennemsnitlige vandføring ved en<br />

vandindvinding formodentlig ikke vil påvirkes nær så meget som medianminimumsvandføringen.<br />

4.7.2 Følsomhedssanalyse<br />

For at vurdere modellens stabilitet og betydningen af følsomheden på af de forskellige<br />

parametre i modellen, er der udført en usikkerhedsanalyse på modellen. Det er gjort ved<br />

at udvælge en række af de parametre i modellen, som vurderes til at have den største<br />

betydning. Disse parametre er blevet varieret indenfor et interval, som det er antaget<br />

svarer til den usikkerhed, der er på deres nøjagtighed. Det er endvidere valgt at måle<br />

disse parametre både før og efter en vandindvinding på 1 mio. m 3 /år for at undersøge,<br />

hvordan ændringen af dem påvirker vandføringen i Hølbækken, tabel 4.6. De valgte<br />

parametre er de hydrauliske ledningsevner på tørv, gytje, sand og de fire kalklag samt<br />

nettonedbøren. De hydrauliske ledningsevner blev ændret med en dekade op og ned,<br />

mens nedbøren varieredes med 10 %. Den hydrauliske ledningsevne for kalk blev<br />

imidlertid sat til maks. 100 m/d i begge modeller, da dette ca. er den øvre grænse ud fra<br />

Bach et al. (2001).<br />

93


Kapitel 4 Opsætning af grundvandsmodellen<br />

Tabel 4.6. Den procentvise ændring af vandføringen i Hølbækken før og efter indvinding ved ændring af<br />

de vigtigste parametre i modellen. Forskellen angiver stigning eller fald i ændring af vandføringen.<br />

94<br />

Model 1 Model 2<br />

Ændring Før (%) Efter (%) Ændring Før (%) Efter (%)<br />

Tørv + 1 dekade 11 12 + 1 dekade 29 -1<br />

- 1 dekade -18 -19 - 1 dekade -35 -27<br />

Gytje + 1 dekade 4 5 + 1 dekade 5 43<br />

- 1 dekade 38 40 - 1 dekade 68 72<br />

Sand + 1 dekade 29 46 + 1 dekade 52 78<br />

- 1 dekade 66 60 - 1 dekade -44 -19<br />

Kalk 1 100 m/d 25 -2 100 m/d -6 3<br />

- 1 dekade 0 1 - 1 dekade 6 6<br />

Kalk 2 + 1 dekade -5 -7 100 m/d -15 -14<br />

- 1 dekade 2 3 - 1 dekade 3 10<br />

Kalk 3 + 1 dekade 12 -23 100 m/d -24 -20<br />

- 1 dekade 6 7 - 1 dekade 49 43<br />

Kalk 4 + 1 dekade -26 -19 + 1 dekade 67 -6<br />

- 1 dekade 34 31 - 1 dekade 32 2<br />

Nedbør + 10 % 18 15 + 10 % 16 22<br />

-10 % -12 -13 -10 % 6 10<br />

Ud fra tabel 4.6 er det muligt at tolke følsomheden på de enkelte parametre. En stor<br />

afvigelse betyder, at parameteren er vigtig i forhold til simuleringen af vandføringen i<br />

Hølbækken. Det skal tages i betragtning af, at nogle lag afhænger af hinanden, så når de<br />

bliver kalibreret sammen, passer de fint i modellen, men hvis kun et af dem ændres, så<br />

kommer der en stor usikkerhed. Dette kommer af, at de andre parametre i modellen ikke<br />

er blevet genkalibreret i forhold til ændringen af parametrene. Det beregnede potentiale<br />

afviger derfor i modellerne mere i forhold til det observerede, når de nævnte parametre<br />

ændres.<br />

Det ses i tabel 4.6, at der er en stor forskel på parametrenes usikkerhed i de to modeller.<br />

Afvigelserne i model 2 er større end i model 1. Model 1 er derfor mere stabil i forhold<br />

ændringen i vandføringen, ved ændring af de forskellige lags hydrauliske ledningsevner.<br />

Den største usikkerhed forekommer ved det nederste kalklag, hvor en hævning af<br />

den hydrauliske ledningsevne giver relativt store afvigelse i model 2 og lidt mindre i<br />

model 1. Den store vandføringsændring opstår sandsynligvis, fordi laget er meget tykt<br />

og i udgangspunktet er sat til en lav hydraulisk ledningsevne. At dette lag i<br />

virkeligheden skulle kunne være en dekade større, virker derfor urimeligt, hvilket også<br />

understøttes af den flowlog, der er optaget for området. Hvis den hydrauliske<br />

ledningsevne bliver sænket i dette lag, har det til gengæld kun lille betydning, hvilket<br />

formodentlig kommer af, at den hydrauliske ledningsevne i forvejen er lav i dette lag.


Nettonedbøren har en vis betydning for usikkerheden i model 2, men mindre betydning<br />

i model 1. Følsomheden på nettonedbøren giver ikke nær så stor afvigelse som de andre<br />

parametre, hvorfor usikkerheden ved data for nettonedbør ikke anses som et problem.<br />

Sandet giver den største usikkerhed, og dette kan ikke umiddelbart forklares. Det skal<br />

imidlertid nævnes, at når den hydrauliske ledningsevne for sandet ændres, stemmer det<br />

observerede potentiale dårligt overens med det beregnede. Dette kan skyldes, at<br />

grundvandsspejlet i området omkring Hølbækken ligger i sandet. Det var her<br />

kalkoverfladen dykkede ned. Dermed har sandet en stor betydning for, hvor meget vand<br />

der når ud i Hølbækken og strømningsforholdene gennem sandet kan dermed påvirke<br />

resten af projektområdet.<br />

4.7.3 Opsummering<br />

I dette afsnit er de usikkerheder, modellen er behæftet med, blevet gennemgået. Der er<br />

blevet argumenteret for, at de generelle fejlkilder er minimeret ved at gennemføre en<br />

grundig kalibrering af de to grundvandsmodeller. Herved er det sikret, at det beregnede<br />

potentiale stemmer godt overens med det observerede. Følsomhedsanalysen viser, at<br />

model 1 er den mest stabile i forhold til ændringerne i de enkelte lags hydrauliske<br />

ledningsevner. Følsomhedsanalysen viser samtidig, at der ikke er store usikkerheder i<br />

modellen, idet de usikkerheder, der i følsomhedsanalysen bliver beregnet til at være<br />

store, alle kan forklares. Grundvandsmodellen anses derfor som repræsentativ for<br />

grundvandsforholdene i projektområdet.<br />

95


5<br />

Konklusion<br />

<strong>Rapporten</strong>s problemformulering var: Hvilken indflydelse vil oprettelsen af<br />

indvindingsboringer i Volsted Plantage få for Hølbækkens vandføring? Dette er blevet<br />

undersøgt ved hjælp af en numerisk grundvandsmodel opstillet for projektområdet. For<br />

at opstille grundvandsmodellen blev projektområdets geologi, hydrologi og<br />

hydrogeologi undersøgt, både ved hjælp af egne undersøgelser og eksisterende data. Det<br />

kom frem, at hoveddelen af grundvandsmagasinet i projektområdet ligger i et kalklag,<br />

der overlejres af en moræne. Hølbækken og Lindenborg Å løber gennem ådale, dannet<br />

ved indsynkning og erosion af kalken, opfyldningen i ådalene udgøres af et gytjelag<br />

med en lav hydrauliske ledningsevne overlejret af et tørvelag, med en højere hydraulisk<br />

ledningsevne. Undersøgelser af Hølbækken viser, at denne i høj grad er grundvandsfødt<br />

og formodentlig har en gennemsnitlig vandføring på ca. 5.100 m 3 /d.<br />

Grundvandsmodellen blev opstillet på baggrund af dette og kalibreret mod pejlinger af<br />

grundvandsspejlet. Der blev opstillet to grundvandsmodeller med forskellige tolkninger<br />

af kalken. Begge modeller viser, at Hølbækkens vandføring gennemsnitligt påvirkes<br />

med ca. 20 %, hvis der indvindes 1.000.000 m 3 /år, som planlagt. Ved ændring af<br />

indvindingsboringernes placering er påvirkningen tilnærmelsesvis den samme. Det vil<br />

ikke være muligt at gennemføre grundvandsindvindingen i Volsted Plantage uden at<br />

påvirke Hølbækkens medianminimumsvandføring med over 10 %. Indvindes der i<br />

stedet 500.000 m 3 /år, påvirkes gennemsnitsvandføringen kun med 10 %.


6<br />

Perspektivering<br />

I denne rapport er der udarbejdet en numerisk grundvandsmodel over Hølbækken, der er<br />

anvendt til at vurdere, hvilken indvirkning en vandindvinding i Volsted Plantage vil<br />

have på Hølbækken. Der har været begrænset tid og ressourcer til rådighed til<br />

udarbejdelsen af denne rapport og derfor er der en række ting, det kunne være<br />

interessant at undersøges yderligere, hvis der var mere tid og flere ressourcer.<br />

Nøjagtigheden af grundvandsmodellen kunne sandsynligvis forbedres, hvis dens<br />

datagrundlag blev udbygget. En af de store usikkerheder i opstillingen af<br />

grundvandsmodellen har været afgrænsningen af kalklaget i området og bestemmelse af<br />

dennes hydrauliske ledningsevne, hvilket i høj grad skyldes manglende kendskab til<br />

kalkens sprækkedannelse. Ved at udføre f.eks. en række flowlog forsøg fordelt over<br />

projektområdet, ville der fremkomme et bedre bud på den vertikale inddeling af kalken,<br />

som muligvis kunne interpoleres til resten af området.. En anden stor usikkerhed i<br />

dataene er variationen i Hølbækkens vandføring. For at bestemme medianminimumsvandføringen<br />

er det nødvendigt at kende variationen over flere år, men det var ikke<br />

mulig at skaffe disse data. I stedet blev der anvendt data fra Lindenborg Å, som<br />

projektområdet er en del af. Et mere nøjagtigt resultat ville fås, ved at måle<br />

vandføringen i selve Hølbækken. Udover disse to store usikkerheder, er der mulighed<br />

for at forbedre nøjagtigheden af de fleste af de anvendte data.<br />

Det kunne også være interessant at arbejde videre med resultatet af grundvandsmodellen<br />

og vurdere, hvad konsekvenserne heraf bliver. Modellen viste, at en vandindvinding på<br />

1.000.000 m 3 /år påvirker Hølbækken for kraftigt i forhold til dens B1 målsætning. En<br />

halvering af vandindvindingen ville give en påvirkning, der er tættere på et acceptabelt


Kapitel 6 Perspektivering<br />

niveau, men det er ikke sikkert, at det er økonomisk rentabelt at lave en vandindvinding<br />

i denne størrelse. Der er mulighed for at anvende afbødende foranstaltninger, som f.eks.<br />

pumpning af vand direkte til Hølbækken, og det kunne også være en mulighed at<br />

placere indvindingsboringerne anderledes for at mindske påvirkningen. En sidste udvej<br />

kunne være at slække på målsætningen for Hølbækken, da drikkevandsindvinding har<br />

højere prioritet end miljøtilstanden, men dette kan kun lade sig gøre, hvis der ikke er<br />

andre muligheder for at opfylde drikkevandsbehovet (Regionplan, 2005: 236).<br />

Udover påvirkningen af vandføringen i Hølbækken, kan vandindvinding i Volsted<br />

Plantage, som nævnt i indledningen, også påvirke områder af økologisk interesse, da<br />

dele af Hølbækkens ådal og Lindenborg Ås ådal er habitatområder i EUs<br />

Habitatdirektiv. En sænkning af grundvandsspejlet kan medføre ændringer af naturen,<br />

eksempelvis frigivelse af okker og ændrede vilkår for vegetationen i området. Der er<br />

desuden en række økonomiske konsekvenser ved en grundvandsindvinding i Volsted<br />

Plantage, idet både det nærliggende dambrug og lodsejerne i området kan blive påvirket<br />

økonomisk. Dambruget ved en påvirkning af vandkvaliteten til deres fisk og lodsejerne<br />

ved en sætning af deres jorde, hvilket kan vanskeliggøre driften. For at vurdere<br />

omkostningerne ved vandindvinding i Volsted Plantage kunne det derfor være<br />

fordelagtigt at undersøge disse.<br />

100


Litteraturliste<br />

AK, 2005<br />

Indvindingsdata, 2006, http://www.aalborg.dk/NR/rdonlyres/62409785-5FAA-4ABC-BA20-<br />

EB0FC9B2E302/4939/Bilag2bskema.pdf, Aalborg Kommune<br />

AK, 2006<br />

Indvindingsdata, 2006, http://www.aalborg.dk/borgerportal/applikationer/kloakspildevand/spildevandsplan/AKoversigtskort.pdf,<br />

Aalborg kommune<br />

Bach et al., 2001<br />

Bach, J; Rikskov, U., Gravesen, P., 2001, Udvidet geologi og Grundvand, Miljøstyrelsen<br />

Bennedsen, 2002<br />

Bennedsen, R., 2002, Kvalitetsplan for vandløb og søer, Nordjyllands Amt, Teknik og miljø<br />

Berthelsen, 1987<br />

Berthelsen, O., 1987, Geologi i Aalborgområdet, råstoffer – fundering – vandindvinding,<br />

Danmarks Geologiske Undersøgelse (DGU)<br />

Brun, 2005<br />

Brun, A., 2005, Diskretisering af modelområdet i til og sted, i: Sonnenborg, T., Henriksen, H. J.,<br />

2005, Håndbog i grundvandsmodellering, GEUS<br />

Burcharth & Jørgensen, 1976<br />

Burcharth, H. F. & Jørgensen, T. W., 1976, Hydrologi, Laboratoriet for hydraulic og<br />

havnebygning Aalborg Universitetscenter<br />

Burt et al., 2002<br />

Burt, T.P., Butturini, A., Clement, J.C., Danielescu, S., Dowrick, D.J., Haycock, N.E., Hefting,<br />

M.M., Hillbricht-Ilkowska, A., Maitre, V., Matheson, F.E., Pinay, G., 2002, Water table<br />

fluctuations in the riparian zone: comparative from a pan-European experiment Journal of<br />

Hydrology 265 129–148, ELSEVIER<br />

Dahm et al., 1998<br />

Dahm, C., Grimm, N., Marmonier, P., Valett, H., Vervier, P., Nutrient dynamics at the interface<br />

between surface waters and groundwaters, Freshwater Biology 1998 nr. 4 : 427-451, Blackwell<br />

Science Ltd.<br />

Danva, 2006<br />

Dansk vand- og spildevandsforening, 2006, http://www.danva.dk/sw160.asp<br />

DDH, 1995<br />

Danmarks Digitale Højdemodel, 1995, metadata på http://www.geodata-info.dk/ds.asp?-<br />

DS=341, KMS<br />

DJF-geodata, 2006<br />

Jordbundskort, Metadata på www.djf-geodata.dk, DJF<br />

DMI 1, 2006<br />

Vejrobservationer, 2006, http://www.dmi.dk/dmi/vejrarkiv?region=1&year=2006&month=3,<br />

Danmarks Meteorologiske Institut<br />

DMI 2, 2006<br />

Vejrobservationer, 2006, http://www.dmi.dk/dmi/vejrarkiv?region=1&year=2006&month=4,<br />

Danmarks Meteorologiske Institut


Litteraturliste<br />

DMU 1, 2006<br />

Ordbog, www.dmu.dk, Danmarks Miljø Undersøgelser, 2006<br />

DMU 2, 2006<br />

Vandføringsdata, http://www2.dmu.dk/1_Om_DMU/2_tvaer-funk/3_fdc_hyd/hyd_db/nordq.asp,<br />

Danmarks Miljø Undersøgelser, 2006<br />

Fetter, 2001<br />

Fetter, C. W., 2001, Applied Hydrogeology, fourth edition, Prentice-Hall, 2001<br />

Finnern et al., 1996<br />

Finnern, H., Grottenthaler, W., Kühn, D., Pälchen, W., Schraps, W. G., Sponagel, H., 1996,<br />

Bodenkundliche Kartieranleitung, Bundesanstalt für Geowissenschaften und Rohstoffe und den<br />

Geologischen Landesämtern in der Bundesrepublik Deutschland<br />

Fobian, 1984<br />

Fobian, A., 1984, Jordbundskortlægning og arealklassificering, Instituttet for teknisk geologi,<br />

Danmarks Tekniske Højskole<br />

Fredericia, 1989<br />

Fredericia, J., 1989, Den hydrogeologiske kortlægning af Nordjyllands Amt, Danmarks<br />

Geologiske Undersøgelse (DGU)<br />

Galsgaard, 1998<br />

Galsgaard, J., 1998, Indføring i sedimentologi, Dansk geoteknisk forening<br />

GEUS 1, 2006<br />

Data tilsendt af, Kelstrup, N., GEUS, 18.05.2006<br />

GEUS 2, 2006<br />

Boringer fra GEUS’s Jupiter database<br />

GEUS, 1994<br />

Prækvartæroverfladens højdeforhold, 1994, metadata på http://www.geodata-info.dk/ds.asp?-<br />

DS=46, GEUS<br />

Hedeselskabet, 2006<br />

Borerapport fra Volsted Plantage, 2004, Borehulslogging og prøvebeskrivelse ved Volsted, nr.<br />

34.2698, NJA<br />

Henriksen et al., 2001<br />

Henriksen, H. J., Sonnenborg, T., Christiansen, H. B., Refsgaard, J. C., Harrar, B., <strong>Rasmussen</strong>,<br />

P., Brun, A., 2001, Retningslinier for opstilling af grundvandsmodel, GEUS og DHI<br />

Henriksen, 2005<br />

Henriksen, H. J., 2005, Dynamisk model, i: Sonnenborg, T., Henriksen, H. J., 2005, Håndbog i<br />

grundvandsmodellering, GEUS<br />

Jensen, 2002<br />

Jensen, J. B., 2002, Parameter and Uncertainty Estimation in Groundwater Modelling,<br />

Hydraulics & Costal Engineering Laboratory, Aalborg University<br />

Jensen, 2005<br />

Jensen, K. H., Matematisk og numerisk grundlag for grundvandsmodellering, i: Sonnenborg, T.,<br />

Henriksen, H. J., 2005, Håndbog i grundvandsmodellering, GEUS<br />

102


Jensen et al., 2005<br />

Jensen, J. B., Langhoff, J. H., Grønvald, P., Mortensen, M., Bliksted, T., Improving model<br />

results on a river valley scale using GSI techniques, NIRAS, Aalborg Utility Company<br />

Langhoff et al., 2005<br />

Langhoff, J. H., Christensen, S., <strong>Rasmussen</strong>, K. R., 2005, Quantification and regionalization of<br />

groundwater–surface water interaction along an alluvial stream Journal of Hydrology xx (2005)<br />

1–17, ELSEVIER<br />

Larkin & Sharp, 1992<br />

Larkin, R. & Sharp, J., 1992, On the relationship between river-basin geomorfology aquifer<br />

hydraulics and ground-water flow direction in alluvial aquifers, Geological Society of America<br />

Bulletin, vol. 104, Side 1608-1620, Geological Society of America<br />

Larsen et al., 1988<br />

Larsen, G., Frederiksen, J., Villumsen, A., Fredericia, J., Gravesen, P., Foged, N., Knudsen, B.,<br />

Bauman, J., 1988, vejledning i Ingeniørgeologisk prøvebeskrivelse, DGF bulletin nr. 1, Dansk<br />

Geoteknisk forening<br />

Larsen, 1989<br />

Larsen, G., 1989, Træk af Danmarks geologi, Dansk geoteknisk forening<br />

Miljøbeskyttelsesloven, 2006<br />

Miljøbeskyttelsesloven, 2006, http://www.retsinfo.dk/_GETDOC_/ACCN/A20010075329-regl,<br />

Miljøministeriet<br />

Miljøstyrelsen, 2000<br />

Forfatter, 2000, Zonering, Detailkortlægning af arealer til beskyttelse af Grundvandsressourcen<br />

Vejledning fra miljøstyrelsen Nr. 3, Miljøstyrelsen<br />

Naturbeskyttelsesloven, 2004<br />

Naturbeskyttelsesloven, 2004, http://www.retsinfo.dk/_GETDOCM_/ACCN/A20040088429-<br />

REGL, Skov og Naturstyrelsen<br />

Nilsson et al., 2003<br />

Nilsson, B., Andersen, H. E., Christensen, S., Dahl, M., Hofmann, C. C., Kronvang, B.,<br />

Langhoff, J. H., Møller, I., I. <strong>Rasmussen</strong>, I., Refsgaard, J. C., 2003, HYdrokemisk interaktion<br />

mellem GRundvand og Overfladevand (HYGRO), Arbejdsrapport fra Miljøstyrelsen Nr. 10,<br />

Miljøstyrelsen<br />

Nilsson et al., 2004<br />

Nilsson, B., Andersen, H. E., Christensen, S., Dahl, M., Hofmann, C. C., Kronvang, B.,<br />

Langhoff, J. H., Møller, I., <strong>Rasmussen</strong>, K. R., Refsgaard, J. C., 2004, Videreudvikling af<br />

ådalstypolog - Grundvand- Overfladevand Interaktion (GOI), Arbejdsrapport fra Miljøstyrelsen<br />

Nr. 16, Miljøstyrelsen<br />

Nilsson et al., 2005<br />

Nilsson, B., Banke, M., Christensen, S., Dahl, M., Engesgaard, P., Fuglsang, A., Kronvang, B.,<br />

Ovesen, N. B., Platen-Hallermund, F. V., Sonnenborg, T., Tornbjerg, H., Wohlfeil-Müller, D.,<br />

2005, Afslutning af ådalstypologi – (Grundvand-Overfladevand Interaktion), Arbejdsrapport fra<br />

Miljøstyrelsen Nr. 17, Miljøstyrelsen<br />

NIRAS, 2005<br />

Skræntfod, 2005, NIRAS<br />

NJA, 2005<br />

Potentiale data, 2005, Nordjyllands Amt<br />

103


Litteraturliste<br />

NJA, 2006<br />

Historiske kort, 2006, http://arealinfo.nja.dk, Nordjyllands Amt<br />

NJA, 2003<br />

Infiltrationskort, 2003, Regional kortlægning af infiltration og sårbarhed på grundlag af<br />

vandbalance for den umættede zone, Nordjyllands Amt<br />

<strong>Rasmussen</strong>, 1975<br />

<strong>Rasmussen</strong>, H. W., 1975, Danmarks Geologi, Gjellerups liniebøger<br />

Refsgaard & Henriksen, 2005<br />

Refsgaard, J. C. & Henriksen, H. J., 2005, Opstilling af hydrogeologisk tolkningsmodel, i:<br />

Sonnenborg, T., Henriksen, H. J., 2005, Håndbog i grundvandsmodellering, GEUS<br />

Refsgaard, 2005<br />

Refsgaard, J. C., 2005, Indledning, i: Sonnenborg, T., Henriksen, H. J., 2005, Håndbog i<br />

grundvandsmodellering, GEUS<br />

Regionplan, 2005<br />

Nordjyllands amtsråd, 2005, Regionplan 2005, Nordjyllands amt – Teknik- og Miljørådet<br />

SNS, 2006<br />

Skov og Naturstyrelsen, 2006,<br />

http://www.skovognatur.dk/Emne/Natura2000/Habitat/Habitatdirektivet/ og<br />

http://www.skovognatur.dk/Emne/Naturbeskyttelse/International/Vandrammedirektiv/<br />

Sonderegger & Weight, 2000<br />

Sonderegger, J. L. & Weight, W. D., 2000, Manual of applied field hydrogeology, McGraw-Hill<br />

Spitz & Moreno, 1996<br />

Spitz, K & Moreno, J., 1996, A Practical Guide to Groundwater and Solute Transport Modeling,<br />

Wiley & Sons<br />

Statistikbanken, 2006<br />

Statistiske data, 2006, www.dst.dk, Danmarks Statistik<br />

Søndergaard, 2003<br />

Søndergaard, K., 2003, Bestemmelse af indvindingsopland til Kongshøj Kildeplads, Aalborg<br />

Universitet<br />

Thomsen, E., 1995<br />

Thomsen, E., 1995, Kalk og kridt i den danske undergrund, i: Nielsen, O. B., 1995, Danmarks<br />

geologi fra Kridt til i dag, Geologisk institut, Århus Universitet<br />

TOP10DK, 2000<br />

Data fra TOP10DK, 2000, Metadata på http://www.geodata-info.dk/ds.asp?DS=103, KMS<br />

Vahman, 2005<br />

Vahman, M., 2005, Modelusikkerhed eller kejserens nye klæder, Vand & Jord, vol 12, 2005.<br />

Wikipedia, 2006<br />

Online leksikon, 2006, www.wikepedia.com<br />

Watertech, 2005<br />

Notat vedrørende TEM-måling i Volsted Plantage, 2005, Watertech<br />

104


Youngs, 2000<br />

Youngs, E. G., 2000, Hydraulic Conductivity of Saturated Soils, I, Smith, Keith A., Soil and<br />

Environmental Analysis: Physical Methods, Marcel Dekker Incorporated<br />

105


Bilag A<br />

Flowlog<br />

Figuren viser den flowlog, der blev lavet den 10.11.04, af Hedeselskaber for<br />

Nordjyllands Amt. Flowloggen blev lavet 1 km nordvest for Hølbækkens udspring. Det<br />

ses af flowloggen, at vandudstrømningen fra lagene stiger langsomt fra grundvandsspejlet,<br />

10 m. u. t. til 30 m. u. t., fra 30-80 m er udstrømningen omtrent den samme, på<br />

nær et lille udsving ved 50 m. I en dybde på 80-90 m svinger udstrømningen meget og<br />

bliver helt ubetydelig under 90 m (Hedeselskabet, 2006).<br />

Bilag A


B1 Vandføringsmålinger<br />

Bilag B<br />

Vandføring<br />

For at bestemme, hvilke områder af Hølbækken, der bidrager til vandføringen og for at<br />

finde den samlede vandføring, er der udført vandføringsmålinger en række steder langs<br />

Hølbækken. Det er vigtigt, at målingerne udføres kort tid efter hinanden, helst samme<br />

dag, da vandføringen varierer over tid, og dette kan forstyrre sammenligningen af de<br />

forskellige målepunkter. Vandføringsmålingerne i Hølbækken blev udført den 27.02 og<br />

den 28.02, mens to af målepunkterne blev målt igen den 20.04 til sammenligning.<br />

Vejret op til målingerne kan ses i tabel B1 og i figur B1.<br />

Vandføringen et bestemt sted i et vandløb kan bestemmes ved at multiplicere arealet af<br />

vandløbets tværsnit med vandets gennemsnitshastighed. Vandets gennemsnitshastighed<br />

findes ved at opdele tværsnittet i et kvadratnet og måle vandhastigheden i hvert af<br />

nettets kvadrater. Vandføringen (Q) bestemmes af ligning 1:<br />

Ligning 1<br />

Hvor:<br />

Q = Vandføring (cm 3 )<br />

Ai = Delareal i (cm 2 )<br />

Vi = Delarealets middelhastighed (cm/s)<br />

Q <br />

n<br />

<br />

i1<br />

Herved fremkommer et hastighedsprofil, der viser, hvor hurtigt vandet strømmer i de<br />

forskellige kvadrater i tværsnittet. Ved at beregne vandføringen i alle nettets kvadrater<br />

og summere resultaterne, findes den samlede vandføring for tværsnittet. Nøjagtigheden<br />

af målingen afhænger af det antal kvadrater, som tværsnittet opdeles i. Et større antal<br />

kvadrater giver en højere nøjagtighed.<br />

For at forbedre nøjagtigheden af måleresultatet, og simplificere beregningerne, kan det<br />

være en fordel at afgrænse vandløbet i siderne med plader, således at der opnås et<br />

tværsnit, der er tilnærmet rektangulært. Endvidere er der en række forhold, der skal<br />

tages i betragtning, når der vælges et målested. Det bør være en retlinet strækning med<br />

en jævn vandstrømning og en jævn bund uden store sten, der kan forårsage turbulens.<br />

Der bør heller ikke være grøde tilstede, som kan påvirke vandføringen. For at mindske<br />

påvirkningen af vandstrømningen under prøvetagningen, måles vandhastigheden ved<br />

brug af en lille propel, der sænkes ned i bækken for enden af en tynd stang.<br />

AiVi<br />

Bilag B


Fejlkilder<br />

Der indføres en vis usikkerhed ved at antage, at hastigheden er ens i hele det anvendte<br />

kvadrat. Vandløbets udformning giver også anledning til usikkerheder, da det ikke er<br />

muligt at bestemme arealet af kvadraterne præcist.<br />

Vurdering<br />

Det er vurderet, at vandføringsmålingerne giver et fornuftigt billede på den reelle<br />

vandføring i Hølbækken, når det beregnede konfidensinterval tages i betragtning.<br />

Tabel B1. Oversigt over det registrerede vejr op til vandføringsmålingerne (Nordjyske, april, 2006; DMI<br />

1, 2006).<br />

Dato Antal vandføringsmålinger<br />

udført<br />

Vejr op til måling (ca. en uge)<br />

27.02 6 (nr. 1, 2, 6, 7, 8, 9) Dagtemperaturer 0,5 – 4 °C, nattefrost.<br />

Noget nedbør i form af sne, ophold på undersøgelsesdagen.<br />

28.02 3 (nr. 3, 4, 5) Noget sne på undersøgelsesdagen, dog 2 timers ophold op til<br />

måling.<br />

20.04 2 (nr. 6, 7) Dagtemperaturer på 8 – 12 °C, ingen nattefrost.<br />

Ingen nedbør på undersøgelsesdagen, nogen nedbør i<br />

dagene op til felttur, med maksimum den 19.04 med 11 mm.<br />

Figur B1. Øverst tv. Nedbør marts 2006. Øverst th. Tryk marts 2006. Nederst tv. Nedbør april 2006.<br />

Nederst april 2006. Alle gælder for Nordjylland (DMI 2, 2006).


B2 Vandføringsstandsmålinger<br />

For at undersøge fluktuationer i vandføringen i Hølbækken, blev der anvendt to divere.<br />

Disse blev sat til at registrere trykniveauet engang i timen over en periode på 52 dage<br />

mellem den 27.02 og den 20.04. Den ene blev lagt under vandoverfladen i Hølbækken<br />

og registrerede trykniveauet her, mens den anden blev lagt på bredden ved atmosfærisk<br />

tryk. Ved at sammenholde disse to målinger findes de forskelle i vandets trykniveau,<br />

som ikke skyldes atmosfæren og dermed er et udtryk for fluktuationer i vandføringen.<br />

Fejlkilder<br />

Mellem den 5.03 til 18.03 er der registreret højere tryk i atmosfæren end i vandløbet.<br />

Ud fra dataene kan det se ud til, at kombinationen af højt tryk og ingen nedbør har<br />

medført et fald i vandføringen i en sådan grad, at diveren er kommet over vandspejlet.<br />

Den 14.03 er der 5 timers periode, hvor trykniveauet i vandløbet er målt til op til 160<br />

cm (ca. 60 cm over atmosfærisk tryk). Udsvinget tilskrives målefejl, da diveren havde<br />

registreret et fald i vandets temperatur til – 8 °C på samme tidsinterval, samtidigt med at<br />

temperaturen målt i luften var ca. -1 °C og der var ikke registreret nogen nedbør. Det<br />

målte trykniveau følger trykmålinger for Nordjylland, taget ved Aalborg Lufthavn,<br />

figur B1.<br />

Vurdering<br />

Resultaterne fra divermålingerne af vandføringen indikerer, at måleudstyret har<br />

foretaget enkelte fejlmålinger, men da disse skiller sig kraftigt fra de øvrige målinger og<br />

resultaterne fra den anden diver, er det nemt at udskille disse. Derfor er trykniveauet i<br />

Hølbækken målt med en høj nøjagtighed, men desværre var det kun muligt at måle over<br />

en kort periode.<br />

Bilag B


Bilag C<br />

Fasebestemmelse<br />

Alle intaktprøverne blev udtaget under grundvandsspejlet. Formålet med forsøget var at<br />

bestemme indholdet af vand, partikler og organisk materiale i intaktprøver fra<br />

feltområdet, både for at finde interne forskelle mellem prøver fra samme lokalitet og for<br />

at finde forskelle mellem de forskellige lag. Endvidere var ønsket at finde fraktionen af<br />

partikler mindre end 0,125 mm, især organiske materiale under denne størrelse, da det<br />

var forventet, at dette ville sige noget om de interne forskelle mellem den hydrauliske<br />

ledningsevne fundet ved samme lokalitet i tillæg til, at denne inddeling ville sige hvor<br />

omsat det organiske materiale var.<br />

Udførelse<br />

Først blev volumen og startvægten af prøverne bestemt, hvorefter de blev spredt ud på<br />

en bakke og lagt til lufttørring i 3 døgn og til sidst vejet igen. Herefter blev der udført<br />

sigteanalyse på prøverne. Tørven kunne ikke sigtes, da den hovedsageligt bestod af<br />

fibre. Gytje- og sandprøverne blev mortet og sigtet, med sigter i størrelser: 8 mm, 2 mm<br />

og 0,125 mm. 8 mm sigten blev sat som grænsen for grus, mens materialer under 2 mm<br />

klassificeres som sand og grus og materiale under 0,125 mm klassificeres som silt og<br />

ler. Efter sigtningen blev prøverne ovntørret, for at bestemme indholdet af det<br />

resterende vand. Til sidst blev prøverne glødet for at bestemme indholdet af organiske<br />

materiale, hvor fraktionerne fra sigtningen blev holdt adskilt.<br />

Fejlkilder<br />

Ved tørring af prøverne kan de mindste korn i silt- og lerfraktionerne klumpe sammen<br />

og danne større partikler, hvilket giver en forkert kornstørrelsesfordelingen ved<br />

sigteanalysen. Dette blev forsøgt undgået ved kun at lufttørre prøverne inden sigtningen,<br />

så der stadig var vand tilstede til at forhindre sammenklumpning. Da der alligevel<br />

forekom sammenklumpning blev prøverne mortet, i et forsøg på at nedbryde de store<br />

klumper til deres naturlige størrelse. På trods af mortningen blev det erfaret, at der<br />

stadig var en sammenklumpning af de små partikler, hvilket har betydet, at fraktionen af<br />

materiale mindre end 0,125 mm i prøverne er blevet undervurderet. Mortningen har<br />

desuden ført til, at noget af det organiske materiale blev knust til mindre dele, hvilket<br />

besværliggør klassificeringen i forhold til hvor nedbrudt det organiske materiale var.<br />

Vurdering<br />

Resultaterne er behæftet med en vis usikkerhed, hvor fraktionen af små partikler er<br />

undervurderet i alle prøverne. Resultaterne anvendes imidlertid ikke direkte, men som<br />

forholdsvise fraktioner mellem forskellige jordtyper. Ved tolkningen af resultaterne<br />

bliver fejlkilderne således undgået.<br />

Bilag C


Bilag D<br />

Hydraulisk ledningsevne<br />

Udførelse<br />

For at bestemme den hydrauliske ledningsevne i de forskellige lag, der er repræsenteret<br />

i Hølbækkens opland, blev der udført gennemstrømningsforsøg på intaktprøver udtaget<br />

fra 4 forskellige lag. Prøverne blev tætnet med filtre, som vandet kunne trænge<br />

igennem, men som stoppede partiklerne. Herefter blev der tilsat vand med et konstant<br />

trykniveau i toppen af prøven og det blev målt, hvor meget der løb igennem pr. tid. Ved<br />

hjælp af Darcys lov, ligning D1, er det muligt at beregne den hydrauliske ledningsevne<br />

i prøven.<br />

Ligning D1. Darcy’s lov til beregning af den hydrauliske ledningsevne (Youngs, 2000).<br />

Hvor:<br />

K : Den hydrauliske ledningsevne (cm/s)<br />

Q : Flowet (ml/s)<br />

L : Prøvens længde (cm)<br />

A : Prøvens tværsnitsareal (cm 2 )<br />

h : Vandsøjlens højde over prøven (cm)<br />

Q L<br />

K <br />

A<br />

h<br />

Luftbobler og lignende uregelmæssigheder kan forstyrre vandstrømmen og skabe<br />

fluktuationer. For at undgå denne fejlkilde, blev der først foretaget målinger efter at en<br />

vandmængde svarende til porevolumen var strømmet igennem. Det gennemstrømmede<br />

vand blev opfanget i en bakke, hvoraf en del fordampede. For at tage højde for denne<br />

fordampning, blev fordampningen fra en fyldt bakke uden vandtilførsel målt. I det<br />

sidste forsøg blev denne fremgangsmåde ændret, hvor udstrømning i stedet blev sat lig<br />

indstrømningen, da det antages, at gennemstrømningen på dette tidspunkt var stabilt.<br />

Forsøgsopstillingen kan ses på figur D1.<br />

Fejlkilder<br />

Ved udtagelse af intaktprøver er det næsten umuligt ikke at påvirke prøverne, især langs<br />

cylinderens rand. Derfor kan prøverne ikke betragtes om værende fuldstændigt intakte,<br />

og der kan forekomme afvigelse i prøvernes egenskaber i forhold til den naturlige<br />

tilstand. Intaktprøver er dog stadig bedre end lejrede prøver, hvor løs jord udtaget i<br />

området, forsøges pakket, så de svarer til de oprindelige forhold. Lokale variationer i<br />

lagene kan også påvirke resultaterne, hvilket er forsøgt mindsket ved at udtage tre<br />

prøver fra hvert lag.<br />

Langs randen af de cylindre, der er anvendt til forsøget, vil strømningsforholdene være<br />

anderledes end i resten af prøven. De kan således fungerer som makroporer, der lader<br />

vandet strømme hurtigere igennem end resten af prøven. For at mindske denne<br />

Bilag D


usikkerhed, er det nødvendigt at anvende større cylindre end der har været til rådighed<br />

til dette forsøg.<br />

Vurdering<br />

Dette forsøg har nogle store fejlkilder. Dog er der højst en faktor 3 mellem resultaterne<br />

fra samme lag, så selvom forsøget ikke kan bruges til at bestemme en nøjagtig værdi for<br />

lagenes hydrauliske ledningsevne, så giver de en indikation af hvilken størrelsesorden,<br />

den befinder sig i.<br />

Figur D1. Billedet viser den opstilling der er blevet anvendt til udførelsen af gennemstrømningsforsøget.


Bilag E<br />

Slugtest<br />

En slugtest udføres ved, at der først bores et hul i jorden, med en bund der ligger under<br />

grundvandsspejlet. Grundvandsspejlet hæves eller sænkes, hvorefter tiden der går, til<br />

grundvandsspejlet indstilles til den naturlige højde, måles.<br />

Der findes flere metoder til udførelse af slugtests, afhængig af om grundvandsmagasinet<br />

er frit eller spændt og alt efter hvilken type boring der skal måles på (Fetter, 2001: 193).<br />

I undersøgelsesområdet er der et frit grundvandsmagasin og uforede håndboringer blev<br />

brugt til udførelse af slugtesten. Derfor blev det valgt at anvende Hvorslev slug-test<br />

metoden som beskrevet i Fetter (2001), som metodisk grundlag for testen. At uforede<br />

håndboringer blev anvendt til undersøgelsen havde udelukkende den effekt, at den<br />

hydrauliske ledningsevne blev bestemt for hele lagfølgen og ikke kun et enkelt lag, som<br />

det kunne være, hvis boringen havde været foret. I det følgende vil fremgangsmåden for<br />

udførelsen af undersøgelsen blive gennemgået.<br />

Hullene til slugtesten blev boret første dag på feltturen og selve slugtesten blev udført<br />

anden dag, for at være sikker på, at grundvandsspejlet havde indstillet sig i det naturlige<br />

niveau. Der var en frygt for, at grundvandet ikke ville indstille sig med det samme, da<br />

materialet blev antaget at være lavpermeabelt. I det borede hul blev slugtesten udført<br />

ved, at der i hullet blev nedsænket en vandfyldt 1,5 l flaske (en slug) under<br />

grundvandsspejlet. Dette førte til, at grundvandsspejlet i borehullet steg med den<br />

mængde vand, der blev fortrængt af den nedsænkede flaske. I løbet af de næste omkring<br />

to timer blev vandstanden målt hvert sekund ved hjælp af en diver, en trykmåler, der<br />

blev placeret i bunden af det borede hul. Efter to timer havde grundvandsstanden<br />

indstillet sig til startniveauet igen og flasken blev taget op af borehullet. Dette førte til,<br />

at vandstanden i hullet faldt med omtrent samme højde som den tidligere var steget,<br />

højden på vandspejlet blev logget i løbet af den tid det tog for vandspejlet af indstille<br />

sig, herefter blev vandstanden i borehullet igen målt i to timer. Der blev derfor i alt målt<br />

to tidsserier over grundvandsspejlets tilbagevenden til naturligt leje. På baggrund af<br />

disse kurver og ligning E1 kan den hydrauliske ledningsevne beregnes.<br />

Ligning E1. Formel til beregning af hydraulisk ledningsevne (Fetter, 2001: 194).<br />

Hvor:<br />

2 L<br />

r ln<br />

r<br />

K<br />

2 Le<br />

t <br />

<br />

<br />

<br />

<br />

K : Den hydrauliske ledningsevne (cm/s)<br />

r : Radius af det borede hul (cm)<br />

Le : Dybden af det borede hul (cm)<br />

t37: Den tid det tager for vandet at falde til 37 % af det oprindelige niveau i hullet (s)<br />

e<br />

37<br />

<br />

<br />

<br />

Bilag E


For at visualisere resultaterne blev forholdet mellem den højde vandet steg til (h) og<br />

udgangshøjden (h0) som en funktion af tiden optegnet, figurerne findes på den vedlagte<br />

Cd-rom. Dette er en eksponentiel funktion og kurven blev derfor optegnet på en<br />

semilogaritmisk skala. I forbindelse med beregningen af den hydrauliske ledningsevne<br />

var de nødvendige parametre i forvejen ikke kendte. Radius af det borede hul var<br />

eksempelvis umulig at måle, den var antageligt den samme ned gennem hele hullet,<br />

fordi dette nogle steder var sunket sammen på grund af materialets løse konsistens.<br />

Radius blev derfor beregnet, hvilket var muligt, da det var kendt hvilken volumen, der<br />

var blevet sænket ned i hullet. Radius over og under vandspejlet varierede i nogle<br />

tilfælde meget, da det blev observeret, at hullet lettere faldt sammen under<br />

grundvandsspejlet. Hullets dybde var heller ikke helt den samme da slugtesten blev<br />

udført som da hullet blev gravet, også fordi hullet var faldet lidt sammen. Som et udtryk<br />

for hullets dybde anvendtes derfor den dybde diveren kom ned på. Tiden der gik til<br />

vandet var faldet til 37 % af det oprindelige niveau blev beregnet på baggrund af den<br />

graf over h/h0-kurven, der blev tegnet op i den logaritmiske skala. Fordelen ved at<br />

bestemme den hydrauliske ledningsevne to gange for hvert hul er, at der derved laves en<br />

kontrol af resultatet.<br />

Fejlkilder<br />

Som nævnt i gennemgangen af resultaterne var der visse problemer i forbindelse med<br />

udførelsen af forsøgene. I dette afsnit vil usikkerheder i denne forbindelse blive<br />

gennemgået. Det var besværligt at få boret huller til slugtesten, der ikke faldt sammen.<br />

Problemet var, at gytje- og tørvlagene nærmest er ”flydende” når de er vandmættede.<br />

Når der så bliver boret, vil de naturligt falde sammen. De borede huller blev derfor ikke<br />

fyldt af vand med det samme, hvilket havde den konsekvens, at der i bunden af hullet<br />

var en trykforskel, der ikke blev udlignet, hvorfor tørven og gytjen strømmede ud i<br />

borehullet og fyldte dette op. At hullet faldt sammen skabte også problemer i<br />

forbindelse med beregningerne, da det ikke var muligt, at anvende den radius og den<br />

dybde der var boret ned til, i beregningerne. Det blev derfor, som tidligere nævnt, valgt<br />

at udregne disse parametre i stedet, hvilket må være behæftet med en hvis usikkerhed.<br />

Det tyder ikke på, at der er usikkerheder på beregningen af radius, idet det gennem<br />

faldet i vandstanden er muligt, ret præcist, at udregne radius. Der er imidlertid mulighed<br />

for, at hullet har været dybere end den dybde, diveren sank ned til og at boringen af<br />

hullet har skabt nogle kanaler, som vandet kunne strømme igennem. Den beregnede<br />

hydrauliske ledningsevne kan derfor være lavere end beregnet i dette forsøg. For at<br />

forhindre, at borehuller falder sammen inden en slugtest udføres, kan hullerne fores med<br />

et rør og rundt omkring røret kan der fyldes op med grus. Hvis det er ønsket kun at<br />

bestemme den hydrauliske ledningsevne for et enkelt lag, kan der fyldes grus langs røret<br />

i den dybde dette lag ligger og et lag af impermeabelt materiale, fx bentonit, langs røret<br />

over og under dette lag for at stoppe vandtilførslen fra disse lag (Sonderegger &<br />

Weight, 2000: 339). Det har der imidlertid ikke været mulighed for at gøre i denne<br />

undersøgelse.<br />

De overstående problemer er blevet forsøgt imødekommet ved, at udføre to slugtests for<br />

hvert af de borede huller. Dette har givet en nyttig baggrund for vurdere, hvorvidt de<br />

udførte slugtests er brugbare i forhold til at sige noget om den hydrauliske ledningsevne<br />

for de borede huller. Det har dermed også givet mulighed for, at de mest troværdige<br />

målinger troværdige kunne udvælges til den videre undersøgelse.


Bilag F<br />

Åbundens hydrauliske ledningsevne<br />

Metode<br />

For at bestemme åbundens hydrauliske ledningsevne blev der anvendt to gennemsigtige<br />

plastikrør, der blev placeret i vandløbsbunden. Det ene rør blev fyldt delvis op med<br />

vand, mens det andet blev ikke fyldt og repræsenterede således koten på vandspejlet i<br />

ådalen på dette sted, figur F1. I løbet af dagen blev det noteret, hvor meget vandspejlet<br />

faldt og derudfra blev der opsat en tidsserie over flowet af vand ud af søjlen gennem<br />

åunden. Dette blev anvendt til at bestemme den hydrauliske ledningsevne, på samme<br />

måde som i gennemstrømningsforsøget. Det andet rør blev anvendt til at måle koten på<br />

vandspejlet, således, at dette kunne anvendes til beregning af åbundens hydrauliske<br />

ledningsevne. Udregningen af åbundens hydrauliske ledningsevne blev foretaget ved, at<br />

anvende Darcys formel, under antagelse af, at den dybde røret var presset ned i<br />

åbunden, var prøvens længde, og overhøjden på vandspejlet var det overtryk, der var på<br />

prøvens overside. Gradienten over prøven blev derfor fundet ved at dividere overhøjden<br />

på vandspejlet med prøvens længde, se ligning F1, der er en omskrivning af Darcys<br />

formel.<br />

Figur F1. Billedet viser den opstilling der blev anvendt til bestemmelse af åbundens hydrauliske<br />

ledningsevne. Røret til venstre er vandfyldt, mens det højre rør måler potentialet.<br />

Bilag F


Fejlkilder<br />

Forsøgene er blevet udført en anelse forskelligt, idet pleksiglasrøret i det første forsøg<br />

ikke blev presset nær så langt ned i åbunden som i det andet forsøg og at vandspejlets<br />

overhøjde var større i det første forsøg end i det andet. Begge disse forhold har været<br />

medvirkende til, at trykfaldet over prøven har været større i det første forsøg, end i det<br />

andet, hvilket kan have været medvirkende til, at denne prøve er blevet en smule<br />

kompakteret og derfor ikke har haft den hydrauliske ledningsevne den naturligt vil have.<br />

Udover dette blev de to forsøg lavet på forskellige dage, hvorfor det er muligt, at<br />

undersøgelsen ikke er lavet under helt de samme forhold.<br />

Udover disse forskelle mellem de to forsøg er der nogle mere generelle fejlkilder i<br />

forbindelse med denne type forsøg, det er værd at diskutere. For det første er der<br />

mulighed for, at installationen af røret kan have skabt nogle randbetingelser, der er<br />

forskellig fra resten af søjlen, fordi der let skabes strømningsveje langs kanten.<br />

Derudover er der mulighed for, at sedimenterne i røret kan være blevet kompakteret en<br />

anelse i forbindelse med, at røret er blevet presset ned. Begge disse problemer blev<br />

forsøgt undgået ved at slibe røret, således, at det havde en skarp æg og derfor let gled<br />

gennem sedimenterne. Det var imidlertid umuligt helt at undgå, at de ovennævnte<br />

fejlkilder i nogen grad påvirkede forsøget, da sedimenterne var meget løse og det derfor<br />

var svært, rent praktisk, at installere rørene uden at forstyrre sedimentet. Ingen af de<br />

nævnte fejlkilder skønnes imidlertid at have påvirket forsøget så meget, at resultaterne<br />

ikke kan anvendes, da resultaterne fra denne undersøgelse stemmer godt overens med<br />

resultaterne fra slugtesten og gennemstrømnings-forsøget.<br />

Ligning F1.<br />

Q<br />

K <br />

h<br />

<br />

A<br />

l<br />

<br />

Hvor:<br />

K: Åbundens hydrauliske ledningsevne (cm/s)<br />

Q: Flowet ud af røret (cm 3 /s)<br />

h: Overhøjden af vandspejlet (cm 2 )<br />

l: Den dybde røret var presset ned i åbunden (cm)<br />

A: Tværsnitsarealet af plastic røret (cm 2 )


Bilag G<br />

Diffus udstrømning<br />

For at måle udstrømningen af vand fra bunden blev der anvendt et seepagemeter, et 1 m<br />

langt og 40 cm bredt pleksiglasrør, der var lukket i den ene ende. I den lukkede ende var<br />

der monteret en fitting med en ventil, hvor der blev fastsat en pose til opsamling af<br />

vand. Hele anordningen var vandtæt for at sikre, at kun den diffuse udstrømning blev<br />

målt. Selve undersøgelsen blev udført ved, at pleksiglasrøret forsigtigt blev presset ned i<br />

åbunden og vandfyldt. Herefter blev en pose til opsamling af vand monteret, og<br />

opstillingen blev efterladt i ca. 7 timer. Derefter blev den mængde vand, der var blevet<br />

opsamlet i posen, målt. Under hele forsøget lå posen dækket af vand, således at trykket i<br />

posen var det samme som i vandløbet, og der derfor ikke ville være en trykforskel<br />

indenfor og udenfor seepagemetret. Den vandmængde, der blev opsamlet i posen, svarer<br />

derfor til den vandmængde, der ville strømme ud i vandløbet ved samme tryk.<br />

Opstillingen kan ses på figur G1 . Selve udregningen af udstrømningen foregik ved, at<br />

den opsamlede vandmængde blev målt og ud fra arealet af seepagemetret blev den<br />

specifikke udstrømning udregnet, se ligning G1 . Den samlede udstrømning fra hele<br />

åbunden i vandløbet er derefter estimeret ved, at gange det samlede areal af åbunden (se<br />

under) med den specifikke udstrømning.<br />

Ventil Fitting<br />

Figur G1. Principskitse af seepagemetret.<br />

Pose til<br />

opsamling<br />

af vand<br />

Beregning af åbundens areal<br />

Beregningen blev gjort ud fra målinger af bredde og dybde tre steder langs bækken. Det<br />

samlede bundareal er beregnet til 2.464 m 2 , tabel G1 . Bundarealet er beregnet ved at<br />

udregne bundens bredde i de tre tværsnit. Bundbredden mellem punkterne er<br />

interpoleret og længden er ganget med bredden. Det er antaget, at bundbredden ved<br />

vandløbets udspring er 0.<br />

Bilag G


Ligning G1.<br />

<br />

<br />

Q<br />

A<br />

S <br />

<br />

t<br />

Hvor:<br />

S: Specifik udstrømning (ml/m 2 * s -1 )<br />

Qud: Udsrømningen målt forsøget (ml)<br />

Asep: Arealet af seepagemetret (m 2 )<br />

t: Den tid forsøget varede (s)<br />

Tabel G1. Input til beregning af åbundens areal.<br />

Lok 1 Lok 2 Lok 3<br />

Bredde (m) 2 2 2,1<br />

Dybde (m) 0,25 0,45 0,20<br />

Længde af tværsnittet (m) 2,06 2,19 2,13<br />

Længde af å før punktet (m) 405 104 124<br />

Længde af å efter punktet (m) 104 124 498<br />

Areal i alt (m 2 ) 2.464<br />

Fejlkilder diffus udstrømning<br />

I forbindelse med installationen er de samme fejlkilder behæftet med dette forsøg som<br />

med installeringen af pleksiglasrørerne til målingen af åbundens hydrauliske ledningsevne.<br />

Forsøgets største fejlkilde må imidlertid siges at være, at der kun er udført en<br />

måling af den diffuse udstrømning og at det ikke er sikkert, at resultatet er repræsentativt<br />

for hele vandløbet. Eksempelvis kunne det tænkes, at resultatet ikke er repræsentativt<br />

for det område, der ligger opstrøms for gytjelaget. Det er muligt, at denne øvre del<br />

af vandløbet har en større diffus udstrømning, da der her ikke er et gytjelag til at<br />

begrænse vandets bevægelse. Dette understøttes også af, at der opstrøms for gytjelaget<br />

er en stor tilvækst i vandføringen i forbindelse med forekomsten af flere kilder.<br />

Samtidig kan det siges, at geologien er så homogen i den nedre del af ådalen, hvor<br />

gytjelaget findes, at resultatet sandsynligvis er gældende for denne del af ådalen.<br />

Målingen anses derfor at give en maksimumsudstrømning fra åbunden i det område af<br />

ådalen, hvor gytjelaget findes.<br />

ud<br />

sep

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!