A Hybrid PCG-Bat Algorithm for 2D Gravity Inversion: Applications for Ore Deposits Exploration and Interpretation of Sedimentary Basins

2021 ◽  
pp. 104497
Author(s):  
Mohamed Gobashy ◽  
Maha Abdelazeem ◽  
Mohamed Abdrabou ◽  
Mohamed H. Khalil
2020 ◽  
Author(s):  
Mohamed Abdrabou ◽  
Maha Abdelazeem ◽  
Mohamed Gobashy

<p>Geophysical data such as gravity data can be inverted to get a subsurface image, which depicts the subsurface distribution of physical property. Consequently, inversion of geophysical data has an effective role for interpreting measured geophysical anomalies in hydrocarbons and mineral applications. Interest about ore deposits exploration and sedimentary basins interpretation is associated with their economic importance. The presence of sedimentary basins gives lower amplitude of gravity anomalies with negative signals, due to the negative density contrast as these sedimentary basins have lower density than that of the neighboring basement rocks. In prospecting ore deposits, studying the spatial distributions of densities in the subsurface is essential of significance.Two dimensional forward modelling strategy can be done via locating the rectangular cells with fixed size directly underneath the location of the observed data points using regular grid discretization. Density vector of the subsurface rectangular cells are obtained via solving the 2D gravity inverse problem by optimizing an objective function (i.e., the differences between observed and inverted residual gravity data sets). In this work, a hybrid algorithm merging a bat (BAT) algorithm with the preconditioned conjugate gradient (PCG) method is suggested as a mean for inverting surface gravity anomalies to obtain the density distribution in the subsurface. Like the hybrid, minimization algorithm has the capability to make use of the advantages of both two techniques. In this hybrid algorithm, the BAT algorithm was utilized to construct an initial solution for the PCG technique. The BAT optimizer acts as a rapid build-up of the model, whereas the second modifies the finer model approximated solution. This modern algorithm was firstly applied on a free-noise synthetic data and to a noisy data with three different levels of random noise, and good results obtained through the inversion. The validity and applicability of our algorithm are applied to real residual gravity anomalies across the San Jacinto graben in southern California, USA, and Sierra Mayor - Sierra Pinta graben, USA and prospecting of the Poshi Cu-Ni deposits, Xinjiang, northwest China. The obtained results are in excellent accordance with those produced by researchers in the published literature.</p><p> </p><p><strong>Keywords: </strong>Gravity data, 2D Inversion, BAT algorithm, Preconditioned Conjugate Gradient, Sedimentary Basins.</p>


1991 ◽  
Vol 33 (8) ◽  
pp. 822-829 ◽  
Author(s):  
D. I. Pavlov ◽  
D. I. Gorzhevskiy ◽  
G. A. Goleva ◽  
M. K. Kalinko ◽  
A. A. Kartsev ◽  
...  

Geophysics ◽  
2017 ◽  
Vol 82 (2) ◽  
pp. G23-G34 ◽  
Author(s):  
João B. C. Silva ◽  
Darcicléa F. Santos

We have developed a novel approach to compute, in an efficient and versatile way, the gravity anomaly produced by an arbitrary, discrete 3D distribution of density contrast. The method allows adjustable precision and is particularly suited for the interpretation of sedimentary basins. Because the gravity field decays with the square of the distance, we use a discrete Green’s operator that may be made much smaller than the whole study area. For irregularly positioned observations, this discrete Green’s operator must be computed just at the first iteration, and because each of its horizontal layers presents a center of symmetry, only one-eighth of its total elements need to be calculated. Furthermore, for gridded data on a plane, this operator develops translation symmetry for the whole region of interest and has to be computed just once for a single arbitrary observation position. The gravity anomaly is obtained as the product of this small operator by any arbitrary density contrast distribution in a convolution-like operation. We use the proposed modeling to estimate the basement relief of a [Formula: see text] basin with density contrast varying along [Formula: see text] only using a smaller Green’s operator at all iterations. The median of the absolute differences between relief estimates, using the small and a large operator (the latter covering the whole basin) has been approximately 9 m for a 3.6 km deep basin. We also successfully inverted the anomaly of a simulated basin with density contrast varying laterally and vertically, and a real anomaly produced by a steeply dipping basement. The proposed modeling is very fast. For 10,000 observations gridded on a plane, the inversion using the proposed approach for irregularly spaced data is two orders of magnitude faster than using an analytically derived fitting, and this ratio increases enormously with the number of observations.


2019 ◽  
Vol 122 (4) ◽  
pp. 397-420 ◽  
Author(s):  
L.D. Ashwal

Abstract On the last page of his 1937 book “Our Wandering Continents” Alex Du Toit advised the geological community to develop the field of “comparative geology”, which he defined as “the study of continental fragments”. This is precisely the theme of this paper, which outlines my research activities for the past 28 years, on the continental fragments of the Indian Ocean. In the early 1990s, my colleagues and I were working in Madagascar, and we recognized the need to appreciate the excellent geological mapping (pioneered in the 1950s by Henri Besairie) in a more modern geodynamic context, by applying new ideas and analytical techniques, to a large and understudied piece of continental crust. One result of this work was the identification of a 700 to 800 Ma belt of plutons and volcanic equivalents, about 450 km long, which we suggested might represent an Andean-type arc, produced by Neoproterozoic subduction. We wondered if similar examples of this magmatic belt might be present elsewhere, and we began working in the Seychelles, where late Precambrian granites are exposed on about 40 of the >100 islands in the archipelago. Based on our new petrological, geochemical and geochronological measurements, we built a case that these ~750 Ma rocks also represent an Andean-type arc, coeval with and equivalent to the one present in Madagascar. By using similar types of approaches, we tracked this arc even further, into the Malani Igneous Province of Rajasthan, in northwest India. Our paleomagnetic data place these three entities adjacent to each other at ~750 Ma, and were positioned at the margins, rather than in the central parts of the Rodinia supercontinent, further supporting their formation in a subduction-related continental arc. A widespread view is that in the Neoproterozoic, Rodinia began to break apart, and the more familiar Gondwana supercontinent was assembled by Pan-African (~500 to 600 Ma) continental collisions, marked by the highly deformed and metamorphosed rocks of the East African Orogen. It was my mentor, Kevin Burke, who suggested that the present-day locations of Alkaline Rocks and Carbonatites (called “ARCs”) and their Deformed equivalents (called “DARCs”), might mark the outlines of two well-defined parts of the Wilson cycle. We can be confident that ARCs formed originally in intracontinental rift settings, and we postulated that DARCs represent suture zones, where vanished oceans have closed. We also found that the isotopic record of these events can be preserved in DARC minerals. In a nepheline syenite gneiss from Malawi, the U-Pb age of zircons is 730 Ma (marking the rifting of Rodinia), and that of monazites is 522 Ma (marking the collisional construction of Gondwana). A general outline of how and when Gondwana broke apart into the current configuration of continental entities, starting at about 165 Ma, has been known for some time, because this record is preserved in the magnetic properties of ocean-floor basalts, which can be precisely dated. A current topic of active research is the role that deep mantle plumes may have played in initiating, or assisting, continental fragmentation. I am part of a group of colleagues and students who are applying complementary datasets to understand how the Karoo (182 Ma), Etendeka (132 Ma), Marion (90 Ma) and Réunion (65 Ma) plumes influenced the break-up of Gondwana and the development of the Indian Ocean. Shortly after the impingement of the Karoo plume at 182 Ma, Gondwana fragmentation began as Madagascar + India + Antarctica separated from Africa, and drifted southward. Only after 90 Ma, when Madagascar was blanketed by lavas of the Marion plume, did India begin to rift, and rapidly drifted northward, assisted by the Marion and Deccan (65 Ma) plumes, eventually colliding with Asia to produce the Himalayas. It is interesting that a record of these plate kinematics is preserved in the large Permian – Eocene sedimentary basins of western Madagascar: transtensional pull-apart structures are dextral in Jurassic rocks (recording initial southward drift with respect to Africa), but change to sinistral in the Eocene, recording India’s northward drift. Our latest work has begun to reveal that small continental fragments are present in unexpected places. In the young (max. 9 Ma) plume-related, volcanic island of Mauritius, we found Precambrian zircons with ages between 660 and 3000 Ma, in beach sands and trachytic lavas. This can only mean that a fragment of ancient continent must exist beneath the young volcanoes there, and that the old zircons were picked up by ascending magmas on their way to surface eruption sites. We speculate, based on gravity inversion modelling, that continental fragments may also be present beneath the Nazareth, Saya de Malha and Chagos Banks, as well as the Maldives and Laccadives. These were once joined together in a microcontinent we called “Mauritia”, and became scattered across the Indian Ocean during Gondwana break-up, probably by mid-ocean ridge “jumps”. This work, widely reported in international news media, allows a more refined reconstruction of Gondwana, suggests that continental break-up is far more complex than previously perceived, and has important implications for regional geological correlations and exploration models. Our results, as interesting as they may be, are merely follow-ups that build upon the prescient and pioneering ideas of Alex Du Toit, whose legacy I appreciatively acknowledge.


2017 ◽  
Vol 6 (1) ◽  
pp. 193-198 ◽  
Author(s):  
Kunal Kishore Singh ◽  
Upendra Kumar Singh

Abstract. Particle swarm optimization (PSO) is a global optimization technique that works similarly to swarms of birds searching for food. A MATLAB code in the PSO algorithm has been developed to estimate the depth to the bottom of a 2.5-D sedimentary basin and coefficients of regional background from observed gravity anomalies. The density contrast within the source is assumed to vary parabolically with depth. Initially, the PSO algorithm is applied on synthetic data with and without some Gaussian noise, and its validity is tested by calculating the depth of the Gediz Graben, western Anatolia, and the Godavari sub-basin, India. The Gediz Graben consists of Neogen sediments, and the metamorphic complex forms the basement of the graben. A thick uninterrupted sequence of Permian–Triassic and partly Jurassic and Cretaceous sediments forms the Godavari sub-basin. The PSO results are better correlated with results obtained by the Marquardt method and borehole information.


Sign in / Sign up

Export Citation Format

Share Document