A simple and accurate evaluation method of broken rock zone thickness (BRZT), which is usually used to describe the broken rock zone (BRZ), is meaningful, due to its ability to provide a reference for the roadway stability evaluation and support design. To create a relationship between various geological variables and the broken rock zone thickness (BRZT), the multiple linear regression (MLR), artificial neural network (ANN), Gaussian process (GP) and particle swarm optimization algorithm (PSO)-GP method were utilized, and the corresponding intelligence models were developed based on the database collected from various mines in China. Four variables including embedding depth (ED), drift span (DS), surrounding rock mass strength (RMS) and joint index (JI) were selected to train the intelligence model, while broken rock zone thickness (BRZT) is chosen as the output variable, and the k-fold cross-validation method was applied in the training process. After training, three validation metrics including variance account for (VAF), determination coefficient (R2) and root mean squared error (RMSE) were applied to describe the predictive performance of these developed models. After comparing performance based on a ranking method, the obtained results show that the PSO-GP model provides the best predictive performance in estimating broken rock zone thickness (BRZT). In addition, the sensitive effect of collected variables on broken rock zone thickness (BRZT) can be listed as JI, ED, DS and RMS, and JI was found to be the most sensitive factor.