blended_head_change: Calculate the transmissivity weighted blended head
In jkennel/sokol: Blended water level analysis for wells

Description Usage Arguments Details Value References See Also Examples

View source: R/blended_head_change.R

'estimate_blended_head' takes a vector of transmissivity and a vector of hydraulic head and calculates the expected blended head for a cross-connected well. The length units for transmissivity and hydraulic head should be the same.

1	blended_head_change(transmissivity, head_1, head_2)

`transmissivity`	the transmissivity distribution for each interval
`head_1`	the initial formation head distribution for each interval
`head_2`	the updated formation head distribution for each interval

'estimate_blended_head' uses the following to calculate the blended head for each head distribution based on Sokol's 1963 equation 6:

h_{blended} = \frac{∑\limits_{i=1}^{N} T_i h_i}{∑\limits_{i=1}^{N} T_i}

where:

* h_{blended} is the blended head in the well * T_i is the transmissivity of the i^th interval * h_i is the formation head of the i^th interval * i is the interval number * N is the total number of intervals

a single value of the change in blended head

Sokol, D. (1963). Position and fluctuations of water level in wells perforated in more than one aquifer. Journal of Geophysical Research, 68(4), 1079-1080.

Other blended head methods: estimate_blended_head(), estimate_formation_head(), estimate_missing(), plot_blended()

transmissivity <- c(1, 1)           # random T distribution
head_1 <- c(1,3)                    # decreasing head with depth
head_2 <- c(1,5)                    # decreasing head with depth

blended_head_change(transmissivity, head_1, head_2)