Edge finite-element modeling using unstructured meshes: the 3D MT example

摘要

Three-dimension (3D) forward modeling is a challenging problem for geometrically complex MT problems. We presented a new edge-based finite-element algorithm cooperated with unstructured mesh for accurately and efficiently simulating 3D MT responses. The curl-curl equation of electric field in frequency was used to deduce the H(curl) variation weak form; and then, Galerkin weighting rule was performed on the linear unstructured edge-based finite-element space to obtain a sparse complex linear equation; finally, BI-CGSTAB solver with ILU preconditioner was employed to offer the unknown electric fields. A local mesh refinement technique in the neighbor of measuring MT stations was used to distinctly update the accuracies of numerical solutions. Two synthetic models validated the powerful performances of our algorithms. We believe our method should have more legs in more complex MT cases.