The present invention discloses a method for inverting the density anomaly of the lithospheric mantle in the orogenic belt adjacent to the craton. The
gravity anomaly caused by the stripping of the crust and the Moho surface is removed by the
empirical relationship between crustal velocity and density and the gravity forward modeling method. The
gravity anomaly forward test of the lithospheric mantle plate in the orogenic belt adjacent to the craton is jointly constrained by the residual
gravity anomaly of the lithospheric mantle and the accumulated
mass anomaly of the residual lithospheric mantle, and a series of qualified model sets are obtained. Finally, the depth, thickness and residual density of the optimal plate model are selected by the interface constraint obtained by the
receiver function. The method of the present invention makes maximum use of known geophysical information and adopts a stable
estimation strategy to obtain reliable lithospheric mantle density anomaly information of the orogenic belt adjacent to the craton, thereby improving the research level of the lithospheric mantle density anomaly in the orogenic belt adjacent to the craton, and playing an important role in promoting the study of continental
lithosphere evolution and mineral resource exploration.