Gravity prospecting profiles or dense gravity arrays are arranged in the earthquake occurrence area to obtain Bouguer gravity anomaly and invert the 2D/3D crust in this area. The physical property distribution difference and tectonic activity of main tectonic units are analyzed to provide a reference for the study of the seismogenic tectonic environment. The data are observed through high-precision gravity data and three-dimensional coordinates of the gravity station, respectively, by relative gravity survey under the control of absolute gravity and GNSS positioning survey.