Gravity profiles interpretation applying a metaheuristic particle optimization algorithm of mineralized bodies resembled by finite elements.

The interpretation of gravity anomalies is crucial for identifying subsurface mineralized targets and understanding the density variations between the targets and the surrounding structures. To confirm the presence of ore and mineral targets, simple geometric bodies are often used. One of the commonly used global metaheuristic algorithms for gravity data analysis is the particle optimization algorithm. In this study, we employed this method to determine the parameters of buried bodies that resemble finite vertical cylinders by inferring gravity anomalies profiles (amplitude coefficient, depth to top, depth to bottom, origin, and length of the target representing the difference between two depths). The algorithm utilizes particle movement to identify the best way to reach the global or optimum solution. The algorithm's performance was evaluated on synthetic-examples with and without noise (5 % and 10 % levels) and also verified on a real dataset for mineral exploration from Canada. The results showed that the algorithm's stability and accuracy were not affected by the presence of noise and multi-models. Moreover, the field case results were consistent with the existing geological information, borehole data, and previously published outcomes.

A new inversion method of estimation of simultaneous near surface bulk density variations and terrain correction across the Bandar Charak (Hormozgan-Iran).

A gravity inversion method based on the Nettleton-Parasnis technique is used to estimate near surface density in an area without exposed outcrop or where outcrop occurrences do not adequately represent the subsurface rock densities. Its accuracy, however, strongly depends on how efficiently the regional trends and very local (terrain) effects are removed from the gravity anomalies processed. Nettleton's method implemented in a usual inversion scheme and combined with the simultaneous determination of terrain corrections. This method may lead to realistic density estimations of the topographical masses. The author applied this technique in the Bandar Charak (Hormozgan-Iran) with various geological/geophysical properties. These inversion results are comparable to both values obtained from density logs in the mentioned area and other methods like Fractal methods. The calculated densities are 2.4005 gr/cm3. The slightly higher differences between calculated densities and densities of the hand rock samples may be caused by the effect of sediment-filled valleys.

Comparison result of inversion of gravity data of a fault by particle swarm optimization and Levenberg-Marquardt methods.

The purpose of this study was to compare the performance of two methods for gravity inversion of a fault. First method [Particle swarm optimization (PSO)] is a heuristic global optimization method and also an optimization algorithm, which is based on swarm intelligence. It comes from the research on the bird and fish flock movement behavior. Second method [The Levenberg-Marquardt algorithm (LM)] is an approximation to the Newton method used also for training ANNs. In this paper first we discussed the gravity field of a fault, then describes the algorithms of PSO and LM And presents application of Levenberg-Marquardt algorithm, and a particle swarm algorithm in solving inverse problem of a fault. Most importantly the parameters for the algorithms are given for the individual tests. Inverse solution reveals that fault model parameters are agree quite well with the known results. A more agreement has been found between the predicted model anomaly and the observed gravity anomaly in PSO method rather than LM method.

Gravity inversion of a fault by Particle swarm optimization (PSO).

Particle swarm optimization is a heuristic global optimization method and also an optimization algorithm, which is based on swarm intelligence. It comes from the research on the bird and fish flock movement behavior. In this paper we introduce and use this method in gravity inverse problem. We discuss the solution for the inverse problem of determining the shape of a fault whose gravity anomaly is known. Application of the proposed algorithm to this problem has proven its capability to deal with difficult optimization problems. The technique proved to work efficiently when tested to a number of models.