Incorporating geologic information into hydraulic tomography: A general framework based on geostatistical approach

Abstract

Hydraulic tomography (HT) has become a mature aquifer test technology over the last two decades. It collects nonredundant information of aquifer heterogeneity by sequentially stressing the aquifer at different wells and collecting aquifer responses at other wells during each stress. The collected information is then interpreted by inverse models. Among these models, the geostatistical approaches, built upon the Bayesian framework, first conceptualize hydraulic properties to be estimated as random fields, which are characterized by means and covariance functions. They then use the spatial statistics as prior information with the aquifer response data to estimate the spatial distribution of the hydraulic properties at a site. Since the spatial statistics describe the generic spatial structures of the geologic media at the site rather than site‐specific ones (e.g., known spatial distributions of facies, faults, or paleochannels), the estimates are often not optimal. To improve the estimates, we introduce a general statistical framework, which allows the inclusion of site‐specific spatial patterns of geologic features. Subsequently, we test this approach with synthetic numerical experiments. Results show that this approach, using conditional mean and covariance that reflect site‐specific large‐scale geologic features, indeed improves the HT estimates. Afterward, this approach is applied to HT surveys at a kilometer‐scale‐fractured granite field site with a distinct fault zone. We find that by including fault information from outcrops and boreholes for HT analysis, the estimated hydraulic properties are improved. The improved estimates subsequently lead to better prediction of flow during a different pumping test at the site.

A Reduced‐Order Successive Linear Estimator for Geostatistical Inversion and its Application in Hydraulic Tomography

“…Although SLE and QLGA share similarities, a major difference between them is that the covariance of unknown parameters in SLE is updated to derive the residual covariance every iteration (Nowak & Cirpka, ; Zha et al, ). To be specific, the collected data are assimilated iteratively with a stochastic linear estimator in SLE to consider the nonlinear relationship between data and parameter as well as the covariance of the parameters.…”

Fusion of Time‐Lapse Gravity Survey and Hydraulic Tomography for Estimating Spatially Varying Hydraulic Conductivity and Specific Yield Fields

“…To be realistic, a three‐dimension variably saturated flow model, which considers the unsaturated hydraulic properties in heterogeneous unconfined aquifers, should be used as advocated by Yeh et al (). Most importantly, the proposed methodology must be tested and validated in the fields as have done by Zha et al (, 2016), Illman et al (), Berg and Illman (), and Zha et al (). In spite of these limitations, the study demonstrates that the proposed method is a promising step toward the development of a useful data fusion concept and method to probe subsurface and to better characterize the spatial variability of hydrogeological parameters in unconfined aquifers.…”

A Non-Uniform Broadcast Fertilization Method and Its Performance Analysis under Basin Irrigation

“…The mean absolute relative error (ARE) [45] is used to evaluate the simulation effect of solute transport in surface water flow between the measured solute concentration and the simulated value. ARE is calculated as,…”