The invention discloses a structural constraint-based normalized gravity-magnetic-electric-seismic joint inversion method. The method includes the following steps that: S1, a normalized joint inversion objective function is established; S2, an initial model is established; S3, forward calculation and Jacobian matrix solving are carried out; S4, a normalized cross-gradient partial derivative and Lagrangian operator are calculated; S5, an iterative model is calculated; and S6, forward calculation is performed on the obtained iterative model, and the fitting difference of calculation results andobservation data is solved, if the fitting difference increases, the maximum number of iterations is reached, or the fitting difference reaches a threshold, loop inversion is stopped, and a normalizedjoint inversion model is outputted, otherwise S3 is executed. With the structural constraint-based normalized gravity-magnetic-electric-seismic joint inversion method of the invention adopted, the limitations of a single method are eliminated, the non-uniqueness of the inversion of geophysical inversion is decreased; and wrong inversion results caused by inversion based on the direct coupling ofphysical property parameters with large differences in order of magnitudes can be avoided. The method does not rely on the physical property relations of rock.