Abstract: Based on the stations of Zhejiang digital seismic network and using 0.1°×0.1° grid points in Zhejiang as the virtual epicenters, we calculated the travel time of Pn, Sn, Pg and Sg seismic phases by 2015yr crust velocity model of Zhejiang province. After adding the random error on the travel time, we used Hyposat software to calculate the hypocenter location by velocity models with different range of disturbances. Comparing the virtual hypocenter with the calculation, we found that when the P-wave velocity deviation reaches within 0.1km/s, the deviation of epicenter is less than 3km. However, even in this case, there are also some improper location results in part of numerical simulations. Comparing with only making use of the Pg and Sg (Pb, Sb) phases, on the condition that Pn, Sn, Pg and Sg (Pb,Sb) phases are simultaneously used, the calculation of focal depth is closer to the virtual given value, and the calculation results are better in consistency.