| Abstract: | Numerical computations using the integral equation method are presented for resistivity and IP responses due to arbitrarily shaped 3-dimensional bodies in a layered earth. The unknown surface charge density distribution is expressed as the solution of Fredholm's integral equation of the second kind. Use of moment method (with pulse basis function and point-collocation) yields the matrix equations for the unknowns. The contributions to Green's function are solved (a) analytically for the primary and (b) by convolution for the secondary contributions resulting in a fast algorithm. The further step of computing potential, apparent resistivity, chargeability etc., for any electrode system, is straightforward. Our results show a good agreement with those from finite difference methods and physical tank experiments. The CPU time is only 138 s on a super-minicomputer for an apparent resistivity pseudo-section, even with 96 elementary cells as used for discretization. A large number of models for different geological situations were studied; some are presented here. |
