In mathematics, statistics, and computer science, particularly in the fields of machine learning and inverse problems, regularization is a process of introducing additional information in order to solve an ill-posed problem or to prevent overfitting. The Tikhonov regularization method is widely used to solve complex problems in engineering. The vertical derivative of gravity can highlight the local anomalies and separate the horizontal superimposed abnormal bodies. The higher the order of the vertical derivative is, the stronger the resolution is. However, it is generally considered that the conversion of the high-order vertical derivative is unstable. In this paper, based on Tikhonov regularization for solving the high-order vertical derivatives of gravity field and combining with the iterative method for successive approximation, the Tikhonov regularization method for calculating the vertical high-order derivative in gravity field is proposed. The recurrence formula of Tikhonov regularization iterative method is obtained. Through the analysis of the filtering characteristics of this method, the high-order vertical derivative of gravity field calculated by this method is stable. Model tests and practical data processing also show that the method is of important theoretical significance and practical value.