Transient Hydraulic Tomography Analysis of Fourteen Pumping Tests at a Highly Heterogeneous Multiple Aquifer–Aquitard System
Illman, Walter A.
Yeh, Tian-Chyi Jim
Mok, Chin Man Bill
Berg, Steven J.
AffiliationUniv Arizona, Dept Hydrol & Atmospher Sci
geostatistical inverse modeling
MetadataShow full item record
CitationZhao, Z., Illman, W. A., Zha, Y., Yeh, T. C. J., Mok, C. M. B., Berg, S. J., & Han, D. (2019). Transient Hydraulic Tomography Analysis of Fourteen Pumping Tests at a Highly Heterogeneous Multiple Aquifer–Aquitard System. Water, 11(9), 1864.
RightsCopyright © 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Collection InformationThis item from the UA Faculty Publications collection is made available by the University of Arizona with support from the University of Arizona Libraries. If you have questions, please contact us at email@example.com.
AbstractHydraulic tomography based on geostatistics has proven to be robust in characterizing subsurface heterogeneity in hydraulic conductivity (K) and specific storage (S-s) through the joint inversion of drawdown records from multiple pumping tests. However, the spatially variable estimates can be smooth or even erroneous for areas where pumping/observation data densities are not high. Previous hydraulic tomography surveys conducted at the North Campus Research Site (NCRS) on the University of Waterloo campus in Waterloo, Canada, revealed that the estimated hydraulic parameters were smooth and the known aquitard was erroneously identified as a high K zone. This was likely the consequence of the site being highly heterogeneous, while only utilizing four pumping tests and not having measurable drawdowns in the low K aquitard for inverse modeling. Here, we investigate whether improved K and S-s estimates could be obtained through the inclusion of additional pumping test data by stressing both aquifer and aquitard zones for a sufficiently long period. Specifically, six additional pumping/injection tests were conducted at the site, and a transient hydraulic tomography analysis with 14 tests was completed. Results reveal that there is a significant improvement to the K and S-s tomograms in terms of the visual correspondence with various geologic units, including its connectivity. More importantly, with the availability of additional data, we found that the inverse model now can better capture the high and low K features for nine boreholes when compared with K values obtained from permeameter tests. The estimated K and S-s tomograms are then used for the forward simulation of one additional pumping test not used for model calibration, revealing reasonable predictions. While encouraging results are obtained by including a large number of pumping tests to the transient hydraulic tomography analysis, stratigraphic boundaries are still smoothed, which is a direct consequence of utilizing a geostatistics-based inversion approach that assumes stationarity in statistical properties. To capture such sharp boundaries, incorporation of additional data types, such as geological and geophysical information, may be necessary when data densities are not sufficiently high.
NoteOpen access journal
VersionFinal published version
SponsorsEnvironmental Security and Technology Certification Program (ESTCP) [ER201212]; National Natural Science Foundation of China (NSFC)National Natural Science Foundation of China ; National Natural Science Foundation of ChinaNational Natural Science Foundation of China