Get capture zone

Get the capture zone for one or more wells

get_capture_zone(wells, aquifer, t_max = 365, wIDs = "all",
  n_particles = 8, buff_m = 20, injection_wells = FALSE, ...)

Arguments

wells	Wells data.frame object, containing well images
aquifer	Aquifer as an aquifer object, with Ksat, porosity, n, and boundaries (which are used to calculate gridded velocities)
t_max	number of days to run particle tracking
wIDs	wells at which to generate capture zones. Set to "all" or numeric value (or vector) containing wID for wells.
n_particles	number of particles to initialize at each well
buff_m	buffer radius (m) around each well at which particles are initialized
injection_wells	if TRUE, particle tracking from injection wells is allowed. if FALSE, particle tracking from injection wells is prohibited.
...	Additional arguments passed to track_particles

Details

Tracking particles are initialized radially around each well at a distance of buff_m. These particles must be initialized outside any grid cells that overlap the well, because particle velocities inside this cell will be incorrect.

Note: get_capture_zone does not work with recharge_type == "D".

Examples

bounds_df <- data.frame(bound_type=c("NF","NF","CH","NF"),m=c(Inf,0,Inf,0),b=c(0,1000,1000,0))
aquifer <- define_aquifer(aquifer_type="confined",Ksat=0.001,n=0.4,h0=0,z0=20,bounds=bounds_df)
wells_df_orig <- wells_example
wells_df_orig[4,"Q"] <- 0.25
wells <- generate_image_wells(define_wells(wells_df_orig), aquifer)
particle_paths <- get_capture_zone(wells, aquifer, t_max = 365*10, wIDs = "all", n_particles = 4)
particle_endpoint <- particle_paths[particle_paths$endpoint,]

library(ggplot2)
ggplot() +
  geom_segment(data=aquifer$bounds,aes(x1,y1,xend=x2,yend=y2,linetype=bound_type)) +
  geom_path(data=particle_paths,aes(x,y,color=as.factor(wID),group=interaction(pID,wID))) +
  geom_point(data=wells[wells$wID==wells$orig_wID,],aes(x,y,color=as.factor(wID)),size=3) +
  geom_point(data=particle_endpoint,aes(x,y,shape=status)) +
  scale_shape_manual(values=c(5,4,3,0,1)) +
  coord_equal()