Seismic response assessment requires reliable information about subsurface conditions, including soil shear wave velocity (Vs). To properly assess seismic response, engineers need accurate information about Vs, an essential parameter for evaluating the propagation of seismic waves. However, measuring Vs is generally challenging due to the complex and time-consuming nature of field and laboratory tests. This study aims to predict Vs using machine learning (ML) algorithms from cone penetration test (CPT) data. The study utilized four ML algorithms, namely Random Forests (RFs), Support Vector Machine (SVM), Decision Trees (DT), and eXtreme Gradient Boosting (XGBoost), to predict Vs. These ML models were trained on 70% of the datasets, while their efficiency and generalization ability were assessed on the remaining 30%. The hyperparameters for each ML model were fine-tuned through Bayesian optimization with k-fold cross-validation techniques. The performance of each ML model was evaluated using eight different metrics, including root mean squared error (RMSE), mean absolute error (MAE), mean absolute percentage error (MAPE), coefficient of determination (R2), performance index (PI), scatter index (SI), A10−I, and U95. The results demonstrated that the RF model consistently performed well across all metrics. It achieved high accuracy and the lowest level of errors, indicating superior accuracy and precision in predicting Vs. The SVM and XGBoost models also exhibited strong performance, with slightly higher error metrics compared with the RF model. However, the DT model performed poorly, with higher error rates and uncertainty in predicting Vs. Based on these results, we can conclude that the RF model is highly effective at accurately predicting Vs using CPT data with minimal input features.