The present invention patent discloses a method for high-order finite element forward modeling of
magnetotellurics based on
hybrid grids, which specifically relates to the field of geophysical electromagnetic three-dimensional forward modeling. The
Delaunay triangulation algorithm is used to divide the earth's surface into triangular grids, and triangular grids are extended on the triangular grids to form
triangular prism units; the areas adjacent to the
triangular prism units are discretized into tetrahedral units, and the underground
electrical model is discretized using the constructed
hybrid grid. Finite element unit analysis is performed through high-order shape functions to discretize the magnetotelluric wave equations and form a magnetotelluric high-order finite element
system; the
hybrid grid finite element
system is solved to obtain the magnetotelluric response
signal. The technical solution of the present invention solves the problem of increased computational cost caused by the use of traditional unstructured networks in MT forward modeling under high-frequency conditions. The hybrid grid can adapt to complex geoelectric models with strong
terrain fluctuations, which requires less computational cost than the use of traditional unstructured units.