Ero sivun ”Mikrobien määrä pohjavedessä” versioiden välillä

Opasnet Suomista
Siirry navigaatioon Siirry hakuun
(→‎Laskenta: vieläkin jokin vika)
Rivi 44: Rivi 44:
<rcode graphics="1" variables="
<rcode graphics="1" variables="
category:Anna arvo vain jos se poikkeaa ylläolevasta oletusarvojen taulukosta|
category:Anna arvo vain jos se poikkeaa ylläolevasta oletusarvojen taulukosta|
name:tilanne|type:selection|options:suoravirtaus;Suoravirtaus;kaivoimu;Kaivoimu|default:suoravirtaus|
name:tilanne|type:selection|options:'suoravirtaus';Suoravirtaus;'kaivoimu';Kaivoimu|default:'suoravirtaus'|
name:kaivon.tuotto|default:|
name:tuloste|type:selection|options:'kulkeutumisaika';Kulkeutumisaika;'mikrobipit.kohde';Mikrobipitoisuus kohteessa;'inaktivaatioaika';Inaktivaatioaika|default:'kulkeutumisaika'|
name:kaivon.tuotto|default:1000|
name:vedenjohtavuus|default:|
name:vedenjohtavuus|default:|
name:paksuus|default:|
name:paksuus|default:|
Rivi 57: Rivi 58:
library(xtable)
library(xtable)


# Ensin valikoidaan tietolähde. Käyttäjän antama input on ensisijainen, jos sitä ei ole, kaivetaan Basesta oletusarvo data.
apparate <- function(objname, source = "Op_fi1755") {
obj <- get(objname)
if(is.null(obj)) obj <- NA
if(is.na(obj) | obj == "") {
temp <- opbase.data(source)
obj <- temp[temp$Parametri == objname & temp$Havainto == "Arvo", ]["Result"]
}
if(is.vector(obj)) obj <- data.frame(Result = obj)
if(is.data.frame(obj)) obj <- new("ovariable", name = objname, data = obj)
return(obj)
}


data <- opbase.data("Op_fi1755")
kaivon.tuotto <- apparate("kaivon.tuotto")
data <- data[data$Havainto == "Arvo", ]
vedenjohtavuus <- apparate("vedenjohtavuus")
paksuus <- apparate("paksuus")
etäisyys <- apparate("etäisyys")
huokoisuus <- apparate("huokoisuus")
mikrobipit.lähde <- apparate("mikrobipit.lähde")
inaktivaatiopit <- apparate("inaktivaatiopit")
inaktivaationopeus <- apparate("inaktivaationopeus")


if(is.null(kaivon.tuotto) | kaivon.tuotto == '')
# Ensin valikoidaan tietolähde. Käyttäjän antama input on ensisijainen, jos sitä ei ole, kaivetaan Basesta oletusarvo data.
kaivon.tuotto <- data[data$Parametri == "kaivon.tuotto", ]["Result"]
if(is.null(vedenjohtavuus) | vedenjohtavuus == '')
vedenjohtavuus <- data[data$Parametri == "vedenjohtavuus", ]["Result"]
if(is.null(paksuus) | paksuus == '')
paksuus <- data[data$Parametri == "paksuus", ]["Result"]
if(is.null(etäisyys) | etäisyys == '')
etäisyys <- data[data$Parametri == "etäisyys", ]["Result"]
if(is.null(huokoisuus) | huokoisuus == '')
huokoisuus <- data[data$Parametri == "huokoisuus", ]["Result"]
if(is.null(mikrobipit.lähde) | mikrobipit.lähde == '')
mikrobipit.lähde <- data[data$Parametri == "mikrobipit.lähde", ]["Result"]
if(is.null(inaktivaatiopit) | inaktivaatiopit == '')
inaktivaatiopit <- data[data$Parametri == "inaktivaatiopit", ]["Result"]
if(is.null(inaktivaationopeus) | inaktivaationopeus == '')
inaktivaationopeus <- data[data$Parametri == "inaktivaationopeus", ]["Result"]
 
if(is.vector(kaivon.tuotto)) kaivon.tuotto <- data.frame(Result = kaivon.tuotto)
if(is.vector(vedenjohtavuus)) vedenjohtavuus <- data.frame(Result = vedenjohtavuus)
if(is.vector(paksuus)) paksuus <- data.frame(Result = paksuus)
if(is.vector(etäisyys)) etäisyys <- data.frame(Result = etäisyys)
if(is.vector(huokoisuus)) huokoisuus <- data.frame(Result = huokoisuus)
if(is.vector(mikrobipit.lähde)) mikrobipit.lähde <- data.frame(Result = mikrobipit.lähde)
if(is.vector(inaktivaatiopit)) inaktivaatiopit <- data.frame(Result = inaktivaatiopit)
if(is.vector(inaktivaationopeus)) inaktivaationopeus <- data.frame(Result = inaktivaationopeus)
 
if(is.data.frame(kaivon.tuotto) kaivon.tuotto <- new("ovariable", name = "kaivon.tuotto", data = kaivon.tuotto)
if(is.data.frame(vedenjohtavuus) vedenjohtavuus <- new("ovariable", name = "vedenjohtavuus", data = vedenjohtavuus)
if(is.data.frame(paksuus) paksuus <- new("ovariable", name = "paksuus", data = paksuus)
if(is.data.frame(etäisyys) etäisyys <- new("ovariable", name = "etäisyys", data = etäisyys)
if(is.data.frame(huokoisuus) huokoisuus <- new("ovariable", name = "huokoisuus", data = huokoisuus)
if(is.data.frame(mikrobipit.lähde) mikrobipit.lähde <- new("ovariable", name = "mikrobipit.lähde", data = mikrobipit.lähde)
if(is.data.frame(inaktivaatiopit) inaktivaatiopit <- new("ovariable", name = "inaktivaatiopit", data = inaktivaatiopit)
if(is.data.frame(inaktivaationopeus) inaktivaationopeus <- new("ovariable", name = "inaktivaationopeus", data = inaktivaationopeus)
 
print(kaivon.tuotto )
print(vedenjohtavuus)
print(paksuus)
print(etäisyys)
print(huokoisuus )
print(mikrobipit.lähde)
print(inaktivaatiopit)
print(inaktivaationopeus)


# Funktio laskee mikrobipitoisuuden pohjavedessä.
# Funktio laskee mikrobipitoisuuden pohjavedessä.


dependencies <- data.frame()
dependencies <- data.frame(Name = c(
"tilanne",
"tuloste",
"kaivon.tuotto",
"vedenjohtavuus",
"paksuus",
"etäisyys",
"huokoisuus",
"mikrobipit.lähde",
"inaktivaatiopit",
"inaktivaationopeus"
))


formula <- function(dependencies, ...) {
formula <- function(dependencies, ...) {
 
ComputeDependencies(dependencies, ...)
cat("Tässä ollaan\n")
kulkeutumisaika <- 0
kulkeutumisaika <- 0
if(tilanne == "suoravirtaus") {
if(tilanne == "suoravirtaus") {
kulkeutumisaika <- etäisyys / (vedenjohtavuus * (paksuus / etäisyys) / huokoisuus) / 3600 / 24
kulkeutumisaika <- etäisyys / (vedenjohtavuus * (paksuus / etäisyys) / huokoisuus) / 3600 / 24
}
}
cat("Tässä ollaan\n")
if(tilanne == "kaivoimu") {
if(tilanne == "kaivoimu") {
kulkeutumisaika <- 2 * pi * huokoisuus * paksuus * etäisyys ^ 2 / kaivon.tuotto # Miksi 2 * ?
kulkeutumisaika <- 2 * pi * huokoisuus * paksuus * etäisyys ^ 2 / kaivon.tuotto # Miksi 2 * ?
}
}
cat("Tässä ollaan\n")
print(mikrobipit.lähde)
print(inaktivaationopeus)
print(kulkeutumisaika)


mikrobipit.kohde <- mikrobipit.lähde * exp(-inaktivaationopeus * kulkeutumisaika)
mikrobipit.kohde <- mikrobipit.lähde * exp(0 - inaktivaationopeus * kulkeutumisaika)
inaktivaatioaika <- log10(inaktivaatiopit / mikrobipit.lähde) / (-inaktivaationopeus / 2.3) # Miksi 2.3?
cat("Tässä ollaan\n")
out <- list(kulkeutumisaika, mikrobipit.kohde, inaktivaatioaika)
inaktivaatioaika <- log10(inaktivaatiopit / mikrobipit.lähde) / (0 - inaktivaationopeus / 2.3) # Miksi 2.3?
names(out) <- c("kulkeutumisaika", "mikrobipit.kohde", "inaktivaatioaika")
cat("Tässä ollaan\n")


return(out)
return(get(tuloste))
}
}


Rivi 134: Rivi 126:
)
)


mikrobit.pohjavesi(
#mikrobit.pohjavesi
tilanne = "suoravirtaus",
mikrobit.pohjavesi <- EvalOutput(mikrobit.pohjavesi, N = 10)
kaivon.tuotto = 1000,
mikrobit.pohjavesi
vedenjohtavuus = 1E-5,
# paksuus = 10,
etäisyys = 100,
huokoisuus = 0.35,  
mikrobipit.lähde = 1E8,
inaktivaatiopit = 10,
inaktivaationopeus = 0.01
)
 
#Koodi laskee mutta väärin. Annetaanko inputiksi vahingossa faktoreita?
 
tilanne = "suoravirtaus",
kaivon.tuotto = NA,
vedenjohtavuus = NA,
paksuus = NA,
etäisyys = NA,
huokoisuus = NA,
mikrobipit.lähde = NA,
inaktivaatiopit = NA,
inaktivaationopeus = NA
) {


</rcode>
</rcode>

Versio 2. joulukuuta 2012 kello 21.31


Kysymys

Mitä mikrobien määrä pohjavedessa kuvaa?

Vastaus

Mikrobien määrä pohjavedessä kuvaa eri mikrobien mitattuja pitoisuuksia pohjavedessä.

Perustelut

Data

Data perustuu kuviteltuun tilanteeseen.

Ero sivun ”Mikrobien määrä pohjavedessä” versioiden välillä(-)
ObsSymboliParametriSeliteArvoYksikköExcel
1KsvedenjohtavuusKyllästyneen tilan vedenjohtavuus1.00E-05m /sD9
2HspaksuusPohjavesivyöhykkeen paksuus5-15mD10
3RetäisyysEtäisyys päästölähteestä100mD11
4nhuokoisuusMaaperän huokoisuus0.35D12
5C0mikrobipit.lähdePäästölähteen mikrobipitoisuus = mikrobipitoisuus akviferissä1.00E+08# /lD13
6CAinaktivaatiopitPitoisuus jolloin mikrobien ajatellaan inaktivoituneen10# /lD14
7µinaktivaationopeusInaktivaatiokerroin0.011 /vrkD15
8tkulkeutumisaikaKulkeutumisaika405.0925926vrkD17: =D11/(D9*(D10/D11)/D12)/3600/24
9Ctmikrobipit.kohdePitoisuus tarkastelukohdassa1.74E+06# /lD18: =D13*EXP((-D15*D17))
10tininaktivaatioaikaInaktivaatioaika1610vrkD19: =LOG10(D14/D13)/(-D15/2.3)

Riippuvuudet

Yksikkö

# /l

Laskenta

Tämä laskenta kuvaa mikrobien kulkeutumista pohjavedessä tai maaperässä.

Yksinkertaisimmillaan mikrobien voidaan arvioida kulkeutuvan päästölähteestä raakavesilähteeseen samalla nopeudella kuin vesi virtaa kyllästyneessä maaperässä. Päästölähteen oletetaan sijaitsevan pohjavesivyöhykkeessä ja virtauksen tapahtuvat lateraalisuunnassa. Päästölähteen virtaamaa ei ole tunnettu, joten pohjavesivyöhykkeen ja päästölähteen ja pohjavesivyöhykkeen mikrobipitoisuudet oletetaan samoiksi. Pumppausta ei ole huomioitu. Tiedot ovat tiedostosta Riskinarviointilaskelmat_inaktivaatio271112.xls.

Anna arvo vain jos se poikkeaa ylläolevasta oletusarvojen taulukosta

tilanne:

tuloste:

kaivon.tuotto:

vedenjohtavuus:

paksuus:

etäisyys:

huokoisuus:

mikrobipit.lähde:

inaktivaatiopit:

inaktivaationopeus:

+ Näytä koodi

Tulokset

Clean Medium Contaminated
Norovirus 0 204.167 2166.667
Rotavirus 0 204.167 2166.667
Murine norovirus 0 204.167 2166.667

Katso myös

Viitteet


Aiheeseen liittyviä tiedostoja

<mfanonymousfilelist></mfanonymousfilelist>

Noudettu kohteesta ”https://fi.opasnet.org/index.php?title=Mikrobien_määrä_pohjavedessä&oldid=16169”