Temperature distribution is the key factor affecting the bonding quality in the glass/glass laser bonding process. In this work, the finite element method was used to establish three-dimensional (3D) numerical analysis model of the temperature field during bonding. Based on the result of the finite element analysis, the crucial parameters and their influences on the temperature distribution were discussed. In order to predetermine the necessary process parameter values for bonding, a nonlinear multiparameter fitting formula was established to predict the maximum temperature. The fitting model was validated experimentally by recording the maximum temperature during laser bonding via an infrared thermal imager.