Hybrid mesh for magnetotelluric forward modeling based on the finite element method

Unstructured tetrahedral grids have been applied in magnetotelluric (MT) forward modeling using the finite element (FE) method because of their adaptability to complex anomalies. However, high-quality results require an extreme refinement of the near-surface area, which leads to excessive meshes and an increased degree of freedom (DoF) of the governing equation of the finite element system. To reduce the computational cost, we have developed a hybrid mesh based on triangular prisms and tetrahedrons. The required elements in the near-surface area are reduced because the quality of the triangular prism is not limited by the element aspect ratio. The deep area is discretized by tetrahedral elements to ensure the flexibility of the unstructured grids. The superiority of this hybrid mesh has been tested on a layered model, the DTM1 model and terrain relief models. The results show that the modeling efficiency has been improved, especially for high-frequency data. The accuracy of the modeling using the hybrid mesh is significantly higher than that of the tetrahedral mesh with a similar DoF. Usage of the hybrid mesh can be easily adapted to complex geoelectric models with strong terrain fluctuations, which requires less computational cost than using conventional unstructured elements.