Hydraulic conductivity imaging from 3‐D transient hydraulic tomography at several pumping/observation densities

Abstract

3‐D Hydraulic tomography (3‐D HT) is a method for aquifer characterization whereby the 3‐D spatial distribution of aquifer flow parameters (primarily hydraulic conductivity, K) is estimated by joint inversion of head change data from multiple partially penetrating pumping tests. While performance of 3‐D HT has been studied extensively in numerical experiments, few field studies have demonstrated the real‐world performance of 3‐D HT. Here we report on a 3‐D transient hydraulic tomography (3‐D THT) field experiment at the Boise Hydrogeophysical Research Site which is different from prior approaches in that it represents a “baseline” analysis of 3‐D THT performance using only a single arrangement of a central pumping well and five observation wells with nearly complete pumping and observation coverage at 1 m intervals. We jointly analyze all pumping tests using a geostatistical approach based on the quasi‐linear estimator of Kitanidis (1995). We reanalyze the system after progressively removing pumping and/or observation intervals; significant progressive loss of information about heterogeneity is quantified as reduced variance of the K field overall, reduced correlation with slug test K estimates at wells, and reduced ability to accurately predict independent pumping tests. We verify that imaging accuracy is strongly improved by pumping and observational densities comparable to the aquifer heterogeneity geostatistical correlation lengths. Discrepancies between K profiles at wells, as obtained from HT and slug tests, are greatest at the tops and bottoms of wells where HT observation coverage was lacking.