Numerical results using the least-squares approach for a gravity model due to horizontal cylinder without and with a 10% of random noise (* K*= 200 mGal×m,

*= 5 m,*z

*= 0,*x

_{o}

*= 1, and profile length = 100 m).*q

Open access peer-reviewed chapter

By Khalid S. Essa and Mahmoud Elhussein

Submitted: January 10th 2020Reviewed: March 16th 2020Published: April 20th 2020

DOI: 10.5772/intechopen.92139

Downloaded: 80

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
- noise
- 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 [33], memory-based hybrid dragonfly [34], differential evolution [35, 36].

Here, a combined formula for both gravity and self-potential [37] 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 (* K*), the polarization angle (

Advertisement## 2. The suggested combined gravity or self-potential formula

g x i z q = K z m x i − x o 2 + z 2 q . E1V x i z θ q = K xcosθ + zsinθ x i − x o 2 + z 2 q . E2J x i = K c x i cosθ n + z p sinθ m x i − x o 2 + z 2 q , E3### 2.1 The least-squares approach

### 2.2 Werner deconvolution approach

J x i x i − x o 2 + z 2 q − Kc x i cosθ n + Kz p sinθ m = 0 , E4J x i x i − e 1 2 + e 2 q − e 3 x i + e 4 = 0 , E5### 2.3 The particle swarm approach

V i k + 1 = c 3 V i k + c 1 rand P best − x i k + 1 + c 2 rand J best − x i k + 1 , E6x i k + 1 = x i k + V i k + 1 , E7Obj = ∑ j = 1 N J j o − J j c 2 N , E8

Firstly, the gravity anomaly formula due to simple geometric shapes is [15, 16, 18]

Secondly, the self-potential anomaly formula for the same simple geometric models is [14]

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:

where * K*is the amplitude coefficient, which depends on the shape of the buried model,

Essa [37] developed this approach, which was relied on solving the problem of finding the depth from the measured data by solving a nonlinear form * F*(

Werner deconvolution was proposed by Werner in 1953 [38]. 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:

where

Eq. (5) is linear form in the four variables _{1}, _{2}, _{3}, and _{4}, so that a mathematically unique solution can be found for them from evaluating the equation at four points by assuming the shape of the buried structure.

The particle swarm was suggested by [39] 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 (* P*) at every iteration until reach the optimum location (

where * i*at the

The five source parameters (* K*,

where * N*is the data points number,

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.

A gravity anomaly of a horizontal cylinder model is generated using the following parameters * K*= 200 mGal×m,

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., * q*= 0.5,

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.

A self-potential anomaly of a horizontal cylinder model is generated using the following parameters * K*= 200 mV×m,

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 (* K*,

Figure 6 shows a gravity anomaly was collected over a manganese deposit near Nagpur, India [45]. 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 * z*= 56.8 m and

Figure 7 demonstrates a self-potential anomaly over a Weiss copper ore body in the Ergani copper district, Turkey [46]. 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 * z*= 36.9 m and

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.

Advertisement## Acknowledgments

The authors would like to thank and express appreciation to Ms. Dolores Kuzelj, Author Service Manager, for her assistance and cooperation in this issue.

There is no conflict of interest.

Share on my website

Save to Mendeley

Read in ReadCube

Share on TwitterShare on FacebookShare on LinkedIn

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

80total chapter downloads

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reportingEdited by Khalid S. Essa

Edited by Khalid S. Essa

First chapter#### Introductory Chapter: Mineral Exploration from the Point of View of Geophysicists

By Khalid S. Essa and Marc Munschy

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us