Modeling techniques commonly exhibit errors of 3 to 10 percent or more in the calculation of apparent resistivities over earth models for which analytic solutions are easily available. A singularity occurs in the solution of any elliptic partial differential equation for which the forcing function is not smooth. The inability to adequately represent in discrete space a discontinuous function (in this case, the delta function describing the introduction of current at a point) commonly results in numerical error near the source of a modeled singularity.Inspection of an integrated finite-difference method for modeling the dc resistivity geophysical technique indicates much of the error encountered is of singular origin. A procedure is herein detailed by which the singularity is mathematically removed from the modeling process and reintroduced as a last step, thus preventing it from contributing to the numerical error. Using this procedure, the average error in apparent resistivity values for a model of a polar-dipole traverse over a nonconducting sphere is reduced by 40 percent. For a dipole-dipole traverse of a two-layer model the error decreases by 75 percent, and in the case of a Wenner profile of a model of a vertically faulted earth, the average error is diminished by 90 percent.

This content is PDF only. Please click on the PDF icon to access.

First Page Preview

First page PDF preview
You do not currently have access to this article.