A joint inversion method is adopted to invert 1D P-wave velocity of Jiangsu Province from bulletin data. Different from traditional method only for the first arrival,we use the first and secondary arrivals simultaneously to invert 1D velocity model. The secondary arrival for inversion can constrain velocity resolution of middle-lower crust. In order to constrain the bulletin data,a procedure is constructed to discriminate the different phases,such as Pn and Pg phase. It shows that this procedure can operate effectively in joint inversion. Synthetic travel-times are calculated with several 1D P-wave velocity models for bulletin data. The result presents that our inversion model performs the least RMS residuals. When the bulletin data have high location ability,the result of our inversion can be used as the initial model for 3D seismic tomography and routine earthquake determination in Jiangsu Province.