Several graphical, and numerical methods have been developed to interpret SP anomalies including curve matching, characteristic points, least-squares, derivative and gradient analysis, nonlinear modeling, and Fourier analysis techniques. Examples of such techniques used for interpreting the self-potential anomalies of horizontal and vertical cylinders, spheres, sheets and geological contacts are classified as follows.

(1) Methods using only a few points on the anomaly curve. These were originally developed by DeWitte (1948), Yüngül (1950), Paul (1965), Paul et al. (1965), Bhattacharya and Roy (1981), Atchuta Rao and Ram Babu (1983). The essential disadvantage of these methods was related to the fact that only a few points are used on the anomaly curve, and hence, the interpreted results are not reliable.

(2) Curve matching techniques. According to Meiser (1962), Satyanarayana Murty and Haricharen (1985), the field curve is compared with sets of theoretical curves either manually or using a computer. This process is cumbersome and the complexity of the method is very high especially when the variables are numerous.

(3) Least-squares methods. Here, the model parameters that give a best fit are derived using an initial presumption, as well as characteristic points and distances from the measured anomaly curve. Examples of these methods are of Abdelrahman and Sharafeldin (1997), [Abdelrahman et al., 1997b] and [Abdelrahman et al., 2003].

The above-mentioned three categories have the disadvantage of using characteristic points, distances, curves, and nomograms for interpretation, which is subject to human error in estimating these few points and distances, especially if they are derived by interpolation procedure. This can consequently lead to errors in derived parameters. Moreover, some of these methods demand knowledge of the anomalous body shape, and based on the interpreter's choice, the results could vary a lot.

(4) Methods using derivative analysis and gradients. Namely those of [Abdelrahman et al., 1997a],[Abdelrahman et al., 1998a], [Abdelrahman et al., 1998b] and [Abdelrahman et al., 2003] belong to this category.

(5) Methods using Fourier analysis and the wave number domain. Particularly those of Atchuta Rao et al. (1982) and Roy and Mohan (1984) belong to this group.

(6) Modeling and Inversion Methods such as those given by Guptasarma (1983), Furness (1992), and Shi and Morgan (1996).

The last three categories are greatly influenced by noise in measured data and can lead to serious errors.

ER -