A fast imaging method for the interpretation of self

Scientific Reports volume 13, Article number: 13548 (2023) Cite this article

248 Accesses

1 Altmetric

Metrics details

We describe a rapid imaging approach for the interpretation of self-potential data collected along profile by some geometrically simple model of cylinders and spheres. The approach calculates the correlation coefficient between the analytic signal (AS) of the observed self-potential measurements and the AS of the self-potential signature of the idealized model. The depth, electric dipole moment, polarization angle, and center are the inverse parameters we aim to extract from the imaging approach for the interpretative model, and they pertain to the highest value of the correlation coefficient. The approach is demonstrated on noise-free numerical experiments, and reproduced the true model parameters. The accuracy and stability of the proposed approach are examined on numerical experiments contaminated with realistic noise levels and regional fields prior to the interpretation of real data. Following that, five real field examples from geothermal systems and mineral exploration have been successfully analyzed. The results agree well with the published research.

Fox1 proposed the self-potential method using a copper electrode and galvanometer instrument to explore a copper-sulfide ore body at Cornwall, England. The self-potential method has advanced2,3,4,5,6,7,8,9,10,11,12,13,14,15 and been widely used for graphite, sulfide, magnetite, uranium and gold prospecting16,17,18,19,20,21,22, mapping paleo-shear zones23,24, archaeological investigations25, geotechnical engineering26, cave discovery27, coal fires detection28,29,30, and monitor water movement31,32,33. The electrical self-potential method has been applied to a broad range of monitoring studies, like landslides or mass movements that happened by cumulative pore pressure in the rock34.

Self-potential is a passive technique to measure self-potential differences that occur naturally in the earth’s subsurface14. The mechanism and origin of the self-potential anomalies have been discussed by several authors9,10,13,14,35,36,37,38,39. Self-potential methods are preferred over other geophysical techniques in measurements that are sensitive to fluid movements through fractured and porous rock, and under natural or an applied hydraulic gradient responding to weak fluid movements38,40. Several approaches for self-potential data forward modeling, inversion and interpretation have been developed10,13,41,42,43,44,45,46,47,48,49,50,51,52,53. These approaches can be grouped into two categories.

Geometry and parameters of the assumed source models. Top, middle and bottom panels present the sphere, semi-infinite vertical cylinder, and infinitely long horizontal cylinder models.

Flowchart showing the workflow of the developed scheme.

Class I is pertinent to the self-potential anomaly of multi-dimensional arbitrary structures including both the two (2D) and three-dimensional (3D) SP modeling and inversion. Multi-dimensional SP inversions can be non-unique, unstable and require large computational time54,55,56,57,58,59,60. Standard methods to retrieve a stable solution of an ill-posed inverse problem are the regularization techniques61.

Class II approximates the collected self-potential anomaly by some geometrically simple models, such as vertical cylinders, horizontal cylinders and spheres. This class offers a fast quantitative interpretation, and the objective is to infer the depth, location and polarization parameter of the interpretive model that best fits the observed data. The study we pursue here belongs to this class. Numerous quantitative methods (graphical and numerical) were established for Class II in order to obtain the shape, depth and polarization parameter of the causative source from the measured self-potential anomaly43,44,45,46,49,51,53,62,63,64,65. The disadvantage of these methods is that they are subjective, and consequently can result in some error in the model parameters66.

Model 1: SP data (noise free). (a) SP anomaly produced from a horizontal cylinder model. (b) Horizontal (HD) and vertical (VD) derivatives of the SP data rendered in (a). (c) Analytic signal (AS) amplitude. (d) Image of the R-parameter (R).

Model 1: noise-free data. Relationship of the R-parameter, shape factor, and depth.

Abdelrahman et al.63 introduced a graphical technique to infer the depth and shape of the buried structures from the second moving average residual self-potential anomalies. Santos67 applied the particle swarm optimization (PSO) scheme to invert the self-potential anomalies by some ideal geometric structures, like spheres, cylinders and inclined sheets. Mehanee66 developed a regularized scheme for the interpretation of self-potential data using the conjugate gradient minimizer in the space of logarithmed and non-logarithmed model parameters. Di Maio et al.68 presented a spectral analysis method for the interpretation of self-potential data by some geometrically simple model based upon the periodogram method (PM), multi Taper method (MTM), and maximum entropy method (MEM) to recover the depth of the anomalous body. Sungkono and Warnana69 applied the black hole algorithm (BHA) to self-potential data considering simple geometric bodies of sphere, horizontal cylinder and inclined sheet to determine the corresponding model parameters. The analytic signal (AS) amplitude of Nabighian70 can provide a key role in self-potential anomaly interpretation, and uses the spatial derivatives of the data, e.g.71,72,73. It is noted that the analytic signal method has been discussed numerously in the published literature for the interpretation of gravity and magnetic data, e.g.71,72,74,75,76,77,78,79,80,81,82.

An imaging methodology is presented in this paper for the interpretation of self-potential data by some idealized bodies (spheres and cylinders). The goal is to retrieve the self-potential profile’s origin point as well as the depth, angle of polarization, and shape index of the anomalous body. The technique estimates the correlation coefficient (the R-parameter) between the AS of the observed self-potential data and the AS of the numerical self-potential response of an assumed elucidative model. The favored elucidative model is that which attains the maximum correlation coefficient.

Model 2: Noisy data. (a) Noisy SP anomaly subject to interpretation. (b) Derivatives. (c) Amplitude. (d) Image.

Model 3: regional anomaly’s impact. (a) Composite SP anomaly of sphere model and first-order regional calculated by expression (6). (b) Derivatives. (c) Amplitude. (d) Image.

The self-potential imaging scheme presented here retrieves the horizontal location of the anomalous body as well as its depth, polarization angle and amplitude coefficient, and has three main benefits. First, the whole self-potential data profile is used when estimating the spatial parameters (horizontal location and depth) of the buried source, which are considered key information in geophysical prospecting. Second, the scheme uses an exact formula for the direct solution. Third, it does not demand a priori information about the subsurface resistivity distribution nor high computational resources. To our best knowledge, the R-parameter imaging method presented in this paper for the interpretation of self-potential data measured along profile by idealized models was not developed before. It is relevant to note that rigorous 3D inversion of self-potential data is computationally expensive, and requires a priori information for the model parameters (3D electric conductivity distribution) we invert for17.

The paper is structurally described as follows. The “Self-potential direct solution” section presents the direct problem (forward modeling solution). In “The method” section, the formulation of the proposed imaging scheme is explained. The “Numerical examples” section validates the method using synthetic models contaminated with a wide range of noise, regional self-potential signatures, and interference anomalies. Real data examples are carefully analyzed and discussed in the “Field examples” section, and finally some findings are reported.

Model 4: noisy composite SP anomaly. (a) SP anomaly (generated by the SP data rendered in Fig. 6a) with 20% random noise. (b) Corresponding derivatives. (c) Amplitude. (d) Image.

The Hi’iaka self-potential anomaly, the Kilauea volcano, Hawaii, USA. Profile of the self-potential measurements (dashed line). Inferred location of the Hi’iaka dike (solid line) (taken from Davis83 with permission from Elsevier).

The self-potential signature (V) of some simple geometrical sources at an observation point (\(x_j\), z) along profile (Fig. 1) is given by, e.g., Yungul53 and Mehanee66:

where \(x_j\) (m) is the coordinate of the measurement station (Fig. 1), j is the index of the measurement station, \(x_{\circ }\) (m) is the origin point of the self-potential profile, \(z_\circ\) and z (m) are the coordinates of the buried body and the observation station, \(\theta\) (degrees) is the polarization angle, q (dimensionless) is the shape factor (q = 1.5, 1 and 0.5 for sphere, horizontal cylinder and semi-infinite vertical cylinder), n is the number of data points and K is the electric dipole moment. It is noted that the unit of K (mV m\(^{2q-1}\)) is function of the shape factor (q)62, and that \(\theta\) is measured clockwise (Fig. 1) and ranges from 0 to − 180\(^{\circ }\) in the above formula8.

The AS expression70 is:

where \(\frac{\partial V}{\partial z}\) and \(\frac{\partial V}{\partial x_{j}}\) are the vertical and horizontal derivatives of the self-potential anomaly.

The analytic signal’s amplitude (\(|A_{s}(x_j,z)|\)) of the self-potential anomaly is given by Nabighian70:

Taking the vertical and horizontal derivatives of formula (1), and by substituting the results in expression (3), we get:

The R-parameter (correlation coefficient) is dependent upon both of the amplitude analytic signal of the actual self-potential data (Aso) observed along profile and that of the analytic signal amplitude of the calculated (theoretical) self-potential data (Ast) produced by an assumed source (for example a sphere):

The analytic signal [Aso] can be assessed numerically using expression (3), whereas the analytic signal [Ast] of an assumed source is calculated analytically from expression (4).

To calculate the parameters of a presumed source, the imaging parameter (R-parameter) is first plotted on a 2D map from which the depth \(z_o\) can be readily read. It is worthy noting that when the parameters of the SP profile (Fig. 1) coincide with the buried anomalous source, the imaging parameter achieves its maximum value (so-called here R-max). Further pertinent detail is provided in the “Numerical examples” section. As is seen from expression (5), the computation of R(\(x_o\), \(z_o\), \(\theta\)) does not need knowledge of the electric dipole moment K, which is calculated from the maximum self-potential response. Using expression (1), the predicted self-potential response is then calculated. Fig. 2 presents a flowchart showing the workflow of the developed scheme, which takes about 2 s on a simple PC to estimate the parameters of the interpretive model that resembles the buried anomaly.

The Hi’iaka self-potential anomaly, the Kilauea volcano, Hawaii, USA. (1) The Hi’iaka SP anomaly profile surveyed in 1973. (2) The Hi’iaka SP anomaly profile surveyed in 1995. (3) The Hi’iaka SP anomaly profile surveyed in 1997. (4) The Hi’iaka SP anomaly profile surveyed in 2012. For each anomaly profile, (a) observed and calculated data. (b) Derivative of the observed SP data illustrated in (a). (c) Analytic signal amplitude. (d) Image.

The Osnabr\(\ddot{{u}}\)ck self-potential anomaly, Germany. Location of the survey area (star) north of Osnbrück, and isolines of Vitrinite Reflectance of the maturity map (taken from Gurk et al.6 with permission from Elsevier). NL: Netherlands, B: Belgium.

The approach proposed here has been examined on synthetic self-potential data generated by various source models (for example, horizontal cylinder, sphere and vertical cylinder). The suggested scheme is first verified on numerical experiments without noise. After that, the data have been contaminated with realistic noise, and interpreted in order to assess the stability of the scheme e.g.8,84. Second, in order to assess the stability of the scheme further, the effect of the regional background (that is embedded into the measured self-potential data) on the results is carefully investigated.

The Osnabr\(\ddot{{u}}\)ck self-potential anomaly, Germany. (a) Observed and calculated SP anomaly profile. (b) Derivatives of the observed SP anomaly. (c) Analytic signal amplitude. (d) Image (R-max = 0.9976 at q = 1.2, \(\theta\) = \(-97^{\circ }\), \(z_o\) = 23 m and \(x_o\) = 250 m).

The self-potential response (Fig. 3a) of an idealized body of a horizontal cylinder shape (K = 2500 mV m, \(z_o\) = 7 m, \(\theta\) = \(- 40^{\circ }\), \(x_o\) = 60 m and profile length = 120 m) is calculated from formula (1). Following the recipe discussed above for the interpretation scheme proposed here, Fig. 3b renders the derivatives (horizontal and vertical) of the self-potential anomaly (Fig. 3a). The corresponding analytic signal amplitude (Fig. 3c) is then calculated from the spatial derivatives (Fig. 3b) using expression (3).

The mosaic surface S (which was gridded into 1-m spaces in the x- and z-directions) extended, respectively, to 120 \(\times\) 12 m in these directions (that is (\(x_{o}\), \(z_o\)) \(\in\) S = (0, 120) \(\times\) (1, 12)), and was used to compute and map the correlation coefficient (R-parameter). Expression (5) is employed to calculate the R-parameter for each possible sources (q = 0.5–1.5), where the largest value (R-max) of the R-parameter is attained at the true assumed source (that is q = 1 and R-max = 1.0) (Table 1). Figure 3d presents the image of the R-parameter composed using expression (5) assuming source of a horizontal cylinder model. The R-parameter’s maximum value is marked by the black dot, which denotes the true model parameters of the buried structure (Fig. 3d).

We reiterate that the R map (Fig. 3d) shows the 2D distribution of the obtained R-parameter values. The R-parameter measures the goodness of fit between the observed and predicted self-potential data, and is not representative of geologic structures. An R value of 1 means that the observed and predicted self-potential data are in perfect fit.

To further assess the developed imaging scheme, a number of shape values have been investigated. The scheme is found stable and can retrieve the true values of the model parameters as can be seen from the results presented in Fig. 4 and Table 1.

Osnabr\(\ddot{{u}}\)ck self-potential anomaly, Germany: Two-dimensional electromagnetic inversion results (taken from Gurk et al.6 with permission from Elsevier).

The Suleymankoy self-potential anomaly, Turkey. (a) Observed and calculated SP anomaly profile. (b) Derivatives of the observed SP anomaly. (c) Analytic signal amplitude. (d) Image (R-max = 0.9985 at q = 0.8, \(\theta\) = \(-130^{\circ }\), \(z_o\) = 27 m and \(x_o\) = 72 m).

We added 20% random noise into the self-potential data (Fig. 3a) of the above-mentioned synthetic model (Fig. 5a). Figure 5b,c show the spatial derivatives of the noisy self-potential anomalous signature (Fig. 5a) and their corresponding AS amplitude. Figure 5d shows the maximum R-parameter value (black dot) with a magnitude of 0.98. The estimated model parameters (K= 2906.30 mV m, \(z_o\) = 7.4 m, \(\theta\) = \(- 42.60^{\circ }\), and \(x_o\) = 60 m for an assumed shape factor q of 1.0) (Fig. 5d) are in good agreement with the actual ones.

It can be concluded from the above analysis that the R-parameter imaging method can produce accurate model parameters when the self-potential data are contaminated with noise.

The Malachite Mine self-potential anomaly, USA. (a) Observed and calculated SP anomaly profile. (b) Derivatives of the observed SP anomaly. (c) Analytic signal amplitude. (d) Image (R-max = 0.9841 at q = 0.67, \(\theta\) = \(-112^{\circ }\), \(z_o\) = 15 m and \(x_o\) = 88 m).

A synthetic self-potential anomaly for a sphere-shaped body (K = 1000 \({mV\,m}^2\), \(z_o\) = 4 m, \(\theta\) = \(- 25^{\circ }\), \(x_o\) = 60 m and profile length = 120 m) combined with a first-order regional anomaly generates the composite self-potential anomaly that is shown in Fig. 6a. The simulation formula of the composite anomaly has the form:

Following the procedures discussed above for Model 1, the derivatives of the composite SP response (Fig. 6a) are depicted in Fig. 6b. Figure 6c,d depict the corresponding analytic signal and distribution of the R-parameter values. Using a q of 1.5 (a sphere model), the pertinent model parameters (K=1155.6 mV m\(^2\), \(z_o\) = 4.3 m, \(\theta\) = \(- 25^{\circ }\), \(x_o\) = 61 m) inferred from the imaging scheme are found in good agreement with the true values.

To assess the accuracy of the developed imaging scheme, the composite SP data profile (Fig. 6a) has been contaminated with 20% random noise (Fig. 7a). Figure 7b,c demonstrate the derivatives and the AS amplitude. The correlation map (Fig. 7d) rendered a maximum value of 0.74, which corresponds to an inverse model of K= 1640.60 mV m\(^2\), \(z_o\) = 5.4 m, \(\theta\) = \(-28^{\circ }\) and \(x_o\) = 61 m for q = 1.5, which matches well the actual source model. This supports that the method has potential in exploration geophysics.

The scheme is analyzed on five published real self-potential data from geothermal systems and mineral exploration in the following sections.

A number of self-potential geophysical surveys were carried out in 1973, 1995, 1997 and 2012 over a basaltic dike intrusion (referred to as the Hi’iaka dike, Hawaii) (Fig. 8). The Hi’iaka dike intruded into the upper part of the Kilauea volcano, which is associated with the Hi’iaka and Pauahi craters eruption alongside the Kilauea rift zone87,88. It made a 100-m long surface fracture that erupted magma southwest of the Hi’iaka crater. Measurements of geophysical surveying recommended that the fracture is continued about 1.5 km under the subsurface in the southwest direction (Fig. 8).

The measurements of the self-potential data over the Hi’iaka dike intrusion commenced in 1973 by Zablocki89 and Zablocki90, and continued in 1995, 1997 and 201283. Localized fluid disruption is blamed for the SP anomaly83,89,90. Davis83 stated that “geothermal reservoirs are found above intrusions of magma such as dikes or dike swarms, which set up hydrothermal circulation generating hot water and steam from which energy can be tapped”. The self-potential anomaly profiles of 1973, 1995, 1997, and 2012 are digitized into 10-m intervals (Fig. 9(1)–(4)). The Hi’iaka SP anomaly profiles have been interpreted by Davis83 using the self-potential inversion approach of Sill38. Davis83 interpreted the profiles by a trapezoidal source (approximated by a dike shaped model) located at different depths ranging from 50 to 190 m, and attributed the increase in depth to the magma cooling and heat loss at the top of the dike.

The aforementioned Hi’iaka self-potential anomaly profiles of 1973, 1995, 1997 and 2012 are interpreted using the scheme developed here (Fig. 9(1)–(4)). The derivatives of the self-potential response of each profile (Fig. 9(1)a,(2)a,(3)a,(4)a), and the corresponding analytic signal amplitude are rendered in Fig. 9(1)b,(2)b,(3)b,(4)b, and in Fig. 9(1)c,(2)c,(3)c,(4)c. The R-parameter values are reported in Fig. 9(1)d,(2)d,(3)d,(4)d with the R-max value for each SP profile. For the 1973, 1995, 1997, 2012 profiles, the estimated model parameters are found (K = -10688 mV m\(^{{2q-1}}\), \(z_o\) = 54 m, \(\theta\)= − 90, \(x_o\) = 300 m, and q = 0.8 with an R-max of 0.99), (K = − 10072 mV m\(^{{2q-1}}\), \(z_o\) = 57 m, \(\theta\)= − 90, \(x_o\) = 280 m, and q = 0.8 with an R-max of 0.98), (K = − 4340 mV m\(^{{2q-1}}\), \(z_o\) = 69 m, \(\theta\)= − 110, \(x_o\)= 310 m, and q = 0.7, with an R-max of 0.98), and (K = − 1718 mV m\(^{{2q-1}}\), \(z_o\) = 177 m, \(\theta\)= -115, \(x_o\) = 340 m, and q = 0.6 with an R-max of 0.98). Table 2 tabulates the recovered model parameters for each profile, and shows that the observed self-potential anomaly is fit by a dike-like model with a shape factor of 0.6–0.8 (that is q = 0.6–0.8). Analysis shows that there is good match between the depths of the interpretive source (trapezoidal model) stated in the published literature and the depths obtained here (Table 3). The match between the observed and calculated self-potential data for each profile is depicted in Fig.9(1)a–(4)a, which is quite good.

It is re-noted that Davis83 reported that the SP anomaly remained strong throughout the measurement duration. However, the SP anomaly of 2012 is about 60% of that of 1973. Therefore, the variation in the recovered depths (54–177 m, Table 3) of the interpretated self-potential profiles (measured in 1973–2012) is not unexpected, and is attribuated to the magma cooling and heat loss at the top of the dike83.

The Bavarian woods self-potential anomaly, Germany. (a) Observed and calculated SP anomaly profile. (b) Derivatives of the observed SP anomaly. (c) Analytic signal amplitude. (d) Image (R-max = 0.99 at q = 0.94, \(\theta\) = \(- 62^{\circ }\), \(z_o\) = 41 m and \(x_o\) = 263 m).

A self-potential anomaly near the Osnabr\(\ddot{{u}}\)ck area (Fig. 10), Northwest Germany6 has been carried out to trace a graphite anomaly that has a quasi-vertical form in the Lias-epsilon shales. Gurk et al.6 found a significant single self-potential anomaly of roughly − 600 mV (Fig. 11a), which supports conductive graphite minerals. The 500-m long self-potential anomaly profile is meshed into 5-m intervals (Fig. 11a).

Figure 11b and c present the derivatives, and the AS amplitude of the SP anomaly, respectively. The R-parameter maximum value (R-max = 0.9976, black dot, Fig. 11d) was determined with the corresponding best interpretive parameters: K = 46527 mV m\(^{{2q-1}}\), \(z_o\) = 23 m, \(\theta\) = \(-97^{\circ }\), \(x_o\) = 250 m, and q = 1.2 (Fig. 11a,d and Table 4). According to the recovered R-parameters, the subsurface anomalous body is approximated by a horizontal cylinder like-structure with a horizontal location of 250 m and an estimated depth to the center of 23 m, which is in good agreement with the interpreted results of Gurk et al.6 and Mehanee8 (Table 5). The variation in the magnitude of the model parameter K is due to the inconsistent use of unit (Table 5) as the interpretive models are not quite identical; they range from thin sheet to quasi-horizontal cylinder.

In order to map the 2D electric conductivity (inverse of resistivity) distribution in the underground, Gurk et al.6 measured a radio magnetotelluric data (apparent resistivities and phase) profile on the initial 400 m of the self-potential profile described earlier. The corresponding 2D inverse results of Gurk et al.6 revealed a prominent conductive anomalous body (Fig. 12), the location and depth of which correlate well with the results inferred from the approach developed here (Fig. 11).

The SP anomaly of Suleymankoy53 was carried out for copper deposits. The mine is characterized by alpine ophiolite containing several copper deposits. The anomaly is gridded at intervals of 2 m long (Fig. 13a). The self-potential anomaly of Suleymankoy is interpreted using the presented R-parameter imaging technique. Figure 13b–d depict the corresponding derivatives, the AS amplitude, and the imaging, which reveals an R-max value of 0.9985. The parameters revealed from interpretation are K = 1898 mV m\(^{{2q-1}}\), \(z_o\) = 27 m, \(\theta\) = \(-130^{\circ }\), \(x_o\) = 76 m and q = 0.8 (Fig. 13a–d and Table 6). The observed and calculated self-potential data are rendered in Fig. 13a. Table 7 presents a comparison between the obtained results and those mentioned in the published literature.

As can be seen from Table 7, the reported depths are in reasonable agreement, whereas the electric dipole moments (K) encountered some variations, which could be attributed to the various approximations employed in the interpretation schemes used in this table and how the parameter K is calculated from these schemes.

The Malachite Mine is an amphibolite belt bounded by schists and gneisses. The self-potential profile over the Malachite Mine is digitized at intervals equal to 2 m (Fig. 14a).

The Malachite SP anomaly is interpreted using the R-parameter imaging technique. Figure 14b,c render the gradients, and the AS amplitude. Upon applying the imaging technique, an R-max of 0.9841 (Fig. 14d) was obtained with the following interpretive parameters: K = 515 mV, \(z_o\) = 15 m, \(\theta\) = \(-112^{\circ }\), \(x_o\) = 88 m, and q = 0.67 (Fig. 14a,d and Table 8). The subsurface structure was approximated by a semi-infinite vertical cylindrical structure with a horizontal spatial location of 88 m. The estimated depth (15 m to the top of the structure) is in good match with the drilling information and previous interpreted works (Table 9).

Table 9 shows that the parameter K encountered some variation; this is attributed to two main reasons. First, the inconsistent use of unit as the interpretive models are not quite identical; they range from vertical cylinder to quasi-vertical cylinder. Second, the nature of the approximations employed in interpretation schemes reported in this table, and how the parameter K is calculated from these schemes.

Figure 15a depicts the self-potential anomaly collected over a graphite ore body at the southern Bavarian woods, Germany95. The self-potential anomaly profile is digitized using a 1-m sampling interval. Several authors have interpreted this anomaly profile. Al-Garani96 interpreted the anomaly as a quasi-vertical cylinder using the neural network inversion with \(z_o\) = 33 m (depth to the top). Mehanee66 analyzed this SP profile with a horizontal cylinder of a depth to the center \(z_o\) of 46 m using a regularized inversion. Gokturkler and Balkaya5 described the anomaly by a horizontal cylinder model using a genetic algorithm (\(z_o\) = 45.03 m; to the center), simulated annealing (\(z_o\) = 47.59 m; to the center) and particle swarm optimization algorithm (\(z_o\) = 47.59 m). Di Maio et al.97 fit the profile by a horizontal cylinder by applying a spectral analysis and tomographic approach with \(z_o\) = 44.9 m (to the center).

We interpreted this self-potential anomaly profile using the R-parameter imaging technique. The analytic signal amplitude anomaly calculated from the horizontal and vertical derivatives (Fig. 15b) is presented in Figure 15c. The R-parameter values are rendered in Fig. 15d. An R-max of 0.9983 pertains to an interpretive model of K = 13821 mV m\(^{{2q-1}}\), \(z_o\) = 41 m, \(\theta\) = \(-62^{\circ }\), \(x_o\) = 263 m, and q = 0.94 (Fig. 15a,d and Table 10). The presented analysis shows that the depth and shape (which resembles a quasi-horizontal cylinder model) are generally in good agreement with the aforementioned results but the results of Al-Garani96, who interpreted the data by a quasi-vertical cylinder (q = 0.7) (Table 11). We do not expect all interpretation methods to yield the same results as each method has its own assumptions and limitation.

As mentioned above, the R map shows the 2D distribution of the obtained R-parameter values. The R-parameter measures the goodness of fit between the observed self-potential data, and the theoretical self-potential data generated from the model parameters (z, K, \(\theta\)) of the interpretive idealized model. The R parameter value does not provide uncertainty estimate for the evolved model parameters.

The non-uniqueness is one of the most challenging issues in geophysical data interpretation e.g.99, where multiple approximative solutions can equally fit the observed data. Joint inversion could help minimize this issue and provide better understanding e.g.84. It is worthy noting that it is very rare to solely use/measure one data kind when it comes to a detailed geophysical prospecting program. In industry, multiple data sets are essential for comprehensive understanding and for maximizing the potential of the underlying exploration program. Usually multiple geophysical data along with geological information are used, inverted and interpreted in an integrated manner (the so-called joint interpretation) to hopefully recover and select a unique inverse model the data of which match the measured data sets, and that fits into the underlying geologic setting of the area under study. May be this is the best we can do in exploration geophysics in order to resolve the non-uniqueness issue of an inverse problem solution.

A rapid imaging scheme has been developed for the interpretation of self-potential data. In about 2 s on a simple PC, the scheme can estimate the parameters of the interpretive model (which is in the context of sphere, horizontal cylinder or vertical cylinder) that resembles the buried structure. The developed scheme uses the amplitude of the AS of the self-potential data undergoing interpretation and the amplitude of the AS of the self-potential data calculated by the assumed interpretive model to construct the corresponding 2D image of the so-called R-parameter. The scheme attains the largest value of the R-parameter when the recovered parameters coincide with the actual ones. It is noted that the R-parameter is independent of the electric dipole moment (K). The analyzed numerical examples demonstrated the stability of the developed scheme, and that its accuracy can be affected by the nearby geological structures. The five field data examples (from geothermal systems and mineral prospecting) analyzed here show that the scheme is capable of producing good results that agree well with those reported in other published research. The developed imaging scheme can have some potential in geothermal investigation and reconnaissance studies.

The self-potential computer code and the datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Fox, R. W. On the electro-magnetic properties of metalliferous veins in the mines of Cornwall. Philos. Trans. R. Soc. Lond. 120, 399–414 (1830).

ADS Google Scholar

Abbas, M. et al. Redox potential distribution of an organic-rich contaminated site obtained by the inversion of self-potential data. J. Hydrol. 554, 111–127 (2017).

Article ADS CAS Google Scholar

Barde-Cabusson, S., Finizola, A. & Grobbe, N. A practical approach for self-potential data acquisition, processing, and visualization. Interpretation 9, T123–T143. https://doi.org/10.1190/INT-2020-0012.1 (2021).

Article Google Scholar

Fedi, M. & Abbas, M. A. A fast interpretation of self-potential data using the depth from extreme points method. Geophysics 78, E107–E116 (2013).

Article ADS Google Scholar

Göktükler, G. & Balkaya, C. Inversion of self-potential anomalies caused by simple geometry bodies using global optimization algorithms. J. Geophys. Eng. 9, 498–507 (2012).

Article Google Scholar

Gurk, M., Bosch, F. P. & Tougiannidis, N. Electric field variations measured continuously in free air over a conductive thin zone in the tilted Lias-epsilon black shales near Osnabr\(\ddot{\text{ u }}\)ck, Northwest Germany. J. Appl. Geophys. 91, 21–30 (2013).

Article ADS Google Scholar

Ikard, S. J. et al. Saline pulse test monitoring with the self-potential method to nonintrusively determine the velocity of the pore water in leaking areas of earth dams and embankments. Water Resour. Res. 48, W04201 (2012).

Article ADS Google Scholar

Mehanee, S. A. Simultaneous joint inversion of gravity and self-potential data measured along profile: Theory, numerical examples, and a case study from mineral exploration with cross validation from electromagnetic data. IEEE Trans. Geosci. Remote Sens. 60, 1–20. https://doi.org/10.1109/TGRS.2021.3071973 (2022).

Article Google Scholar

Patella, D. Introduction to ground surface self-potential tomography. Geophys. Prospect. 45, 653–681 (1997).

Article ADS Google Scholar

Patella, D. Self-potential global tomography including topographic effects. Geophys. Prospect. 45, 843–863 (1997).

Article ADS Google Scholar

Rao, K., Jain, S. & Biswas, A. Global optimization for delineation of self-potential anomaly of a 2D inclined plate. Nat. Resour. Res. 30, 175–89 (2021).

Article Google Scholar

Revil, A. & Jardani, A. The self-potential method theory and applications in environmental geosciences (Cambridge University Press, 2013).

Book Google Scholar

Revil, A., Ehouarne, L. & Thyreault, E. Tomography of self-potential anomalies of electrochemical nature. Geophys. Res. Lett. 28, 4363–4366 (2001).

Article ADS CAS Google Scholar

Sato, M. & Mooney, H. M. The electrochemical mechanism of sulfide self-potentials. Geophysics 15, 226–249 (1960).

Article Google Scholar

Zhu, Z. et al. Autonomous underwater vehicle based marine multicomponent self-potential method: Observation scheme and navigational correction. Geosci. Instrum. Methods Data Syst. 10, 1–25 (2020).

ADS Google Scholar

Biswas, A., Mandal, A., Sharma, S. P. & Mohanty, W. K. Delineation of subsurface structure using self-potential, gravity and resistivity surveys from South Purulia Shear Zone, India: Implication to uranium mineralization. Interpretation 2, T103–T110 (2014).

Article Google Scholar

Essa, K., Mahanee, S. & Smith, P. D. A new inversion algorithm for estimating the best fitting parameters of some geometrically simple body to measured self-potential anomalies. Explor. Geophys. 39, 155–163 (2008).

Article ADS Google Scholar

Essa, K. S. A particle swarm optimization method for interpreting self potential anomalies. J. Geophys. Eng. 16, 463–477 (2019).

Article Google Scholar

Kawada, Y. & Kasaya, T. Marine self-potential survey for exploring sea floor hydrothermal ore deposits. Sci. Rep. 7, 13552. https://doi.org/10.1038/s41598-017-13920-0 (2017).

Article ADS CAS PubMed PubMed Central Google Scholar

Kawada, Y. & Kasaya, T. Self-potential mapping using an autonomous underwater vehicle for the Sunrise deposit, Izu-Ogasawara arc, southern Japan. Earth Planets Space 70, 142–156 (2018).

Article ADS Google Scholar

Rani, K., Guha, A., Pal, S. K. & Kumar, K. V. Broadband reflectance, emittance spectroscopy and self potential geophysical survey for targeting gold sulphide lode deposit in Bhukia, Rajasthan, India. Geocarto Int.https://doi.org/10.1080/10106049.2018.1508310 (2018).

Article Google Scholar

Sultan, A. S., Salah, A. M., Fernando, M. S. & Ahmed, S. H. Geophysical exploration for gold and associated minerals, case study: Wadi El Beida area, south eastern desert. Egypt. J. Geophys. Eng. 6, 345–56 (2009).

Article Google Scholar

Biswas, A. A review on modeling, inversion and interpretation of self-potential in mineral exploration and tracing paleo-shear zones. Ore Geol. Rev. 91, 21–56 (2017).

Article Google Scholar

Mehanee, S. Tracing of paleo-shear zones by self-potential data inversion: Case studies from the KTB, Rittsteig, and Grossensees graphite-bearing fault planes. Earth Planets Space 67, 14 (2015).

Article ADS Google Scholar

Wynn, J. C. & Sherwood, S. I. The self-potential (SP) method: An inexpensive reconnaissance archaeological mapping tool. J. Field Archaeol. 11, 195–204 (1984).

Google Scholar

Hunter, L., & Powers, M. Geophysical investigations of earthen dams: An overview, 21th SAGEEP, pp. 1083–1096 (2008).

Vichabian, Y. & Morgan, F. D. Self-potentials in cave detection. Lead. Edge 21, 866–871. https://doi.org/10.1190/1.1508953 (2002).

Article Google Scholar

Karaoulis, M., Revil, A. & Mao, D. Localization of a coal seam fire using combined self-potential and resistivity data. Int. J. Coal Geol. 128–129, 109–118 (2014).

Article Google Scholar

Shao, Z., Wang, D., Wang, Y. & Zhong, X. Theory and application of magnetic and self-potential methods in the detection of the Heshituoluogai coal fire, China. J. Appl. Geophys. 104, 64–74 (2014).

Article ADS Google Scholar

Shao, Z. et al. Experimental study of the self-potential anomaly caused by coal fires. J. Appl. Geophys. 145, 124–132 (2017).

Article ADS Google Scholar

DesRoches, A. J., Butler, K. E. & MacQuarrie, K. T. B. Surface self-potential patterns related to transmissive fracture trends during a water injection test. Geophys. J. Int. 212, 2047–2060 (2018).

Article ADS CAS Google Scholar

Roubinet, D., Linde, N., Jougnot, D. & Irving, J. Streaming potential modeling in fractured rock: Insights into the identification of hydraulically active fractures. Geophys. Res. Lett. 43, 4937–4944 (2016).

Article ADS Google Scholar

Soueid Ahmed, A. S., Jardani, A., Revil, A. & Dupont, J. P. Specific storage and hydraulic conductivity tomography through the joint inversion of hydraulic heads and self-potential data. Adv. Water Resour. 89, 80–90 (2016).

Article ADS Google Scholar

Heinze, T., Limbrock, J. K., Pudasaini, S. P. & Kemna, A. Relating mass movement with electrical self-potential signals. Geophys. J. Int. 216, 55–60 (2019).

Article ADS Google Scholar

Eppelbaum, L. V. Review of processing and interpretation of self-potential anomalies: Transfer of methodologies developed in magnetic prospecting. Geosciences 11, 194. https://doi.org/10.3390/geosciences11050194 (2021).

Article ADS CAS Google Scholar

Ishido, T. Electrokinetic mechanism for theW-shaped self-potential profile on volcanoes. Geophys. Res. Lett. 31, L15616. https://doi.org/10.1029/2004GL020409 (2004).

Article ADS Google Scholar

Kilty, K. T. On the origin and interpretation of self-potential anomalies. Geophys. Prospect. 32, 51–62 (1984).

Article ADS Google Scholar

Sill, W. R. Self-potential modeling from primary flows. Geophysics 48, 76–86 (1983).

Article ADS Google Scholar

Stoll, J., Bigalke, J. & Grabner, E. W. Electrochemical modelling of self-potential anomalies. Surv. Geophys. 16, 107–120 (1995).

Article ADS Google Scholar

Fournier, C. Spontaneous potentials and resistivity surveys applied to hydrogeology in a volcanic area: Case history of the Chaine des Puys (Puy-de-Dome, France). Geophys. Prospect. 37, 647–668 (1989).

Article ADS Google Scholar

Abdelrahman, E. M., Saber, H. S., Essa, K. S. & Fouda, M. A. A least-squares approach to depth determination from numerical horizontal self-potential gradients. Pure Appl. Geophys. 161, 399–411 (2004).

Article ADS Google Scholar

Abdelrahman, E. M., Essa, K. S., El-Araby, T. M. & Abo-Ezz, E. R. A least-squares depth-horizontal position curves method to interpret residual SP anomaly profile. J. Geophys. Eng. 3, 252–259 (2006).

Article Google Scholar

Abdelrahman, E. M. et al. New least-square algorithm for model parameters estimation using self-potential anomalies. Comput. Geosci. 34, 1569–1576 (2008).

Article ADS Google Scholar

Babu, R. H. V. & Rao, A. D. A rapid graphical method for the interpretation of the self-potential anomaly over a two-dimensional inclined sheet of finite depth extent. Geophysics 53, 1126–1128 (1988).

Article ADS Google Scholar

Banerjee, B. Quantitative interpretation of self-potential anomalies of some specific geometric bodies. Pure Appl. Geophys. 90, 138–152 (1971).

Article ADS Google Scholar

Bhattacharya, B. B. & Roy, N. A note on the use of nomograms for self-potential anomalies. Geophys. Prospect. 29, 102–107 (1981).

Article ADS Google Scholar

Essa, K. S. Self potential data interpretation utilizing the particle swarm method for the finite 2D inclined dike: Mineralized zones delineation. Acta Geod. Geophys. 55, 203–221 (2020).

Article Google Scholar

Essa, K. S., Diab, Z. E. & Mehanee, S. Self-potential data inversion utilizing the Bat optimizing algorithm (BOA) with various application cases. Acta Geophys. 71, 567–586 (2023).

Article ADS Google Scholar

Fitterman, D. V. Calculations of self-potential anomalies near vertical contacts. Geophysics 44, 195–205 (1979).

Article ADS Google Scholar

Mehanee, S., Essa, K. S. & Smith, P. A rapid technique for estimating the depth and width of a two-dimensional plate from self-potential data. J. Geophys. Eng. 8, 447–456 (2011).

Article Google Scholar

Rao, A. D. & Babu, R. H. V. Quantitative interpretation of self potential anomalies due to two-dimensional sheet-like bodies. Geophysics 48, 1659–1664 (1983).

Article ADS Google Scholar

Soueid Ahmed, A., Jardani, A., Revil, A. & Dupont, J. P. SP2DINV: A 2D forward and inverse code for streaming potential problems. Comput. Geosci. 59, 9–16 (2013).

Article ADS Google Scholar

Yüngül, S. Interpretation of spontaneous polarization anomalies caused by spherical ore bodies. Geophysics 15, 237–246 (1950).

Article ADS Google Scholar

Colangelo, G., Lapenna, V., Perrone, A., Piscitelli, S. & Telesca, L. 2D self-potential tomographies for studying groundwater flows in the Varco d’Izzo landslide (Basilicata, southern Italy). Eng. Geol. 88, 274–286. https://doi.org/10.1016/j.enggeo.2006.09.014 (2006).

Article Google Scholar

Hämmann, M., Maurer, H. R., Green, A. G. & Horstmeyer, H. Self-potential image reconstruction: Capabilities and limitations. J. Environ. Eng. Geophys. 2, 21–35 (1997).

Article Google Scholar

Hermans, T., Nguyen, F., Robert, T. & Revil, A. Geophysical methods for monitoring temperature changes in shallow low enthalpy geothermal systems. Energies 7, 5083–5118. https://doi.org/10.3390/en7085083 (2014).

Article Google Scholar

Mao, D. et al. Resistivity and self-potential tomography applied to groundwater remediation and contaminant plumes: Sandbox and field experiments. J. Hydrol. 530, 1–14 (2015).

Article ADS CAS Google Scholar

Minsley, B. J., Sogade, J. & Morgan, F. D. Three-dimensional self-potential inversion for subsurface DNAPL contaminant detection at the Savannah River Site, South Carolina. Water Resour. Res. 43, W04429. https://doi.org/10.1029/2005WR003996 (2007).

Article ADS Google Scholar

Oliveti, I. & Cardarelli, E. 2D approach for modeling self-potential anomalies: Application to synthetic and real data. Boll. Geof. Teor. Appl. 58, 415–430 (2017).

Google Scholar

Rittgers, J. B. et al. Self-potential signals generated by the corrosion of buried metallic objects with application to contaminant plumes. Geophysics 78, 65–82 (2013).

Article Google Scholar

Tikhonov, A. N. & Arsenin, V. Y. Solutions of ill-posed problems (Wiley, 1977).

MATH Google Scholar

Abdelrahman, E. M., Essa, K. S., Abo-Ezz, E. R. & Soliman, K. S. Self-potential data interpretation using standard deviations of depths computed from moving-average residual anomalies. Geophys. Prospect. 54, 409–423 (2006).

Article ADS Google Scholar

Abdelrahman, E. M., Soliman, K., Essa, K. S., Abo-Ezz, E. R. & El-Araby, T. M. A least squares minimization approach to depth determination from numerical second horizontal self-potential anomalies. Explor. Geophys. 40, 214–221 (2009).

Article ADS Google Scholar

Essa, K., & Mehanee, S. A rapid algorithm for self-potential data inversion with application to mineral exploration. In Presented at the 19th International Geophysical Conference and Exhibition, Australian Society of Exploration Geophysicists, 18–22 November, Perth, Australia (2007).

Murty, B. V. S. & Haricharan, P. Nomogram for the spontaneous potential profile over sheet-like and cylindrical two-dimensional sources. Geophysics 50, 1127–1135. https://doi.org/10.1190/1.1441986 (1985).

Article ADS Google Scholar

Mehanee, S. An efficient regularized inversion approach for self-potential data interpretation of ore exploration using a mix of logarithmic and non-logarithmic model parameters. Ore Geol. Rev. 57, 87–115 (2014).

Article Google Scholar

Santos, F. A. M. Inversion of self-potential of idealized bodies’ anomalies using particle swarm optimization. Comput. Geosci. 36, 1185–1190 (2010).

Article ADS Google Scholar

Di Maio, R., Piegari, E. & Rani, P. Source depth estimation of self-potential anomalies by spectral methods. J. Appl. Geophys. 136, 315–325 (2017).

Article Google Scholar

Sungkono Warnana, D. D. Black hole algorithm for determining model parameter in self-potential data. J. Appl. Geophys. 148, 189–200 (2018).

Article ADS Google Scholar

Nabighian, M. N. The analytic signal of two-dimensional magnetic bodies with polygonal cross-section: Its properties and use for automated anomaly interpretation. Geophysics 37, 507–517 (1972).

Article ADS Google Scholar

Essa, K. S., Mehanee, S., Soliman, K. & Diab, Z. E. Gravity profile interpretation using the R-parameter imaging technique with application to ore exploration. Ore Geol. Rev. 126, 103695 (2020).

Article Google Scholar

Mehanee, S., Essa, K. S. & Diab, Z. E. Magnetic data interpretation using a new R-parameter imaging method with application to mineral exploration. Nat. Resour. Res. 30, 77–95 (2021).

Article Google Scholar

Nabighian, M. N. Additional comments on the analytic signal of two-dimensional magnetic bodies with polygonal cross section. Geophysics 39, 507–517 (1974).

Article Google Scholar

Blakely, R. J. Potential theory in gravity and magnetic applications (Cambridge University Press, 1996).

Google Scholar

Everett, M. E. Near-surface applied geophysics (Cambridge University Press, 2013).

Book Google Scholar

Fedi, M. DEXP: A fast method to determine the depth to the sources of potential fields. Geophysicshttps://doi.org/10.1190/1.2144416 (2005).

Article Google Scholar

Fedi, M., Florio, G. & Quarta, T. A. Multiridge analysis of potential fields: Geometric method and reduced Euler deconvolution. Geophysics 74, L53–L65 (2009).

Article ADS Google Scholar

Macleod, I., Jones, K. & Dai, T. 3-D analytic signal in the interpretation of total magnetic field data at low magnetic latitudes. Explor. Geophys.https://doi.org/10.1071/EG993679 (1993).

Article Google Scholar

Milano, M., Fedi, M. & Fairhead, J. D. Joint analysis of the magnetic field and total gradient intensity in central Europe. Solid Earth 10, 697–712 (2019).

Article ADS Google Scholar

Mushayandebvu, M., van Drielz, P., Reid, A. & Fairhead, J. Magnetic source parameters of two-dimensional structures using extended Euler deconvolution. Geophysics 66(3), 814–823 (2001).

Article ADS Google Scholar

Paoletti, V. et al. Insights into the structure and surface geology of Isla Socorro, Mexico, from airborne magnetic and gamma-ray surveys. Surv. Geophys. 37, 601–623 (2016).

Article ADS Google Scholar

Salem, A., Ravat, D., Gamey, T. J. & Ushijima, K. Analytic signal approach and its applicability in environmental magnetic investigations. J. Appl. Geophys. 49, 231–244 (2002).

Article ADS Google Scholar

Davis, P. M. Geothermal evolution of an intruded dike in the rift zone of Kilauea volcano, Hawaii from VLF and self-potential measurements. J. Volcanol. Geotherm. Res. 302, 64–80 (2015).

Article ADS CAS Google Scholar

Sarkar, K., Mukesh, M. & Singh, U. K. Nature-inspired stochastic hybrid technique for joint and individual inversion of DC and MT data. Sci. Rep. 13, 2668. https://doi.org/10.1038/s41598-023-29040-x (2023).

Article ADS CAS PubMed PubMed Central Google Scholar

Srivastava, S. & Agarwal, B. N. P. Interpretation of self-potential anomalies by enhanced local wave number technique. J. Appl. Geophys. 68, 259–268 (2009).

Article ADS Google Scholar

Agarwal, B. & Srivastava, S. Analyses of self-potential anomalies by conventional and extended Euler deconvolution techniques. Comput. Geosci. 35, 2231–2238 (2009).

Article ADS Google Scholar

Klein, F. W., Koyanagi, R. Y., Nakata, J. S. & Tanigawa, W. R. The seismicity of Kilauea’s magma system. In Volcanism in Hawaii (eds Decker, R. W. et al.) 1019–1185 (US Geological Survey Professional Paper, 1987).

Google Scholar

Tilling, R. I. et al. The 1972–1974 Mauna Ulu eruption, Kilauea Volcano: An example of quasi-steady-state magma transfer. In Volcanism in Hawaii (eds Decker, R. W. et al.) 405–469 (U.S. Geological Survey Professional Papers, USA, 1987).

Google Scholar

Zablocki, C.J. (1976). Mapping thermal anomalies on an active volcano by the self-potential method, Kilauea, Hawaii. In Proceeding 2nd U.N. Symposium on the Development and Use of Geothermal Resources, San Francisco, California, May 1975, pp. 1299-1309.

Zablocki, C. J. Applications of the VLF induction method for studying some volcanic processes of Kilauea volcano, Hawaii. J. Volcanol. Geotherm. Res. 3, 155–195 (1978).

Article ADS Google Scholar

Huff, L.C. Comparison of geological, geophysical, and geochemical prospecting methods at the Malachite Mine, Jefferson County, Colorado. U.S. Geol. Surv. Bull. 161–179 (1963).

Dobrin, M. B. Introduction to geophysical prospecting (Mc Graw Hill, 1976).

Google Scholar

Tlas, M. & Asfahani, J. A best-estimate approach for determining self-potential parameters related to simple geometric shaped structures. Pure Appl. Geophys. 164, 2313–2328 (2007).

Article ADS MATH Google Scholar

Abdelrahman, E. M., El-Araby, T. M. & Essa, K. S. Shape and depth determinations from second moving average residual self-potential anomalies. J. Geophys. Eng. 6, 43–52 (2009).

Article Google Scholar

Meiser, P. A method of quantitative interpretation of self-potential measurements. Geophys. Prospect. 10, 203–218 (1962).

Article ADS Google Scholar

Al-Garani, M. A. Interpretation of spontaneous potential anomalies from some simple geometrically shaped bodies using neural network inversion. Acta Geophys. 58, 143–162 (2010).

Article ADS Google Scholar

Di Maio, R., Piegari, E., Rani, P. & Avella, A. Self-Potential data inversion through the integration of spectral analysis and tomographic approaches. Geophys. J. Int. 206, 1204–1220 (2016).

Article ADS Google Scholar

Essa, K. S. A new algorithm for gravity or self-potential data interpretation. J. Geophys. Eng. 8, 434–446 (2011).

Article Google Scholar

Milano, M., Varfinezhad, R., and Fedi, M. Joint inversion of DC resistivity and potential field data under different model weighting functions , EGU General Assembly 2022, Vienna, Austria, 23-27 May, EGU22-4190, (2022). https://doi.org/10.5194/egusphere-egu22-4190.

Download references

Department of Geophysics, Faculty of Science, Cairo University, Giza, 12613, Egypt

Salah A. Mehanee, Khalid S. Essa, Khaled S. Soliman & Zein E. Diab

School of Earth Sciences, The University of Adeliede, Adelaide, SA, Australia

Salah A. Mehanee

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

S.M. and K.E. proposed the idea and developed the algorithm. K.S. analyzed the results with S.M. and K.E., and assisted in the interpretation of real data. S.M. wrote the paper, and K.E. and K.S. reviewed the manuscript and provided feedback. Z.D. prepared some of the figures, developed the code with S.M., and assisted in the interpretation of real data. All authors read and approved the final paper.

Correspondence to Salah A. Mehanee.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Mehanee, S.A., Essa, K.S., Soliman, K.S. et al. A fast imaging method for the interpretation of self-potential data with application to geothermal systems and mineral investigation. Sci Rep 13, 13548 (2023). https://doi.org/10.1038/s41598-023-39672-8

Download citation

Received: 16 November 2022

Accepted: 28 July 2023

Published: 20 August 2023

DOI: https://doi.org/10.1038/s41598-023-39672-8

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.