SP2DINV: A 2 D forward and inverse code for streaming potential problems

Autor: A. Soueid Ahmed, Jean-Paul Dupont, Abderrahim Jardani, André Revil
Přispěvatelé: Morphodynamique Continentale et Côtière (M2C), Université de Caen Normandie (UNICAEN), Normandie Université (NU)-Normandie Université (NU)-Institut national des sciences de l'Univers (INSU - CNRS)-Université de Rouen Normandie (UNIROUEN), Normandie Université (NU)-Centre National de la Recherche Scientifique (CNRS), Department of Geophysics [Golden CO], Colorado School of Mines
Jazyk: angličtina
Rok vydání: 2013
Předmět:
Zdroj: Computers & Geosciences
Computers & Geosciences, Elsevier, 2013, 59, pp.9-16. ⟨10.1016/j.cageo.2013.05.008⟩
ISSN: 0098-3004
1873-7803
DOI: 10.1016/j.cageo.2013.05.008⟩
Popis: The self-potential method corresponds to the passive measurement of the electrical field in response to the occurrence of natural sources of current in the ground. One of these sources corresponds to the streaming current associated with the flow of the ground water. We can therefore apply the self-potential method to recover non-intrusively some information regarding the ground water flow. We first solve the forward problem starting with the solution of the ground water flow problem, then computing the source current density, and finally solving a Poisson equation for the electrical potential. We use the finite-element method to solve the relevant partial differential equations. In order to reduce the number of (petrophysical) model parameters required to solve the forward problem, we introduced an effective charge density tensor of the pore water, which can be determined directly from the permeability tensor for neutral pore waters. The second aspect of our work concerns the inversion of the self-potential data using Tikhonov regularization with smoothness and weighting depth constraints. This approach accounts for the distribution of the electrical resistivity, which can be independently and approximately determined from electrical resistivity tomography. A numerical code, SP2DINV, has been implemented in Matlab to perform both the forward and inverse modeling. Three synthetic case studies are discussed.
Databáze: OpenAIRE