Abstract：In this paper, an approach for horizon matching with 3D multi-wave data is proposed. Simulated annealing algorithm is used to determine horizons on P- and PS-wave data by adjusting wavelet, time-depth and correlation coefficient in windows with target horizons on multi-offset and multi-wave seismic synthetics. By adjusting the P-wave and S-wave velocity ratio, converted-wave data can be compressed in the P-wave time domain. Therefore the comparison between P-wave and PS-wave in the same domain can be achieved. Results obtained by this method are accurate and objective with very low human intervention. Application results in a study area of Daqing Oilfield show that horizon calibration accuracy with P- and PS-waves is very good and multi-well calibration consistency are improved.