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

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

Tomography of the Darcy velocity from self-potential measurements

An algorithm is developed to interpret self-potential (SP) data in terms of distribution of Darcy velocity of the ground water. The model is based on the proportionality existing between the streaming current density and the Darcy velocity. Because the inverse problem of current density determination from SP data is underdetermined, we use Tikhonov regularization with a smoothness constraint based on the differential Laplacian operator and a prior model. The regularization parameter is determined by the L-shape method. The distribution of the Darcy velocity depends on the localization and number of non-polarizing electrodes and information relative to the distribution of the electrical resistivity of the ground. A priori hydraulic information can be introduced in the inverse problem. This approach is tested on two synthetic cases and on real SP data resulting from infiltration of water from a ditch

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

International audienceThe 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

A new formulation to compute self-potential signals associated with ground water flow

International audienceThe 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 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 is also easily extendable to non-viscous laminar flow problems (high Reynolds number ground water flow in cracks for example) and to unsaturated conditions with applications to the vadose zone. We demonstrate here that 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 the finite element simulations performed with the finite element code Comsol Multiphysics 3.3 and field observations. Finally, this formulation seems also very promising for the inversion of the geometry of ground water flow from the monitoring of self-potential signals

Application of acoustic noise and self-potential localization techniques to a buried hydrothermal vent (Waimangu Old Geyser site, New Zealand)

International audienceA seismo-acoustic and self-potential survey has been performed in the hydrothermal area of the old Waimangu Geyser (New Zealand), which was violently erupting a century ago. Nowadays, no surface activity is visible there. We set-up an array of 16 geophones and recorded a high and steady acoustic ambient noise. We applied the matched field processing (MFP) approach to the acoustic data to locate the sources responsible for the ambient noise. The white noise constraint processor reveals the presence of a unique and well-focused acoustic source at a depth of 1.5 m below the seismic array. For this very shallow source, the application of MFP enabled the determination of both the source location and the dispersion curve of seismic velocity. The study was completed by self-potential (SP) measurements on several profiles around the acoustic noise source, which displayed a large positive anomaly above it. The results of the SP inversion gave an electric streaming current density source very close to the acoustic one. Both sources likely belong to a shallow hydrothermal structure interpreted as a small convective cell of boiling water beneath an impermeable layer. The joint application of these methods is a promising technique to recognize hydrothermal structures and to study their dynamics