RESINVM3D: A 3D resistivity inversion package

Abstract

We have developed an open source 3D, MATLAB based, resistivity inversion package. The forward solution to the governing partial differential equation is efficiently computed using a second-order finite volume discretization coupled with a preconditioned, biconjugate, stabilized gradient algorithm. Using the analytical solution to a potential field in a homogeneous half space, we evaluate the accuracy of our numerical forward solution and, subsequently, develop a source correction factor that reduces forward modeling errors associated with boundary effects and source electrode singularities. For the inversion algorithm we have implemented an inexact Gauss-Newton solver, with the model update being calculated using a preconditioned conjugate gradient algorithm. The inversion uses a combination of zero and first order Tikhonov regularization. Two synthetic examples demonstrate the usefulness of this code. The first example considers a surface resistivity survey with 3813 measurements. The discretized model space contains 19,040 cells. For this example, the inversion package converges in approximately [Formula: see text] on a [Formula: see text] Pentium 4, with [Formula: see text] of RAM. The second example considers the case of borehole based data acquisition. For this example there were 4704 measurements and 13,200 model cells. The inversion for this example requires [Formula: see text] of computational time.

Covariance-constrained difference inversion of time-lapse electrical resistivity tomography data

“…However, it has to be noted that a uniform distribution is inconsistent with the choice of an l 2 -norm for the calculation of the data misfit, which assumes a Gaussian distribution of the error. Still, this is often neglected in practice (Pidlisecky et al, 2007). In this case, the effect is expected to be similar for SCI and CCI.…”

Image Fusion for Travel Time Tomography Inversion

“…Using the abnormal regional model and programming ideas in reference [19], modifying the program according to improved Wilkinson iteration algorithm. Regional location, launch or receiving point coordinates are given in Figure 5.…”

“…High frequency coefficient: (19) Among them, if ; otherwise Low frequency coefficient: (20) Where, if ; otherwise .…”

Geophysical Methods for Monitoring Temperature Changes in Shallow Low Enthalpy Geothermal Systems

“…The subsurface is divided in cells in the three dimensions and the model is inverted considering current flow in all directions [69]. When the collected data results from electrodes situated on a plane, the inversion generally uses a 2.5D scheme, where the resistivity distribution is considered as constant in the direction perpendicular to the section [70].…”