Lp-Norm Inversion of Gravity Data Using Adaptive Differential Evolution

As a popular population based heuristic evolutionary algorithm, differential evolution (DE) has been widely applied in various science and engineering problems. Similar to other global nonlinear algorithms, such as genetic algorithm, simulated annealing, particle swarm optimization, etc., the DE algorithm is mostly applied to resolve the parametric inverse problem, but has few applications in physical property inversion. According to our knowledge, this is the first time DE has been applied in obtaining the physical property distribution of gravity data due to causative sources embedded in the subsurface. In this work, the search direction of DE is guided by better vectors, enhancing the exploration efficiency of the mutation strategy. Besides, to reduce the over-stochastic of the DE algorithm, the perturbation directions in mutation operations are smoothed by using a weighted moving average smoothing technique, and the Lp-norm regularization term is implemented to sharpen the boundary of density distribution. Meanwhile, in the search process of DE, the effect of Lp-norm regularization term is controlled in an adaptive manner, which can always have an impact on the data misfit function. In the synthetic anomaly case, both noise-free and noisy data sets are considered. For the field case, gravity anomalies originating from the Shihe iron ore deposit in China were inverted and interpreted. The reconstructed density distribution is in good agreement with the one obtained by drill-hole information. Based on the tests in the present study, one can conclude that the Lp-norm inversion using DE is a useful tool for physical property distribution using gravity anomalies.