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

被引：62

作者：

Karaoulis, M. ^{[1
]}

Revil, A. ^{[1
,2
]}

Tsourlos, P. ^{[3
]}

Werkema, D. D. ^{[4
]}

Minsley, B. J. ^{[5
]}

机构：

[1] Colorado Sch Mines, Dept Geophys, Golden, CO 80401 USA

[2] Univ Savoie, Equipe Volcans, INSU CNRS LGIT UMR 5559, Le Bourget Du Lac, France

[3] Aristotle Univ Thessaloniki, Dept Geophys, GR-54006 Thessaloniki, Greece

[4] US EPA, ORD, NERL, ESD,CMB, Las Vegas, NV 89193 USA

[5] USGS, Crustal Geophys & Geochem Sci Ctr, Denver Fed Ctr, Denver, CO 80225 USA

来源：

COMPUTERS & GEOSCIENCES | 2013年 / 54卷

关键词：

Resistivity; Induced polarization; Inversion; Time-lapse; SPECTRAL INDUCED POLARIZATION; SHALY SANDS; ELECTRICAL CONDUCTIVITIES; OCCAMS INVERSION;

D O I：

10.1016/j.cageo.2013.01.008

中图分类号：

TP39 [计算机的应用];

学科分类号：

081203 ; 0835 ;

摘要：

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 (R) code package and a graphical user interface (GUI). (C) 2013 Elsevier Ltd. All rights reserved.

引用

页码：164 / 170

页数：7

共 40 条

[11] 4D active time constrained resistivity inversion [J].

Karaoulis, M. C. ;

Kim, J. -H. ;

Tsourlos, P. I. .

JOURNAL OF APPLIED GEOPHYSICS, 2011, 73 (01) :25-34

[12]

Kemna A., 2000, TOMOGRAPHIC INVERSIO

[13]

Kim J.H., 2005, Proc. of Near Surface 2005, the 11th Euro. Meeting of EEGS, pA006

[14] 4-D inversion of DC resistivity monitoring data acquired over a dynamically changing earth model [J].

Kim, Jung-Ho ;

Yi, Myeong-Jong ;

Park, Sam-Gyu ;

Kim, Jae Gon .

JOURNAL OF APPLIED GEOPHYSICS, 2009, 68 (04) :522-532

[15]

LaBrecque D.J., 2009, Journal of Environmental and Engineering Geophysics, DOI DOI 10.4133/JEEG6.2.83

[16] The effects of noise on Occam's inversion of resistivity tomography data [J].

LaBrecque, DJ ;

Miletto, M ;

Daily, W ;

Ramirez, A ;

Owen, E .

GEOPHYSICS, 1996, 61 (02) :538-548

[17] A case study of resistivity and self-potential signatures of hydrothermal instabilities, Inferno Crater Lake, Waimangu, New Zealand [J].

Legaz, A. ;

Vandemeulebrouck, J. ;

Revil, A. ;

Kemna, A. ;

Hurst, A. W. ;

Reeves, R. ;

Papasin, R. .

GEOPHYSICAL RESEARCH LETTERS, 2009, 36

[18] Inversion of 2D spectral induced polarization imaging data [J].

Loke, MH ;

Chambers, JE ;

Ogilvy, RD .

GEOPHYSICAL PROSPECTING, 2006, 54 (03) :287-301

[19] LEAST-SQUARES DECONVOLUTION OF APPARENT RESISTIVITY PSEUDOSECTIONS [J].

LOKE, MH ;

BARKER, RD .

GEOPHYSICS, 1995, 60 (06) :1682-1690

[20] Application of time-lapse ERT imaging to watershed characterization [J].

Miller, Carlyle R. ;

Routh, Partha S. ;

Brosten, Troy R. ;

McNamara, James P. .

GEOPHYSICS, 2008, 73 (03) :G7-G17

← 1 2 3 4 →