combined_kf_calculator.R

get_realLeakage <- function (area_welllocs = 2500*500, #meter^2
                             area_model = 5400*4000,  #meter^2 
                             kf_welllocs = 1E-7, #meter/day
                             kf_natural = 1E-6 #meter/day
) {
  (area_welllocs * kf_welllocs + (area_model - area_welllocs) * kf_natural)/area_model 
}

boreholes_number <- 500
area_welllocs <- 2500*500 #meter^2
area_model <- 5400*4000  #meter^2 
kf_welllocs <- 1E-7
kf_natural <- c(1E-7, 1E-8, 1E-9, 1E-10, 1E-11)

get_realLeakage(kf_welllocs = kf_welllocs, kf_natural = kf_natural)

plot(log10(kf_natural), 
     log10(get_realLeakage(kf_welllocs = kf_welllocs, 
                           kf_natural = kf_natural)), 
     type = "b", 
     xlab = "Log10 geogener kf-Wert (m/s)", 
     ylab = "Log10 mittlerer kf_Wert (m/s)", 
     xlim=c(-11, -6), ylim=c(-11, -6), 
     main = sprintf("Brunnenfeld-Fl\u00E4che: %d m2 (Anzahl: %d, mittlere Brunnenfeldfl\u00E4che: %1.2f m2)\nGesamtfl\u00E4che: %d m2", 
                    area_welllocs,
                    boreholes_number, 
                    area_welllocs/boreholes_number, 
                    area_model), 
     col = "blue",
     pch = 1,
     las = 1)
text(-10.5, -8.4, 
     labels=bquote(paste("kf"[Brunnenfeld], "= 1E-7 m/s")), 
     col = "blue",
     font=2)
points(log10(kf_natural), 
       log10(get_realLeakage(kf_welllocs = 5E-7, 
                             kf_natural = kf_natural)), 
       type="b", 
       col = "darkgreen",
       pch = 2,
       lty = 2)
text(-10.5, -7.7, 
     labels=bquote(paste("kf"[Brunnenfeld], "= 5E-7 m/s")), 
     col = "darkgreen",
     font=2)
points(log10(kf_natural), 
       log10(get_realLeakage(kf_welllocs = 5E-6, 
                             kf_natural = kf_natural)), 
       type="b", 
       col = "red",
       pch = 3,
       lty = 3)
text(-10.5, -6.7, 
     labels=bquote(paste("kf"[Brunnenfeld], "= 5E-6 m/s")), 
     col = "red",
     font=2)

# text(log10(min(kf_natural))+2, 
#      log10(max(get_realLeakage(kf_boreholes = kf_boreholes, kf_natural = kf_natural))), 
#      labels = "area_boreholes * kf_boreholes + (area_model - area_boreholes) * kf_natural)/area_model")