A method may include obtaining a P-wave velocity model and velocity ratio data regarding a geological region of interest. The method may further include generating, based on the P-wave velocity model and the velocity ratio data, an initial S-wave velocity model regarding the geological region of interest. The method may further include determining various velocity boundaries within the initial S-wave velocity model using a trained model. The method may further include updating the initial S-wave velocity model using the velocity boundaries, an automatically-selected cross-correlation lag value based on various seismic migration gathers, and a migration-velocity analysis to produce an updated S-wave velocity model. The method further includes generating a combined velocity model for the geological region of interest using the updated S-wave velocity model and the P-wave velocity model.