Fast 3D time-domain airborne EM forward modeling using random under-sampling

Abstract The high sampling rate of airborne electromagnetic (AEM) systems can create huge data volumes, causing major challenges for three-dimensional (3D) electromagnetic modeling and inversions and constraining the practical applicability of AEM. Rapid and accurate forward modeling is the key to efficient 3D inversion. Here a new strategy for improving the efficiency of 3D time-domain AEM forward modeling is presented, based on compressed-sensing theory combined with the finite-volume method. A Poisson disk random under-sampling method was used to obtain random survey stations and to calculate AEM responses via the finite-volume method at these sampling stations. Coefficients from the curvelet transform with anisotropic characteristics for sparse expression were used to analyze the sparsity of AEM data, with the AEM responses of all survey stations being obtained by reconstruction. Key factors influencing the accuracy of the reconstruction are also discussed, including the sampling method and rate, the sparse transformation, and the noise. To test the effectiveness of the method, reconstructed results were compared with those from traditional schemes for forward modeling using all stations. Results indicate that the new strategy requires only 20% ~ 50% of samples to achieve high-precision reconstruction of the AEM signal even for complex geoelectrical structures, and the efficiency of 3D forward modeling can be enhanced by a factor of two to five, laying a foundation for fast AEM inversions of large datasets and complex models.

[1]  Lawrence Carin,et al.  Bayesian Compressive Sensing , 2008, IEEE Transactions on Signal Processing.

[2]  Xiaoming Huo,et al.  Uncertainty principles and ideal atomic decomposition , 2001, IEEE Trans. Inf. Theory.

[3]  Xiaobo Zhang,et al.  A survey on one-bit compressed sensing: theory and applications , 2018, Frontiers Comput. Sci..

[4]  K. Yee Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[5]  Glenn A. Wilson,et al.  2.5D inversion of airborne electromagnetic data , 2006 .

[6]  Boreom Lee,et al.  Online compressive covariance sensing , 2019, Signal Process..

[7]  E.J. Candes,et al.  An Introduction To Compressive Sampling , 2008, IEEE Signal Processing Magazine.

[8]  Xiaohong Chen,et al.  Seismic data reconstruction based on CS and Fourier theory , 2013, Applied Geophysics.

[9]  Jean-Luc Starck,et al.  Sparse Solution of Underdetermined Systems of Linear Equations by Stagewise Orthogonal Matching Pursuit , 2012, IEEE Transactions on Information Theory.

[10]  Emmanuel J. Candès,et al.  Quantitative Robust Uncertainty Principles and Optimally Sparse Decompositions , 2004, Found. Comput. Math..

[11]  G. Vignoli,et al.  New regional stratigraphic insights from a 3D geological model of the Nasia sub-basin, Ghana, developed for hydrogeological purposes and based on reprocessed B-field data originally collected for mineral exploration , 2019, Solid Earth.

[12]  Douglas W. Oldenburg,et al.  Three dimensional inversion of multisource time domain electromagnetic data , 2013 .

[13]  David L Donoho,et al.  Compressed sensing , 2006, IEEE Transactions on Information Theory.

[14]  Laurent Demanet,et al.  Fast Discrete Curvelet Transforms , 2006, Multiscale Model. Simul..

[15]  Gitta Kutyniok,et al.  1 . 2 Sparsity : A Reasonable Assumption ? , 2012 .

[16]  F. Herrmann,et al.  Simply denoise: Wavefield reconstruction via jittered undersampling , 2008 .

[17]  B. Zhang,et al.  3D time-domain airborne electromagnetic inversion based on secondary field finite-volume method , 2018, GEOPHYSICS.

[18]  Michael S. Zhdanov,et al.  3D inversion of airborne electromagnetic data , 2012 .

[19]  R. Abma,et al.  3D interpolation of irregular data with a POCS algorithm , 2006 .

[20]  Cui Ziming,et al.  Application of Wavelet Transform Data Reconstruction in Seismic Data Interpretation , 2018 .

[21]  D. Beaty,et al.  The Carbonate-Hosted Lisheen Zn-Pb-Ag Deposit, County Tipperary, Ireland , 2002 .

[22]  Felix J. Herrmann,et al.  Randomized sampling and sparsity: Getting more information from fewer samples , 2010 .

[23]  I Frerichs,et al.  Compressive sensing in electrical impedance tomography for breathing monitoring , 2019, Physiological measurement.

[24]  M. Lustig,et al.  Compressed Sensing MRI , 2008, IEEE Signal Processing Magazine.

[25]  Christian Jutten,et al.  Decoding real-field codes by an iterative Expectation-Maximization (EM) algorithm , 2008, 2008 IEEE International Conference on Acoustics, Speech and Signal Processing.

[26]  Yunhe Liu,et al.  Footprint for frequency-domain airborne electromagnetic systems , 2014 .

[27]  Mike E. Davies,et al.  Iterative Hard Thresholding for Compressed Sensing , 2008, ArXiv.

[28]  Michael P. Poland,et al.  Toward continuous 4D microgravity monitoring of volcanoes , 2008 .

[29]  Michael S. Zhdanov,et al.  3D inversion of airborne electromagnetic data using a moving footprint , 2010 .

[30]  Douglas W. Oldenburg,et al.  3-D inversion of airborne electromagnetic data parallelized and accelerated by local mesh and adaptive soundings , 2014 .

[31]  Bhaskar D. Rao,et al.  Sparse signal reconstruction from limited data using FOCUSS: a re-weighted minimum norm algorithm , 1997, IEEE Trans. Signal Process..