3-D gravity inversion for the basement relief reconstruction through modified success-history-based adaptive differential evolution

Abstract

SUMMARY A gravity inversion procedure using the success-history-based adaptive differential evolution (SHADE) algorithm is presented to reconstruct the 3-D basement relief geometry in sedimentary basins. We introduced exponential population size (number) reduction (EPSR) to reduce the computational cost and used self-adaptive control parameters to solve this highly nonlinear inverse problem. Model parametrization was carried out by discretizing the sedimentary cover via juxtaposed right prisms, each placed below each observation point. Resolvability characteristics of the 3-D inverse problem were revealed through some cost function topography landscapes. The fine-tuned control parameter namely, population number allowed us to get best benefit from the algorithm. Additionally, a stabilizing function as a relative constraint was used to avoid undesired effects originated from the ill-posedness of the problem. In the synthetic data cases, the strategy we propose outperformed the linear population number reduction strategy which has won various IEEE–CEC competitions so far. Thorough uncertainty assessments via probability density function and principal component analysis demonstrated the solidity of the obtained inverse models. In the real data case, residual gravity anomalies of two well-known major grabens of Aegean Graben System (Türkiye), calculated thanks to the finite element method, were inverted. It was determined that the inverse solutions obtained for these basement reliefs, whose depths are still controversial, are statistically reliable. Moreover, these depths were found to be less than the depths reported in most previous studies. We conclude that the SHADE using EPSR strategy that we propose is a powerful alternative inversion tool for highly nonlinear geophysical problems.