3D Airborne EM Forward Modeling Based on Finite-Element Method with Goal-Oriented Adaptive Octree Mesh

• Published in: Remote sensing (Basel, Switzerland) Vol. 15; no. 11; p. 2816
• Main Authors: Han, Xue, Ni, Jianfu, Yin, Changchun, Zhang, Bo, Huang, Xin, Zhu, Jiao, Liu, Yunhe, Ren, Xiuyan, Su, Yang
• Format: Journal Article
• Language: English
• Published: Basel MDPI AG 29.05.2023
• Subjects: 3D forward modeling; Accuracy; Adaptive algorithms; adaptive octree mesh; airborne EM; Algorithms; Comparative analysis; Computer applications; Computing costs; Efficiency; Electric fields; Errors; Exact solutions; Finite element method; frequency domain; Geological structures; Geological surveys; Half spaces; Hilbert space; Mathematical functions; Mesh generation; Methods; Modelling; Numerical analysis; Octree encoding; Octrees; Remote sensing; surveys; China
• ISSN: 2072-4292; 2072-4292
• DOI: 10.3390/rs15112816

The finite-element (FE) method for three-dimensional (3D) airborne electromagnetic (AEM) modeling can flexibly simulate complex geological structures at high accuracy. However, it has low efficiency and high computational requirements. To solve these problems, one needs to generate meshes more reasonably. In view of this, we develop an adaptive octree meshing scheme for frequency-domain AEM modeling. The octree meshes have the characteristics of regularity and flexibility, while the adaptive algorithm can effectively refine the mesh locally. In our adaptive mesh generation, the posterior errors and weighted coefficients are used to construct the final weighted posterior errors. We verify the accuracy of our method by comparing its results with semi-analytical solutions for a half-space model. Furthermore, we use the spectral-element (SE) method and our method to calculate EM responses for an abnormal block model and compare their computational costs. The results show that our adaptive scheme has obviously technical advantages over SE method for AEM modeling with multiple frequencies and multiple survey stations. Finally, we calculate a model with complex geological structures to verify the feasibility of our algorithm in complex geological circumstances.