Numerical results using the least-squares approach for a gravity model due to horizontal cylinder without and with a 10% of random noise (
A combined gravity and/or self-potential anomaly formula is utilized to estimate the model parameters of the buried geologic structures represented by simple geometric. The simple geometric shapes (spheres, cylinders, and sheets) are not really found but often applied to reduce the nonuniqueness in interpreting the gravity and self-potential data. Numerous approaches through the combined formula such as least squares, Werner deconvolution, and the particle swarm optimization method are used. The application of these methods was demonstrated by applying a synthetic gravity and self-potential example without and with 10% random noise to compare their efficiency in estimating the model parameters of the buried structures. Besides, they were applied to two field data for mineral exploration. The appraised model parameter values from each method were compared together and with those published in literature.
- gravity and/or self-potential
- model parameters
- mineral exploration
Minerals exploration is vital in many countries to increase the income of their people and their economy relies upon discovering minerals. The minerals or ores mined have different variety according to its important in the economy. Geophysical passive method such as gravity and self-potential play an important role in discovering these minerals or ores [1, 2, 3, 4, 5]. The gravity method based on measuring the variations in the Earth’s gravitational field resulting from the density differences between the subsurface rocks while the self-potential method depended on the electrical potential that develops on the earth’s surface due to flow of the natural electrical current on the subsurface [6, 7]. The interpretation of gravity and self-potential data falls on the main two categories as follows: the first category depends on three-dimensional and two-dimensional data elucidation [8, 9, 10, 11, 12, 13], the second category is depending using the simple geometric-shaped model such as spheres, cylinders, and sheets which are playing a vital role in interpreting the subsurface structures to reach the priors information that help in more investigations [14, 15, 16, 17, 18, 19, 20]. In addition, methods depend on the global optimization algorithms such as genetic algorithm [21, 22, 23, 24], particle swarm [25, 26], simulated annealing [27, 28, 29, 30, 31, 32], flower pollination , memory-based hybrid dragonfly , differential evolution [35, 36].
Here, a combined formula for both gravity and self-potential  is applied to construct this chapter. Moreover, this formula is used to calculate the buried model parameters, for example in case of self-potential data, the parameters are the electric dipole moment or the amplitude coefficient (
2. The suggested combined gravity or self-potential formula
Secondly, the self-potential anomaly formula for the same simple geometric models is 
In Refs. [1, 37], Eqs. (1) and (2) were used to join together to produce a combined gravity or self-potential formula for the simple geometric structures such as a semi-infinite vertical cylinder, a dike, a horizontal cylinder, and a sphere (Figure 1) as follows:
2.1 The least-squares approach
Essa  developed this approach, which was relied on solving the problem of finding the depth from the measured data by solving a nonlinear form
2.2 Werner deconvolution approach
Werner deconvolution was proposed by Werner in 1953 . This approach is used to estimate mainly the origin location and the depth of the buried structures. Werner proposed to transform the equation of unknown parameters into a rational function. Eq. (3) can be rewritten in linear form follow:
Eq. (5) is linear form in the four variables
2.3 The particle swarm approach
The particle swarm was suggested by  and has many various applications, for example, in geophysics [40, 41, 42]. For more detail in this approach, you find it many published literature [43, 44]. The model parameters values of the unknowns are relied upon the objective function, so that every problem can be resolved. In this approach, the particles represent the parameter which we are invert. In the beginning, each particle has a location and velocity. After that each particle changes its location (
where is the velocity of the particle
The five source parameters (
3. Synthetic example
To test the ability of each suggested approach in assessing the buried model parameters for the simple geometric shapes such as spheres, cylinders, and sheets. Two synthetic examples are suggested for these interpretation. First one is belonging to use the gravity data and second is applying the self-potential data.
3.1 Gravity anomaly model
A gravity anomaly of a horizontal cylinder model is generated using the following parameters
First step: using the least-squares approach to interpret the gravity anomaly yielding from the above mentioned parameters for different s-values for the three suggested shape bodies, i.e.,
Moreover, a 10% random noise added to the synthetic gravity data mentioned above (Figure 2) to test the efficiency of the suggested approaches in interpreting the gravity data. Also, the three approaches are used for this data as mentioned in Table 1 (the least-squares approach results), Figure 3 (Werner deconvolution results), and Table 2 (the particle swarm results). Finally, the estimated parameters are in all case are in good agreement with the true parameters.
3.2 Self-potential anomaly model
A self-potential anomaly of a horizontal cylinder model is generated using the following parameters
4. Field examples
The three suggested approaches have been organized to inspect the gravity and self-potential anomalies due to three simple bodies of various structures, e.g., sheets, cylinders, and spheres. Two mineral field examples from India and Turkey have been interpreted to study the reliability of the suggested approaches. The relevant model parameters (
4.1 Gravity anomaly of manganese ore body
Figure 6 shows a gravity anomaly was collected over a manganese deposit near Nagpur, India . This gravity profile has a length of 333 m and digitized with an interval of 27 m. This gravity anomaly is subjected to the three interpretation approaches as discussed earlier. Firstly, the interpreted results due to applying the least-squares approach are shown in Table 5 for various s-values. Besides, the use of Werner deconvolution approach, the interpretive results obtained are
4.2 Self-potential anomaly of manganese ore body
Figure 7 demonstrates a self-potential anomaly over a Weiss copper ore body in the Ergani copper district, Turkey . The Weiss self-potential anomaly profile has a length of 144 m and digitized with an interval of 7.7 m. This anomaly has subjected to the three interpretation approaches as discussed earlier. Firstly, the interpreted results due to applying the least-squares approach are shown in Table 7 for various s-values. Also, the applying of the Werner deconvolution approach, the interpretive results obtained are
The three geophysical approaches (the least-squares approach, Werner deconvolution approach, and the particle swarm approach) discussed here to interpret gravity or self-potential data using a combined formula for the simple geometric models (spheres, cylinders, and dikes) are stable and give a good results. The stability of these approaches has been confirmed and tested applying two synthetic examples with a 10% and without random noise and two field data for mineral explorations. The estimated parameters in all cases demonstrated the importance of these approaches in interpreting the gravity or self-potential data.
The authors would like to thank and express appreciation to Ms. Dolores Kuzelj, Author Service Manager, for her assistance and cooperation in this issue.
Conflict of interest
There is no conflict of interest.