3-D time-domain induced polarization tomography: a new approach based on a source current density formulation. Issue 1 (18th December 2017)

• Record Type: Journal Article
• Authors: Soueid Ahmed, A; Revil, A
• Is Part Of: Geophysical journal international. Volume 213:Issue 1(2018:Apr.)
• Page Start: 244
• Page End: 260
• Publication Date: 2017-12-18
• Subjects: Electrical properties -- Hydrogeophysics -- Electrical resistivity tomography (ERT) -- Inverse theory -- Numerical modelling -- Tomography; Geophysics -- Periodicals; 550
• DOI: 10.1093/gji/ggx547
• Languages: English
• ISSNs: 0956-540X

Abstract

Summary: Induced polarization (IP) of porous rocks can be associated with a secondary source current density, which is proportional to both the intrinsic chargeability and the primary (applied) current density. This gives the possibility of reformulating the time domain induced polarization (TDIP) problem as a time-dependent self-potential-type problem. This new approach implies a change of strategy regarding data acquisition and inversion, allowing major time savings for both. For inverting TDIP data, we first retrieve the electrical resistivity distribution. Then, we use this electrical resistivity distribution to reconstruct the primary current density during the injection/retrieval of the (primary) current between the current electrodes A and B. The time-lapse secondary source current density distribution is determined given the primary source current density and a distribution of chargeability (forward modelling step). The inverse problem is linear between the secondary voltages (measured at all the electrodes) and the computed secondary source current density. A kernel matrix relating the secondary observed voltages data to the source current density model is computed once (using the electrical conductivity distribution), and then used throughout the inversion process. This recovered source current density model is in turn used to estimate the time-dependent chargeability (normalized voltages) in each cell of the domain of interest. Assuming a Cole–Cole model for simplicity, we can reconstruct the 3-D distributions of the relaxation time τ and the Cole–Cole exponent c by fitting the intrinsic chargeability decay curve to a Cole–Cole relaxation model for each cell. Two simple cases are studied in details to explain this new approach. In the first case, we estimate the Cole–Cole parameters as well as the source current density field from a synthetic TDIP data set. Our approach is successfully able to reveal the presence of the anomaly and to invert its Cole–Cole parameters. In the second case, we perform a laboratory sandbox experiment in which we mix a volume of burning coal and sand. The algorithm is able to localize the burning coal both in terms of electrical conductivity and chargeability.