Hydraulic conductivity field characterization from the joint inversion of hydraulic heads and self- potential data

Pumping tests can be used to estimate the hydraulic conductivity field from the inversion of hydraulic head data taken intrusively in a set of piezometers. Nevertheless, the inverse problem is strongly underdetermined. We propose to add more information by adding self-potential data taken at the ground surface during pumping tests. These self-potential data correspond to perturbations of the electrical field caused directly by the flow of the groundwater. The coupling is electrokinetic in nature that is due to the drag of the excess of electrical charges existing in the pore water. These self-potential signals can be easily measured in field conditions with a set of the nonpolarizing electrodes installed at the ground surface. We used the adjoint-state method for the estimation of the hydraulic conductivity field from measurements of both hydraulic heads and self potential during pumping tests. In addition, we use a recently developed petrophysical formulation of the streaming potential problem using an effective charge density of the pore water derived directly from the hydraulic conductivity. The geostatistical inverse framework is applied to five synthetic case studies with different number of wells and electrodes and thickness of the confining unit. To evaluate the benefits of incorporating the self-potential data in the inverse problem, we compare the cases in which the data are combined or not. Incorporating the self-potential information improves the estimate of hydraulic conductivity field in the case where the number of piezometers is limited. However, the uncertainty of the characterization of the hydraulic conductivity from the inversion of the self-potential data is dependent on the quality of the distribution of the electrical conductivity used to solve the Poisson equation. Consequently, the approach discussed in this paper requires a precise estimate of the electrical conductivity distribution of the subsurface and requires therefore new strategies to be developed for the joint inversion of the hydraulic and electrical conductivity distributions.