The invention discloses a magnetotelluric three-dimensional forward modeling method based on a
spherical coordinate system, and belongs to the technical field of geophysical exploration. The method comprises the following steps: a, establishing a magnetotelluric
control equation; b, dividing the
spherical coordinate system into a plurality of small inverted quadrangular
prism grid units; c, setting spherical coordinate
model parameters, and constructing a spherical coordinate staggered grid unit with a discrete
magnetic field; d,
cycling a single frequency, and performing numerical dispersionon the magnetotelluric
control equation in the spherical coordinate staggered grid unit; e, solving a
system of linear equations; and f, calculating impedance via a spherical coordinate
tensor impedance formula, and then substituting the impedance into a Kania
apparent resistivity calculation formula to figure out the
apparent resistivity and phase at a measuring point. By adoption of the magnetotelluric three-dimensional forward modeling method disclosed by the invention, the interference of the earth curvature on the magnetotelluric three-dimensional deep detection can be effectively overcome, calculation errors caused by the earth curvature are avoided, the magnetotelluric three-dimensional forward modeling method is suitable for the numerical
simulation of magnetotelluric three-dimensional forward modeling, and can be better matched with an
earth model, and the accuracy is high.