LAUSR

The conventional, L2-norm-based, regularization term in electromagnetic (EM) inversions implements smooth constraints on model complexity in the space domain, which can smoothen the boundaries of complex underground structures. To improve the resolution of 3-D frequency-domain airborne EM (AEM) inversions, we propose a new algorithm for sparse-regularized inversion based on the shearlet transform. Unlike traditional methods that invert the model parameters in the space domain, we first transform the 3-D resistivity model into the frequency domain and then invert the sparse coefficients using an L1-norm measure to ensure the sparseness of the solution. Finally, we transform the shearlet coefficients back to the space domain to update the model. The shearlet transform has inherent multiscale and multidirectional properties, making it capable of effectively extracting complex geometries such as curved boundaries. We adopt the finite-difference method and the iteratively reweighted least-squares scheme for our 3-D AEM modeling and inversions and apply the “moving footprint” technique to speed up the inversion. Tests using synthetic data show that sparse-regularized inversion based on the shearlet transform can obtain more-focused inversion results than conventional smoothness-constrained inversions based on the L2-norm. Tests using field survey data also reveal that the new method can achieve more realistic underground structures.