Inversion of electromagnetic induction data using a novel wavelet-based and scale-dependent regularization term

TL;DR

This paper introduces a wavelet-based inversion method with a scale-dependent regularization for electromagnetic data, allowing adaptive and sparse conductivity profile reconstructions that better match non-smooth geological features.

Contribution

A novel wavelet transform-based inversion scheme with sparsity constraints, improving stability and adaptability for non-smooth conductivity profiles in electromagnetic data inversion.

Findings

Successfully applied to FDEM dataset

Produces more accurate, less smoothed conductivity profiles

Flexible choice of wavelet basis for different profiles

Abstract

The inversion of electromagnetic induction data to a conductivity profile is an ill-posed problem. Regularization improves the stability of the inversion and, based on Occam's razor principle, a smoothing constraint is typically used. However, the conductivity profiles are not always expected to be smooth. Here, we develop a new inversion scheme in which we transform the model to the wavelet space and impose a sparsity constraint. This sparsity constrained inversion scheme will minimize an objective function with a least-squares data misfit and a sparsity measure of the model in the wavelet domain. A model in the wavelet domain has both temporal as spatial resolution, and penalizing small-scale coefficients effectively reduces the complexity of the model. Depending on the expected conductivity profile, an optimal wavelet basis function can be chosen. The scheme is thus adaptive. Finally, we apply this new scheme on a frequency domain electromagnetic sounding (FDEM) dataset, but the scheme could equally apply to any other 1D geophysical method.