No abstract

We present a wavelet finite-element method (WFEM) based on B-spline wavelets on the interval (BSWI) for three-dimensional (3D) frequency-domain airborne EM modeling using a secondary coupled-potential formulation. The BSWI, which is constructed on the interval (0, 1) by joining piecewise B-spline polynomials between nodes together, has proved to have excellent numerical properties of multiresolution and sparsity and thus is utilized as the basis function in our WFEM. Compared to conventional basis functions, the BSWI is able to provide higher interpolating accuracy and boundary stability. Furthermore, due to the sparsity of the wavelet, the coefficient matrix obtained by BSWI-based WFEM is sparser than that formed by general FEM, which can lead to shorter solution time for the linear equations system. To verify the accuracy and efficiency of our proposed method, we ran numerical experiments on a half-space model and a layered model and compared the results with one-dimensional (1D) semi-analytic solutions and those obtained from conventional FEM. We then studied a synthetic 3D model using different meshes and BSWI basis at different scales. The results show that our method depends less on the mesh, and the accuracy can be improved by both mesh refinement and scale enhancement.

Summary Geomagnetic depth sounding (GDS) is a geophysical electromagnetic (EM) method that studies the deep structure and composition of the earth by using long-period EM signals from geomagnetic observatories and satellites. In this paper, a three-dimensional (3D) anisotropic GDS modeling algorithm is developed. The curl-curl equation is discretized using the edge-based finite-element (FE) method on unstructured tetrahedral grids. In order to solve the computationally demanding problem of EM modeling on a global scale, the complex linear system is first separated into the equivalent real linear systems, and then the real system is iteratively solved by the flexible generalized minimum residual (FGMRES) method with a block diagonal preconditioner. This will greatly reduce the condition number of the linear system and thus speed up the solution process. We verify the accuracy of the proposed algorithm by comparing our results with the existing methods. After that, we design a subduction zone model to simulate the EM responses under isotropic and anisotropic environments, respectively. The numerical results show the high efficiency of the proposed algorithm and the response differences between isotropic and anisotropic models. This research can provide theoretical and technical support for the high-accuracy and efficient inversion of GDS data for the geo-dynamic study.