Afprøvning af reaktiv jernvæg til grundvandsrensning

Bilag 6
Opbygning og kalibrering af grundvands- og stoftransportmodellerne

Grundvandsmodellen

I forhold til de tidligere modeller er modelområdet udvidet mod nordvest for at minimere grænsebetingelsernes indflydelse på modelberegningerne. Det reviderede modelområde er gengivet på figur 1.

Figur 1
Modelområde

Ved modelberegningerne er grundvandspotentialet fastholdt langs den nordvestlige rand, idet afværgeforanstaltningerne ved Vapokon ikke vurderes at have indflydelse på grundvandspotentialet i randområdet. Modellens sydlige rand udgøres af Holmebækken. Da der vurderes at være hydraulisk kontakt mellem Holmebækken og det primære grundvandsmagasin i området er grundvandspotentialet langs modellens sydlige rand ligeledes fastholdt. Endelig er modellens øst-, vest- og nordlige rande placeret parallelt med grundvandets strømningsretning således, at der ikke sker vandudveksling på tværs af randen.

For at øge modellens følsomhed er cellestørrelsen nedsat til 2 m i begge horisontale retninger, og antallet af lag er øget til 5, jf. figur 2. De 5 modellag repræsenterer følgende lagfølge:
Modellag 1 repræsenterer det terrænnære fyld og moræneler.
Modellag 2 repræsenterer det underliggende moræneler med sandlinser i området omkring den reaktive væg. I den øvrige del af området repræsenterer modellag 2 toppen af det sandlag, der udgør det primære magasin i området.
Modellag 3 - 5 repræsenterer de sandlag, der udgør det primære magasin. Bunden af modellag 5 er således fastlagt ud fra overgangen til det nedre moræneler i området. De vertikale afgrænsninger mellem modellag 3 - 5 er så vidt muligt fastlagt ud fra ændringer i kornstørrelse, lerindhold og hydraulisk ledningsevne.

Se her!

Figur 2
Vertikalt snit gennem modelområdet.

De anvendte modelparametre for hvert enkelt lag fremgår af tabel 1. Da modellen udelukkende er anvendt til stationære beregninger, er magasintallet ikke defineret.

Tabel 1:
Anvendte modelparametre

Modellag

Top kote

Bund kote

Geologisk bedømmelse

Porøsitet

Magasin type

Hydraulisk ledningsevne m/s

1

21,5- 25,1

17,2- 22,8

Fyld og moræneler

0,4

Frit

1x10

2
Området omkring den reaktive væg

 

Den øvrige del af området

 

 

17,2- 22,8

 

 

16,4- 20,6

Moræneler med sandstriber
 
  
Sand usorteret

0,4
 
 
 
  
0,35

 

 

Frit/ artesisk

1x10-7
1x10-5
 
 
 
1,5x10-4

3

16,4- 20,6

15,3- 18,6

Sand fint-mellem

0,35

Frit/ artesisk

1,5x10-4

4

15,3- 18,6

13,4- 15,7

Sand mellem-groft.

0,35

Artesisk

3x10-4

5

13,4- 15,7

9,5- 14,3

Sand mellem

0,35

Artesisk

2x10-4


Den horisontale hydrauliske ledningsevne for det primære magasin er reduceret væsentligt i forhold til den hydrauliske ledningsevne anvendt i de tidligere modeller på 8×10-4 m/s. Reduktionen er foretaget på baggrund af de udførte slugtests og sigteprøver samt vandbalance vurderinger og revurdering af de tidligere udførte prøvepumpninger. Baggrund for reduktionen i de horisontale hydrauliske ledningsevner er beskrevet nærmere i appendix A.

Der findes ingen målinger af den vertikale hydrauliske ledningsevne. I modellen er anvendt værdier på ca. 1/10 af den horisontale hydrauliske ledningsevne, hvilket erfaringsmæssigt er gældende for materiale aflejret af smeltevand.

Nedsivningen er på baggrund af en overordnet vandbalance for drikkevandsområdet ved Søndersø gengivet i /8/ og /9/ skønnet til 40 mm/år i de ubefæstede dele af modelområdet. I de befæstede dele af området er nedsivningen skønnet til ca. 3 mm/år svarende til udsivningen fra utætte kloakker o. lign. Fordelingen mellem befæstede og ubefæstede arealer fremgår af figur 3.

Figur 3
Befæstede og ubefæstede arealer

Stoftransportmodellen

Den opstillede grundvandsmodel er udvidet med en stoftransportmodel, der kan simulere transport og spredning af stof opløst i grundvandet, men ikke sorption og nedbrydning. Sorption og nedbrydning vurderes at have størst betydning ved transport over længere afstande eller tid. Udeladelsen af disse vurderes derfor ikke at have væsentlig betydning for de beregnede stofkoncentrationen i drænvandet, idet drænet er placeret relativt tæt på kildeområdet, og de højeste koncentrationer i drænvandet forventes umiddelbart efter etableringen af drænet.

Stofspredningen i modellen sker ved dispersion med nedenstående dispersiviteter (tabel 2)

Tabel 2:
Dispersiviteter anvendt i modellen

Modellag

Langsgående dispersivitet

Tværgående horisontal dispersivitet

Tværgående vertikal dispersivitet

1

1,0

0,001

0,05

2

0,5

0,001

0,005

3

0,5

0,001

0,005

4

0,5

0,001

0,005

5

0,5

0,001

0,0005


Kalibrering af modellerne

Den opstillede grundvandsmodel er kalibreret med henblik på at opnå en rimelig overensstemmelse mellem de målte grundvandspotentialer i det primære magasin og de modellerede potentialer. Disse fremgår af figur 4.

Figur 4.
Kalibreret grundvandspotentiale

Modellen er ligeledes kalibreret med henblik på at opnå en rimelig overensstemmelse mellem beregnet/målte og modellerede værdier af:
indstrømning til modelområdet
partikelhastigheden ved den reaktive væg

Disse fremgår af tabel 3. Heraf ses at der er god overensstemmelse mellem beregnet/målt og modelleret partikelhastighed ved den reaktive væg, mens indstrømningen til modelområdet er ca. dobbelt så stor i modellen, som den beregnede værdi. Dette medfører, at der i modellen sker en mindre overestimering af vandmængder og partikelhastighed. Denne overestimering vurderes dog at ligge indenfor usikkerheden på de anvendte modelparametre.

Tabel 3:
Overensstemmelse mellem beregnede/målte og modellerede værdier

 

Beregnet/målt værdi

Modelleret værdi

Indstrømning til modelområdet

27 - 76 m3/d

141 m3/d

Partikelhastighed ved reaktiv væg

25 - 123 m/år

69 m/år


Den opstillede stoftransportmodel er kalibreret med henblik på at opnå en rimelig overensstemmelse mellem målte og modellerede koncentration, udbredelse og mængde af chlorerede komponenter. Kalibreringen er foretaget ved en justering af dispersiviteterne og kildestyrken. Kildestyrken er ved kalibreringen nedsat til 60 kg/år, mens de kalibrerede dispersiviteter fremgår af tabel 2.

Af figur 5 og 6 fremgår udbredelsen af chlorerede komponenter i henholdsvis toppen af magasinet (modellag 3) og bunden af magasinet (modellag 5) efter kalibrering med en simuleringsperiode på 10 år. I tabel 4 er de tilsvarende simulerede koncentrationer af chlorerede opløsningsmidler i en række boringer sammenlignet med de målte koncentrationer.

Ved kalibreringen er der lagt størst vægt på at opnå overensstemmelse mellem simulerede og målte koncentrationer i toppen af magasinet omkring Vapokon-grunden, hvor drænet er placeret.

Stoftransportmodellen kan som tidligere nævnt ikke simulere tilførelse af forureningskomponenter ved opløsning af fri fase fra pools beliggende på bunden af magasinet. Det medfører, at modellen underestimerer forureningskoncentrationerne i bunden af magasinet umiddelbart under kildeområdet. Modellen tager endvidere ikke højde for den naturlige nedbrydning af de chlorerede komponenter. Det betyder, at modellen overestimerer koncentrationerne i området syd for den reaktive væg, hvor den naturlige nedbrydning er mest udtalt. Disse uoverensstemmelser vurderes dog ikke at have indflydelse på beregningen af stofkoncentrationen i drænvandet, idet der ved beregningen anvendes de faktiske koncentrationer i magasinet.

Figur 5.
Kalibreret udbredelse af chlorerede komponenter i toppen af magasinet (moddellag 3)
   

Figur 6.
kalibreret udbredelse af chlorerede komponenter i bunden af magasinet (moddellag 5)

 
Tabel 4:

Målte og simulerede koncentrationer

Boring

Koncentration af chlorerede komponenter (µg/l)

Målt

Simuleret i model lag

2

3

4

5

B 6

43

179

207

 

 

E 5

582

 

240

71

 

E 7.1

19

 

 

 

154

E 8.2

4195

3575

3910

 

 

E 8.1

348

 

 

1026

247

E 9

56

 

297

257

 

E 10

13

1410

1317

 

 

E 12

1.3

 

 

36

9

G 21.2

32660

 

4527

 

 

G 21.1

10

 

 

 

0.4

G 22.2

0.8

 

0.04

0.00

 

G 22.1

1.5

 

 

 

0.00

G 23.2

33921

 

23409

99

 

G 23.1

1677

 

 

 

0.8

G 29.3

10041

 

5290

 

 

G 29.2

4173

 

 

458

 

E 4.1

1654

 

 

458

 

G 29.1

27

 

 

 

 

G 33.2

16075

3019

 

 

 

G 33.1

61

 

 

 

249

G 34.2

47352

4871

5357

 

 

G 34.1

173

 

 

1170

255

 

Appendix A

Revurdering af den hydrauliske ledningsevne i det primære magasin ved Vapokon.

Ved skitseprojekteringen af reaktiv væg i Søndersø er der anvendt en k-værdi på 0,0008 m/s og en naturlig strømningshastighed på 387 m/år. Disse værdier virker meget høje, og der skal derfor her gives en vurdering af rimeligheden.

Vandbalance vurdering

I det følgende er der foretaget et skøn af den hydrauliske ledningsevne på baggrund af et skøn over grundvandstilstrømningen til Vapokon og de tidligere konstaterede grundvandsgradienter og magasintykkelser.

Grundvandstilstrømningen (q) til Vapokon er lig med den grundvandsdannelse (Ig), som finder sted mellem Vapokon og grundvandsskellet. Dvs.

q = Ig * L

hvor L er afstanden fra grundvandsskellet langs en strømningslinie til Vapokon.

På baggrund af et potentialekort for grundvandsstanden /1/ vurderes L til ca. 5000 m (se figur 1) - usikkerheden på denne vurdering skønnes til +/- 1000 m.

Ifølge Kortlægning og klassificering af grundvandsressourcen i Fyns Amt, /1/, er den totale grundvandsdannelse Ig i Søndersø 33 mm/år. Usikkerheden på dette skøn vurderes til +/- 10 mm/år.

Grundvandstilstrømningen til Vapokon kan således beregnes til 165 m3/år/m på tværs af strømningsretningen, med et usikkerhedsinterval på [92; 258 m3/år/m].

Magasinets tykkelse (m) er, jf. /2/, mellem 6 og 10 meter på Vapokon med et gennemsnit på ca. 7,5 m. Heraf følger, at filterhastigheden ( v = q/m) bliver lig 22 m/år med et usikkerhedsinterval på [9; 43] m/år.

Partikelhastigheden af grundvandet er defineret som filterhastigheden divideret med den effektive porøsitet, som er sat til 0,3. Partikelhastigheden bliver således 73 m/år, med et usikkerhedsinterval på [30; 143] m/år.

Ifølge /2/ (figur 5.3 og 5.5) er grundvandets gradient i det primære magasin 5,4 promille, mens gradienten i /3/ (figur 6.10) er 4,7 promille.

Med en middelgradient på 5 promille og en grundvandstilstrømning på 165 m3/år/m fås en transmissivitet T på 10-3 m2/s med et usikkerhedsinterval (q=92 til 258) på [1,7×10-3; 0,5×10-3] m2/s.

Ved en magasintykkelse (m) på mellem 6 og 10 meter på Vapokon, med et gennemsnit på ca. 7,5 m, bliver den hydrauliske ledningsevne (k = T/m) 1,3×10-4 m/s, med et usikkerhedsinterval på [2,8×10-4; 0,5×10-4] m/s.

Vurdering af sigteanalyser

RAMBØLL har foretaget 3-4 sigteanalyser af sandmagasinet ved den reaktive væg. Den hydrauliske ledningsevne er vurderet ved hjælp af Takis formel: k = 0,01 * d210 , hvor d10 er kornstørrelsen i mm fra sigteanalysen. Herved findes en hydraulisk ledningsevne på mellem 1,0×10-4 og 1,7×10-4 m/s med en usikkerhed på ca. 50%.

Vurdering af prøvepumpning i 1996

Der er tidligere foretaget en prøvepumpning /2/ med 5 m3/h fra boring PB1 fra den 21/05/96, kl. 18:54, til den 26/05/96, kl. 14:34. Efterfulgt af en tilbagepejling indtil den 28/05/96, kl. 11:00. Ved PB1 er sænkningen ikke større, end at vandstanden står over top af sandmagasinet under hele prøvepumpningen. Det ses, at der er stor forskel mellem sænkning og stigning - så stor at den ikke kan forklares med fortsat sænkning på grund af prøvepumpningen. Generelt er påvirkningen af grundvandsmagasinet ikke særlig stor, hvilket øger usikkerheden på tolkningen.

Der er ikke taget hensyn til sæsonvariationer ved analyse af prøvepumpningen. Dette kan give anledning til øget usikkerhed, idet der tilsyneladende er en generel sænkning af vandspejlet i perioden. I /3/ er der lagt vægt på resultatet fra stigningen, selvom denne kun har varet i ca. 45 timer, mens sænkningen har varet 116 timer. Resultat for pumpeboringen er givet ved en transmissivitet T = 0,0067 m2/s for stigningsperioden og T = 0,00373 m2/s for sænkningsperioden. En retolkning af sænkningen tyder på T= 0,00245 m2/s. Da datagrundlaget er mere sikkert for sænkningerne, er der her taget udgangspunkt i disse, hvilket betyder, at de beregnede transmissiviteter bliver næsten halveret. Det skal bemærkes, at boringerne B7 og B8, i modsætning til de andre boringer i prøvepumpningen, står i den frie del af grundvandsmagasinet, hvilket her kan vanskeliggøre tolkningen.

Vurdering af prøvepumpning på rensningsanlægget

I forbindelse med dimensionering af en grundvandssænkning på rensningsanlægget Snavevej 27 blev der, jf. /3/, foretaget en prøvepumpning i efteråret 1997 med en varighed på 50 minutter. Ved prøvepumpningen blev der overvejende fundet transmissiviteter på mellem 0,002 til 0,003 m2/s. Kun ved boring E7 er der fundet en transmissivitet på 0,008 m2/s, hvilket forklares med, at grundvandsmagasinet er 10 meter tykt, mens det ved de øvrige boringer er ca. 7 meter. Der er dog andre muligheder, idet grundvandsmagasinet er frit eller overgår til at blive frit ved E7, mens det er artetisk ved de øvrige boringer. Dette kan forklare, hvorfor sænkningstakten er relativ lav ved E7.

Vurdering af slugtest

I forbindelse med forundersøgelsen er der udført slugtest på boringerne G29, G33 og G34. Tolkningen af disse slugtests fremgår af nedenstående tabel.

Boring

Filter-interval
(m u.t.)

Hydraulisk ledningsevne
(m/s)

G29

4-5
5-6
6-7
7-8
8-9
9-10
10-11
11-12

1.69·10-4
6.98·10-5
6.66·10-5
3.74·10-6
1.68·10-4
9.48·10-5
3.54·10-7
1.02·10-7

G33

2.5-3.5
3.5-4.5
4.5-5.5
5.5-6.5
6.5-7.5
7.5-8.5

2.43·10-5
1.33·10-4
3.12·10-5
6.61·10-5
3.86·10-5
3.01·10-5

G34

7.2-9.2
2.2-4.2

2.09·10-4
2.50·10-4


Det skal bemærkes, at boring G29 er placeret lige ved siden af boring E4, og det derfor må forventes, at de har samme hydrauliske ledningsevne.

Konklusion

I den gamle model er der anvendt en hydraulisk ledningsevne på 8×10-4 svarende til en partikelhastighed på 387 m/år ved den reaktive væg, mens vandbalance vurderingen tyder på en hastighed på mellem 30 og 143 m/år med et gennemsnit på 73 m/år.

I nedenstående tabel er der lavet en samstilling af de hydrauliske ledningsevner, som er fremkommet ved anvendelse af de forskellige metoder.

Det ses, at alle de nye metoder giver hydraulisk ledningsevner på mellem 0,00005 og 0,0005 m/s. Undtagelserne er boring E7, B7, og B8, som ikke vurderes at være repræsentative for en naturlig grundvandsstrømning uden oppumpning.

Det konkluderes derfor, at en gennemsnitlig hydraulisk ledningsevne på 0,00025 m/s må anses for at være rimelig for det primære magasin i området omkring Vapokon.

Metode

Transmissivitet
  
m2/s

Magasin tykkelse
m

Hydraulisk ledningsevne
m/s

Tidligere modeller

0,008

 

0,0008

Vandbalance vurdering

0,001

7,5

0,00005-
0,00028

Sigtekurve vurdering

 

 

0,0001 -
0,00017

HOH-prøvepumpning
Boring AF1
Boring E8
Boring E7
Boring AF2
Boring E9
Boring E10

0,00242
0,00293
0,00828
0,00204
0,00202
0,00202

7
7,3
8,4
7
7
7,5

0,0003
0,0004
0,0010
0,0003
0,0003
0,0003

B&K-prøvepumpning
Boring PB1
Boring B5
Boring B7
Boring B8
Boring E4-2

 
0,00245 / 0,00672
0,00310/ 0,00575
/ 0,00573
/ 0,00613
0,00301 / 0,00669

 
7,4
9,5
7,5
8,8
5,7

 
0,0003
0,0003
0,0005

Slugtests
Boring G29
Boring G33
Boring G34

 

 

0,0001
0,00005
0,0001

 

Referencer

/1/ Forslag til udpegning af områder med særlige drikkevandsinteresser. Fyns Amt. Maj 1997.
[Tilbage]
 
/2/ Vapokon Petrokemisk Værk A/S. Supplerende forureningsundersøgelser indenfor Vapokon grunden. Birch & Krogboe for Fyns Amt. 1996-07-09.
[Tilbage]
 
/3/ Skitseprojekt for reaktiv væg. Vapokon, Søndersø. Fyns Amt. Juli 1998.
[Tilbage]