Specific storage and hydraulic conductivity tomography through the joint inversion of hydraulic heads and self-potential data

• Published in: Advances in water resources Vol. 89; pp. 80 - 90
• Main Authors: Ahmed, A. Soueid, Jardani, A., Revil, A., Dupont, J.P.
• Format: Journal Article
• Language: English
• Published: Elsevier Ltd 01.03.2016; Elsevier
• Subjects: Aquifers; Computational fluid dynamics; Earth Sciences; Fluid flow; Geophysics; Hydraulics; Inverse problems; Inversion; Inversions; Joint inversion; Pumping tests; Resistivity; Sciences of the Universe; Self-potential; Tomography; Transient hydraulic tomography; Inversion; Self-potential; Pumping tests; Joint inversion; Transient hydraulic tomography
• ISSN: 0309-1708; 1872-9657
• DOI: 10.1016/j.advwatres.2016.01.006

•Self-potential data can be jointly inverted with head data.•The two types of data have complementary sensitivity map.•The reconstruction of the hydraulic conductivity and specific storage is improved.•The low-cost of the self-potential method makes it attractive for pumping tests. Transient hydraulic tomography is used to image the heterogeneous hydraulic conductivity and specific storage fields of shallow aquifers using time series of hydraulic head data. Such ill-posed and non-unique inverse problem can be regularized using some spatial geostatistical characteristic of the two fields. In addition to hydraulic heads changes, the flow of water, during pumping tests, generates an electrical field of electrokinetic nature. These electrical field fluctuations can be passively recorded at the ground surface using a network of non-polarizing electrodes connected to a high impedance (> 10MOhm) and sensitive (0.1mV) voltmeter, a method known in geophysics as the self-potential method. We perform a joint inversion of the self-potential and hydraulic head data to image the hydraulic conductivity and specific storage fields. We work on a 3D synthetic confined aquifer and we use the adjoint state method to compute the sensitivities of the hydraulic parameters to the hydraulic head and self-potential data in both steady-state and transient conditions. The inverse problem is solved using the geostatistical quasi-linear algorithm framework of Kitanidis. When the number of piezometers is small, the record of the transient self-potential signals provides useful information to characterize the hydraulic conductivity and specific storage fields. These results show that the self-potential method reveals the heterogeneities of some areas of the aquifer, which could not been captured by the tomography based on the hydraulic heads alone. In our analysis, the improvement on the hydraulic conductivity and specific storage estimations were based on perfect knowledge of electrical resistivity field. This implies that electrical resistivity will need to be jointly inverted with the hydraulic parameters in future studies and the impact of its uncertainty assessed with respect to the final tomograms of the hydraulic parameters.