Abstract: By using 919 P arrivals of 45 earthquakes in Shandong province and its neighboring area recorded in the Shandong seismic network, the 1D velocity model of Shandong area is obtained with the Kissling method.The earthquakes are relocated through simultaneous inversion combining epicenter and 1D P-wave velocity structure.After relocation, the travel time residual is reduced and the earthquake location accuracy of is improved.