Abstract The numerical solution of Green s function for the potential in 2-D arbitrary inhomo-geneous media with axial symmetry has been given by use of efficient half-analytical, half-numerical hybrid method. Then the logging responses of various kinds of the DC resistivity log with axisym-metric excitation have been obtained by using surface integral equation method to match the boundary conditions on the electrodes of the logging sonde. Comparing the results with that obtained by other methods, one can see good p...