This invention provides a numerical
simulation method for
magnetotellurics under the Lorentz specification. The Maxwell equations are transformed into a set of governing equations consisting only of vector potentials using the Lorentz specification. A two-dimensional
Fourier transform is performed on the quadratic
vector potential governing equations in the horizontal direction, converting the three-dimensional governing equations into one-dimensional governing equations, reducing computational load and storage requirements, and improving computational efficiency. The
Fourier transform can employ the standard
Fourier transform algorithm, the offset sampling Fourier transform
algorithm, or the non-uniform sampling Fourier transform
algorithm. Then, the one-dimensional equations are solved using the one-dimensional
finite element method with quadratic interpolation shape functions, resulting in three pentagonal equation sets, which are then solved using the pursuit method. This invention offers fast computation speed and small
memory footprint, providing a new method for efficient and high-precision numerical
simulation of large-scale
magnetotellurics.