Estimation of the variation in the integrated hydraulic conductivity around the Tunåsen infiltration facility : An investigation of the effects of uncertainties in existing information

Estimation of the variation in the integrated hydraulic conductivity around the Tunåsen infiltration facility - An investigation of the effects of uncertainties in existing information Almost half of Sweden’s drinking water volume is produced from groundwater. The main fraction of this is extracted from eskers, some of which allow for very large extraction rates. Despite this, the groundwater volume is not sufficient in some areas. This has led to an extensive use of artificial recharge. In Uppsala, the total volume added through artificial recharge is divided between four infiltration facilities. The largest of these facilities is Tunåsen basin infiltration facility, which is situated along the Uppsala esker. For Tunåsen as well as for other recharge facilities, it is of great importance that the infiltrated water is allowed a sufficient transport time before it is extracted further down the flow path. This is in order to ensure that the water obtains characteristics similar to those of natural groundwater. A way to estimate the transport time is to model the groundwater flow in the area between the infiltration facility and the extraction site. However, due to the fact that groundwater flow in eskers often is very complex, modeling such a system requires data of high quality and resolution. The aim of this master’s thesis was to document the variation in estimated values of integrated horizontal hydraulic conductivity and to investigate how the variation was affected by uncertainties connected to existing information for the area between the Tunåsen infiltration facility and the extraction sites in Storvad and Galgbacken. The thesis was done as a pilot study for Uppsala vatten och avfall AB’s upcoming construction of a high-resolution model of the Uppsala and Vattholma eskers. The investigation began with a literature review, from which existing information was compiled. An integrated horizontal hydraulic conductivity was calculated based on two types of information: descriptions of the layers in borehole profiles and information from sieve analyses. The obtained values were then compared to a rough estimate based on information about groundwater flow and hydraulic gradient. The compilation of existing material showed that the largest uncertainties were related to borehole profiles. This was due to the use of different ways to describe soil layers and to lack of information of a soil layer’s true grain size composition. The calculations of integrated hydraulic conductivity showed the largest variation when based on descriptions of borehole profiles. The calculations that were made based on sieve analyses as well as the rough estimate based on groundwater flow and hydraulic gradient showed a smaller variation. However, without the existence of measurements of the actual hydraulic conductivity, it is impossible to say for certain how the documented uncertainties affect the variation. In the upcoming work with the construction of the high-resolution model it is therefore of great importance that such information is obtained, for example by performing pumping tests. Skattning av den integrerade hydrauliska konduktivitetens variation kring Tunåsens infiltrationsanläggning - En utredning av påverkan från möjliga osäkerheter i befintlig information I Sverige utgörs nästan hälften av den totala dricksvattenvolymen av grundvatten. Den största delen av grundvattnet utvinns ur rullstensåsar, vilka i vissa fall tillåter mycket stora uttag. Trots detta är grundvattenvolymen på vissa håll inte tillräcklig. Detta har lett till att den naturliga grundvattenmängden i sex procent av landets grundvattentäkter kompletteras genom konstgjord infiltration. I Uppsala sker konstgjord infiltration på fyra platser, där Tunåsens infiltrationsanläggning belägen på Uppsalaåsen är den största. För Tunåsen såväl som för andra områden där konstgjord infiltration används, är det av stor betydelse att det infiltrerade vattnet har en tillräckligt lång transporttid mellan infiltrations- och uttagspunkt för att möjliggöra att det får grundvattenliknande egenskaper innan det pumpas upp. Ett sätt att uppskatta denna transporttid är genom att modellera grundvattenflödet i området mellan infiltrationsanläggningen och den eller de grundvattentäkter där vattnet pumpas upp. Till följd av att grundvattenströmningen i rullstensåsar ofta är komplex, ställer dock en modellering av ett sådant område stora krav på både kvalitet och upplösning hos den information från vilken modellen byggs upp. Syftet med detta examensarbete var att dokumentera variationen i en skattad integrerad horisontell hydraulisk konduktivitet och utvärdera hur denna påverkades av osäkerheter relaterade till den idag befintliga informationen för området mellan Tunåsens infiltrationsanläggning och Storvads respektive Galgbackens grundvattentäkt. Detta gjordes som en förstudie till Uppsala vatten och avfall ABs kommande arbete med upprättandet av en noggrann modell över hela Uppsala- och Vattholmaåsen. Arbetet började med en litteraturstudie och en informationsinsamling, genom vilka den idag tillgängliga informationen sammanställdes. Utifrån den sammanställda informationen beräknades sedan värden på hydraulisk konduktivitet baserat på beskrivningarna i de upprättade lagerföljderna och baserat på information från tidigare utförda siktanalyser. De beräknade värdena jämfördes därefter med värden framtagna genom en överslagsberäkning baserad på information om grundvattenflöde och hydraulisk gradient. Sammanställningen av befintlig information visade på att de mest betydande osäkerheterna var relaterade till jordlagerdata, där de kunde kopplas till såväl användandet av olika beskrivningssätt i lagerföljderna som till brister i informationen om verklig kornstorleksfördelning. Skattningarna av integrerad horisontell hydraulisk konduktivitet resulterade i störst variation när de baserades på jordlagerföljderna, då det högsta värdet i många fall var 107 gånger så stort som det lägsta. De beräkningar som utfördes på siktanalyser visade på en mer begränsad variation, vilket även var fallet för överslagsberäkningarna. Utan mätningar på den verkliga hydrauliska konduktiviteten går det dock inte att säkert avgöra hur mycket de funna osäkerheterna påverkar variationen. Inför utvecklingen av den kommande modellen är det därför nödvändigt att komplettera den idag befintliga informationen med ny data, till exempel genom att provpumpningar utförs