IP4DI: A software for time-lapse 2D/3D DC-resistivity and induced polarization tomography

详细信息    查看全文

• 作者：M. Karaoulis^a ; ^{marios.karaoulis@gmail.com} ; A. Revil^a ; ^b ; ^{arevil@mines.edu} ; P. Tsourlos^c ; ^{tsourlos@geo.auth.gr} ; D.D. Werkema^d ; ^{werkema.d@epamail.epa.gov} ; B.J. Minsley^e ; ^{bminsley@usgs.gov}
• 关键词：Resistivity ; Induced polarization ; Inversion ; Time-lapse
• 刊名：Computers and Geosciences
• 出版年：2013
• 出版时间：April, 2013
• 年：2013
• 卷：54
• 期：Complete
• 页码：164-170
• 全文大小：771 K

文摘

We propose a 2D/3D forward modelling and inversion package to invert direct current (DC)-resistivity, time-domain induced polarization (TDIP), and frequency-domain induced polarization (FDIP) data. Each cell used for the discretization of the 2D/3D problems is characterized by a DC-resistivity value and a chargeability or complex conductivity for TDIP/FDIP problems, respectively. The governing elliptic partial differential equations are solved with the finite element method, which can be applied for both real and complex numbers. The inversion can be performed either for a single snapshot of data or for a sequence of snapshots in order to monitor a dynamic process such as a salt tracer test. For the time-lapse inversion, we have developed an active time constrained (ATC) approach that is very efficient in filtering out noise in the data that is not correlated over time. The forward algorithm is benchmarked with simple analytical solutions. The inversion package IP4DI is benchmarked with three tests, two including simple geometries. The last one corresponds to a time-lapse resistivity problem for cross-well tomography during enhanced oil recovery. The algorithms are based on MATLAB^? code package and a graphical user interface (GUI).