R/estimate_blended_head.R
estimate_blended_head.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.
estimate_blended_head(transmissivity, head)
| transmissivity | the transmissivity distribution for each interval |
|---|---|
| head | the head distribution for each interval |
a single value of the blended head
`estimate_blended_head` uses the following to calculate the blended head 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:
blended_head_change(),
estimate_formation_head(),
estimate_missing(),
plot_blended()
transmissivity <- abs(rnorm(10, sd = 100)) # random T distribution head <- sort(rnorm(10)) # decreasing head with depth estimate_blended_head(transmissivity, head)#> [1] -0.4531166