LAUSR

Three-dimensional (3D) airborne electromagnetic (AEM) inversions have been restricted by the modeling efficiency resulted from the complex geology in exploration areas and massive amount of data collected by AEM systems. In order to improve the modeling efficiency, we develop an algorithm that combines the hexahedral vector finite-element (FE) with octree meshes, in which the boundary conditions are imposed via an algebraic constraint to ensure the continuity of the FE solution. This makes the division with hexahedral meshes more flexible for complex geology like rugged topography or underground structures, so that we can reduce the number of elements while maintaining the accuracy. After formulating the forward problem, we check the accuracy of our algorithm by taking a homogeneous half-space model and comparing the results of our octree method with the semi-analytical solutions. Furthermore, we demonstrate the efficiency of our octree method by comparing with the traditional FE method using tetrahedral meshes. Finally, we subdivide a complex topography constructed using two-dimensional (2D) Gaussian rough surface and calculate the EM responses with and without anomaly embedded. The results show that the EM responses are overwhelmed by the earth topography. We carry out the topographic correction by taking a method based on the ratio of EM responses with and without anomaly. The experiments show that after topographic correction to AEM data, the response of anomaly becomes more distinguishable, so that the anomaly can be clearly identified. Furthermore, we also calculate the EM response for a realistic model — the Ovoid Zone ore body located at Voisey’s Bay, Labrador, Canada, to verify the flexibility and practicality of our algorithm.