Hydraulic Tomography Estimates Improved by Zonal Information From the Clustering of Geophysical Survey Data

Abstract Hydraulic tomography (HT) has been demonstrated as a robust approach to characterize subsurface heterogeneity through the inverse modeling of multiple pumping data. However, smooth or even erroneous tomograms can result when insufficient observations are involved in the inversion. In this study, the feasibility of integrating geophysical survey data into HT analysis is investigated. First, k ‐means clustering is utilized to extract zonal information from borehole geophysical logs, and a new type of spatial constraints containing geological knowledge is proposed to obtain improved hydrostratigraphic boundaries along boreholes. Next, zonation models are constructed by applying clustering‐based zone geometry and populating zonal estimates of hydraulic conductivity ( K ) from analyzing pumping data. Afterwards, zonation models are treated as the initial guess of spatial variability in the geostatistical inversion of HT analysis. Additionally, local K measurements can be utilized to further improve HT estimates. Comparative cases of HT analyses are designed for a numerical sandbox experiment to highlight the HT performance integrated with geophysical surveys, in which the geostatistical inversion is initialized with: (a) a homogeneous K field; (b) zonation models built by the clustering of disparate geophysical surveys with/without spatial constraints; and (c) zonation improved by incorporating local K measurements. Based on ln K field comparisons and validation through predictions of drawdowns and tracer plume migration from independent tests not used in the calibration effort, we find that integration of geophysical surveys into HT analysis by clustering with spatial constraints is demonstrated as an effective approach, and local K measurements can further improve HT estimates.