Multiscale 3-D Stochastic Inversion of Frequency-Domain Airborne Electromagnetic Data

In mineral, environmental, and engineering explorations, we frequently encounter geological bodies with varied sizes, depths, and conductivity contrasts with surround rocks and try to interpret them with single survey data. The conventional three-dimensional (3-D) inversions significantly rely on the size of the grids, which should be smaller than the smallest geological target to achieve a good recovery to anomalous electric conductivity. However, this will create a large amount of unknowns to be solved and cost significant time and memory. In this paper, we present a multi-scale (MS) stochastic inversion scheme based on shearlet transform for airborne electromagnetic (AEM) data. The shearlet possesses the features of multi-direction and multi-scale, allowing it to effectively characterize the underground conductivity distribution in the transformed domain. To address the practical implementation of the method, we use a compressed sensing method in the forward modeling and sensitivity calculation, and employ a preconditioner that accounts for both the sampling rate and gradient noise to achieve a fast stochastic 3-D inversion. By gradually updating the coefficients from the coarse to fine scales, we obtain the multi-scale information on the underground electric conductivity. The synthetic data inversion shows that the proposed MS method can better recover multiple geological bodies with different sizes and depths with less time consumption. Finally, we conduct 3-D inversions of a field dataset acquired from Byneset, Norway. The results show very good agreement with the geological information.