man/examples/example.R

In teldosas/wellsolve: Minimize the Total Pumping Cost from an Aquifer to a Water Tank, via a Pipe Network

library(wellsolve)

projected <- "+proj=utm +zone=34"
x0 = 650000
y0 = 4400000

# the first value is the tank and the rest are wells
xc <- c(0,100,180,100,700,800,900)
xc <- xc + x0 - 10000

yc <- c(0,0,0,80,0,0,900)
yc <- yc + y0 - 10000

# xstart is a wild guess about the result.
# This helps the algorithm search for the result in the correct area
xstart <- c(0.250,0.1,0.25,0.1,0.1,0.01)

nw=7 #number of wells + 1 (the tank)

#                   6
#                  /
#                 /
#       3        /
#      /        /
#    /        /
#  /        /
# T---1----2------4----5
#
# T=Tank 1-6=Wells no. 1-6
links <- data.frame(matrix(ncol=2, nrow=(nw-1)))
links[1:(nw-1),1] <- 1:(nw-1)
links[1,2] <- 0
links[2,2] <- 1
links[3,2] <- 0
links[4,2] <- 2
links[5,2] <- 4
links[6,2] <- 2

Dlist <- list()
Dlist[1:6] <- 0.05


# number of extra wells
nwextra <- 2

xextra <- c(500,800)
xextra <- xextra + x0 - 10000

yextra <- c(500,300)
yextra <- yextra + y0 - 10000

qextra <- c()
qextra[1] <- 0
qextra[2] <- 0

result <- wellsolve(projected, T=0.0025 ,R=3000, ff=0.03, qtot<-0.5, r0=0.200,
                  Dlist, qextra, nw, xc, yc, links,
                  nwextra, xextra, yextra, xstart)

# Case A
Dlist[1:6] <- 2
result <- wellsolve(projected, T=0.0025 ,R=3000, ff=0.03, qtot<-0.5, r0=0.200,
                    Dlist, qextra, nw, xc, yc, links,
                    nwextra, xextra, yextra, xstart)

# Case LC1
result <- wellsolve(projected, T=0.0025 ,R=3000, ff=0, qtot<-0.5, r0=0.200,
                    Dlist, qextra, nw, xc, yc, links,
                    nwextra, xextra, yextra, xstart)

# Case B
Dlist[1:6] <- 0.4
result <- wellsolve(projected, T=0.0025 ,R=3000, ff=0.03, qtot<-0.5, r0=0.200,
                    Dlist, qextra, nw, xc, yc, links,
                    nwextra, xextra, yextra, xstart)

# Case C
qextra[1] <- 0.1
qextra[2] <- 0.4
result <- wellsolve(projected, T=0.0025 ,R=3000, ff=0.03, qtot<-0.5, r0=0.200,
                    Dlist, qextra, nw, xc, yc, links,
                    nwextra, xextra, yextra, xstart)