R/blended_head_change.R
blended_head_change.Rd`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.
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 |
a single value of the change in blended head
`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{\sum\limits_{i=1}^{N} T_i h_i}{\sum\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
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)#> [1] 1