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: |
Ohm's law
Mathematical optimization Partial differential equation 010504 meteorology & atmospheric sciences [SDE.MCG]Environmental Sciences/Global Changes Mathematical analysis 010502 geochemistry & geophysics 01 natural sciences Finite element method Streaming current Tikhonov regularization symbols.namesake symbols Tensor Electrical resistivity tomography Computers in Earth Sciences Poisson's equation [SDU.ENVI]Sciences of the Universe [physics]/Continental interfaces environment 0105 earth and related environmental sciences Information Systems Mathematics |
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 |
Externí odkaz: |