Abstract: Multi-wave seismic exploration offers distinct advantages in lithology identification, reservoir description, and gas-bearing property prediction. However, due to the different propagation mechanisms of longitudinal waves and transverse waves, the same geological interface often exhibits different kinematic and dynamic characteristics in the data of longitudinal waves (PP waves) and converted transverse waves (PS waves), which makes it difficult to achieve high-precision matching between PP waves and PS waves, especially in the northwest Sichuan region, especially when the signal inside the tight sandstone of the Shaximiao Formation is weak, the matching difficulty is even greater. By constructing pseudo shear wave seismic data similar to PS waves, the PP wave and PS wave data were matched on a large scale, and dynamic sliding time windows were used for local automatic fine matching, achieving precise matching of longitudinal and transverse waves in the target layer. On this basis, using converted shear wave data to calculate relative shear wave impedance, an initial three-dimensional lithofacies body is constructed and used as a constraint to perform phase controlled geostatistical inversion on the converted shear waves, forming a high-precision three-dimensional lithofacies body. Furthermore, constrained by high-precision three-dimensional lithofacies bodies, phase controlled geostatistical inversion is carried out on longitudinal wave data to obtain phase controlled longitudinal wave impedance inversion results. Based on the inversion results of phase controlled transverse wave impedance, the ratio of longitudinal and transverse wave velocities is calculated to achieve effective prediction of reservoir and gas bearing properties. This step-by-step phased multi wave joint prediction method improves the accuracy of reservoir and gas content prediction. Drilling results show that the predicted values of reservoir thickness and gas content have small errors compared to the measured values, proving the effectiveness of the method.