Articles | Volume 11, issue 5
Hydrol. Earth Syst. Sci., 11, 1661–1671, 2007
https://doi.org/10.5194/hess-11-1661-2007
Hydrol. Earth Syst. Sci., 11, 1661–1671, 2007
https://doi.org/10.5194/hess-11-1661-2007

  16 Oct 2007

16 Oct 2007

Forward Modeling and validation of a new formulation to compute self-potential signals associated with ground water flow

A. Bolève1,3, A. Revil1,2, F. Janod3, J. L. Mattiuzzo3, and A. Jardani2,4 A. Bolève et al.
  • 1CNRS- LGIT (UMR 5559), University of Savoie, Equipe Volcan, Chambéry, France
  • 2Colorado School of Mines, Dept. of Geophysics, 1500 Illinois street, Golden, CO, 80401, USA
  • 3SOBESOL, Savoie Technolac, BP 230, F-73375 Le Bourget du Lac Cedex, France
  • 4CNRS, University of Rouen, Département de Géologie, Rouen, France

Abstract. The classical formulation of the coupled hydroelectrical flow in porous media is based on a linear formulation of two coupled constitutive equations for the electrical current density and the seepage velocity of the water phase and obeying Onsager's reciprocity. This formulation shows that the streaming current density is controlled by the gradient of the fluid pressure of the water phase and a streaming current coupling coefficient that depends on the so-called zeta potential. Recently a new formulation has been introduced in which the streaming current density is directly connected to the seepage velocity of the water phase and to the excess of electrical charge per unit pore volume in the porous material. The advantages of this formulation are numerous. First this new formulation is more intuitive not only in terms of establishing a constitutive equation for the generalized Ohm's law but also in specifying boundary conditions for the influence of the flow field upon the streaming potential. With the new formulation, the streaming potential coupling coefficient shows a decrease of its magnitude with permeability in agreement with published results. The new formulation has been extended in the inertial laminar flow regime and to unsaturated conditions with applications to the vadose zone. This formulation is suitable to model self-potential signals in the field. We investigate infiltration of water from an agricultural ditch, vertical infiltration of water into a sinkhole, and preferential horizontal flow of ground water in a paleochannel. For the three cases reported in the present study, a good match is obtained between finite element simulations performed and field observations. Thus, this formulation could be useful for the inverse mapping of the geometry of groundwater flow from self-potential field measurements.